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Automatic  Detection  and  Recognition 
of  the  First  Arrival  Phase  of  Seismic 
Event  Signals  Contaminated  by  Noise. 


Kathleen  Ann  Alden 


ABSTRACT 


A  historical  prospective  on  the  methods  that  have  been  used 
to  automatically  detect  event  signals  and  pick  first  arrival 
times  is  developed.  Then  the  data  set  used  to  characterize  and 
test  the  detection  techniques  is  described.  Next  the  features 
which  discriminate  seismic  signals  effectively  from  the  back¬ 
ground  noise  are  characterized. 

Two  automatic  detection  methods  are  investigated.  (1)  The 
AR  (8) -spectral  estimates  of  the  signal  and  noise  are  used  to 
develop  the  AR  ( p ) -spectral  estimate  for  a  synthetic  waveform. 
(2)  Several  non-parame trie  tests  are  employed  in  a  time  domain 
detector  to  discriminate  event  signals  from  background. 

The  non-parametric  detector  chosen,  RANK  2700,  employs  a 
modified  rank  sum  test  to  locate  the  seismic  event  and  pick  its 
first  arrival  time.  Errors  in  the  automatic  first  arrival  pick 
for  152  of  the  event  traces  in  the  data  set  are  used  to  analyze 
the  performance  of  the  RANK  2700  detector. 


INTRODUCTION 


Since  the  advent  of  the  computer,  work  has  been  done 
in  many  fields  to  develop  techniques  to  harness  the 
computer  to  do  large  tedious  jobs  more  swiftly  and  with 
fewer  fluctuations  in  performance  than  its  human 
counterpart.  One  of  these  fields  is  the  identification 
and  classification  of  seismic  signals  from  digital  seismic 
event  records. 

There  are  thousands  of  seismic  events  occurring  every 
day  and  hundreds  of  seismic  stations.  It  is  not  practical 
to  evaluate  all  these  seismic  records  at  the  speed 
necessary  to  keep  up  with  the  influx  of  data.  Thus,  many 
studies  have  been  conducted  to  implement  various 
techniques  to  (1)  determine  first  arrival  times,  (2) 
classify  the  seismic  events  correctly  and  (3)  locate  the 
origins  of  the  seismic  events. 

It  is  becoming  increasingly  uneconomical  to  pick 
first  arrivals  by  hand  and  a  computer  can  be  used  to 
identify  first  arrivals  more  consistently  than  would  be 
picked  by  hand  on  an  oscillogram.  At  the  same  time, 
digital  recorders  are  becoming  more  common  on  even  low- 
cost  seismic  systems,  and  it  can  be  expected  that  in  the 
future  computer  techniques  will  become  more  attractive. 


Historical  Prospective  on  Methods  Used  To 
Automatically  Detect  Events  and 
Pick  First  Arrival  Times 

Many  scientists,  including  R.V.  Allen  (1978),  K.R. 
Anderson  (1978),  R.  Blandford  (1983),  K.S.Fu  (1982),  J.E. 
Gaby  (1983),  and  H.H.  Liu  (1981)  have  tried  to  devise  an 
effective  automatic  pattern  recognition  system  for  seismic 
signals  during  the  last  twenty  years.  Computer  techniques 
for  picking  first  events  have  yet  to  gain  widespread 
acceptance  (P.J.  Hatherly,  1982).  Effective  results  have 
been  elusive  due  to  the  nature  of  the  seismic  signals  and 
the  methods  used  to  model  and  predict  the  observed  values. 

Both  statistical  methods;  1 . e ., maximum  likelihood 
estimator,  maximum  entropy  spectra,  etc.  (C.H.  Chen, 
1981),  and  structural  methods;  i.e.,  pattern  recognition 
schemes  that  use  shape  features  such  as  slope,  radius  of 
curvative,  period  and  amplitude,  (J.E.  Gaby  and  K. 
Anderson,  1983),  have  been  employed  to  characterize  and 
identify  seismic  events. 

Statistical  classification  algorithms  can  be  grouped 
into  one  of  two  types,  parametric  or  non-parametr ic . 
Parametric  algorithms  assume  a  particular  class 
statistical  distribution,  commonly  the  normal 

distribution,  and  then  estimate  the  parameters  of  that 
distribution,  such  as  the  mean  and  variance,  to  use  in 
algorithm  classification.  Non-parametric  algorithms  make 
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no  assumptions  about  the  class  distributions.  Non- 
parametric  techniques  are  sometimes  termed  robust  because 
they  work  well  for  a  wide  variety  of  class  distributions, 
if  the  class  signatures  (mean,  variance,  etc.)  are 
reasonably  distinct.  Parametric  techniques  usually  yield 
good  results  under  the  same  conditions  as  the  non- 
parametric  techniques  if  the  signatures  of  the  classes  are 
reasonably  distinct,  even  if  the  assumed  class 
distribution  is  invalid. 

The  most  effective  example  of  the  statistical  methods 
is  94515  correct  recognition  of  regional  and  teleseismic 
events  using  maximum  entropy  spectra  and  spectral  ratio 
( energy>0 . 5Hz/energy<0 . 5Hz ) ,  (C.H.  Chen,  1982).  But  even 
94515  is  not  a  good  enough  performance  record  for  an 
automatic  pattern  recognition  system.  A  system  cannot  be 
considered  automatic  if  it  requires  human  watchdogs  to 
monitor  its  progress  and  correct  mistakes. 

Part  of  the  problem  of  correctly  characterizing  and 
classifying  seismic  signals  is  in  obtaining  an  accurate 
first  arrival  time.  There  is  always  ambiguity  associated 
with  measuring  the  first  arrival  time  from  seismograms, 
whether  it  is  done  by  seismologist  or  machine,  since  these 
signals  of  finite  bandwidth  are  of  unknown  shape  and 
contaminated  by  noise.  Such  ambiguity  can  be  reduced  by 
combining  the  processes  of  picking  arrivals  in  an 
iterative  fashion  (C.H.  Chen,  1982).  Inaccurate  first 
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arrivals  distort  the  structural  and  statistical 
characterizations  of  a  seismic  event  that  are  used  for 
classification. 

Effective  results  have  been  elusive  due  to  the  nature 
of  the  seismic  signal  which  contains  events  contaminated 
by  noise.  In  many  fields  of  signal  processing:  for 
example,  the  development  of  speech  recognition  systems, 
the  automatic  analysis  of  signals  requires  the  recognition 
of  specific  features  in  the  signal.  (Gaby  and  Anderson, 
1983)  Since  more  a  priori  knowledge  exists  about  the 
characteristic  shapes  of  words,  it  is  easier  to 
automatically  Identify  a  word  and  classify  it  correctly. 
The  paucity  of  a  priori  information  available  on  the 
morphologies  (shapes)  of  seismic  signals  inhibits 
automatic  pattern  recognition  (Gaby  and  Anderson,  1983). 

The  lack  of  a  priori  knowledge  characterizing  seismic 
signals  has  prompted  the  recent  application  of  pattern 
recognition  techniques  to  find  and  develop  methods  to 
discriminate  and  classify  seismic  signals.  The  methods 
used  to  find  the  broad  characteristics  of  a  seismic  event 
include  storing  the  trace  in  a  binary  tree  structure  and 
using  affinity  techniques,  first  developed  for  image 
processing,  to  combine  small  segments  of  the  signal  and 
store  the  signal  within  the  tree  structure  at  different 
levels  of  complexity  (Gaby  and  Anderson,  1983).  In  other 
words,  this  method  would  use  features  such  as  period, 


slope  and  amplitude  to  divide  the  event  trace  into 
segments  and  then  try  to  associate  those  segments  with  the 
sequence  of  arrival  phases  that  make  up  a  seismic  event. 

Augmented  transition  networks  ( ATN ) ,  originally 
created  to  provide  a  formal  environment  to  develop 
grammatical  rules  describing  finite  state  grammers,  have 
also  been  used  to  develop  seismic  signal  structural 
characteristics  interactively  (Anderson,  1981). 

All  the  techniqes  for  seismic  pattern  recognition 
schemes  require  that  the  first  arrival  time  of  the 
waveform  be  known.  The  seismologist  uses  the  morphology 
(shape)  of  the  seismic  signal  to  identify  the  first 
arrival  time  correctly  along  with  changes  in  period  and 
amplitude.  Determining  the  correct  first  arrival  time  is 
important  in  classifying  the  signal  and  locating  the 
origin  of  the  seismic  event.  For  an  automatic  pattern 
recognition  system  to  replace  the  seismologist  it  must 
identify  first  arrival  times  correctly  and  classify  the 
event  with  better  than  the  9485  current  success  rate  of 
techniques  previously  attempted. 

Automatic  processes  that  have  been  employed  to 
determine  first  arrival  times  rely  on  statistics  to 
distinguish  between  two  populations  (signal  and  noise) 
occurring  on  a  seismic  trace.  Since  the  density 
distribution  functions  for  both  the  signal  and  noise  are 
not  known  a  priori  they  have  traditionally  been  estimated 
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background  noise.  Robinson( 1967 )  wrote  subroutines  to 
pass  the  proposed  first  arrival  if  there  was  a  significant 
prediction  error  at  the  first  arrival  and  for  two  terms 
after  it.  Significance  was  established  by  making 
comparisions  with  the  prediction  errors  within  the  noise. 
A  test,  which  predicts  the  values  at  the  first  arrival 
from  the  previous  values  using  the  technique  of  linear 
least-square  prediction,  was  first  used  by  Wadsworth  et 
al ( 1953 )  to  identify  seismic  events  on  the  basis  of  a 
prediction  failure.  (P.J.  Hatherly,  1982) 

Most  techniques  to  distinguish  noise  from  the  first 


arrival  signal  use  parametric  statistics  which  assume  the 
signal  has  a  known  distribution.  This  is  an  incorrect 
assumption  for  the  background  noise  preceding  a  seismic 
event.  In  general,  background  noise  is  a  non-stationary 
process  due  to  seasonal  changes  and  atmospheric 

varic  :ions,  with  an  unknown  distribution.  Over  a  short 
time  interval,  less  than  100  seconds,  background  noise  can 
usually  be  considered  stationary  except  in  the  case  where 
it  precedes  a  seismic  event.  When  a  seismic  event  occurs 
the  mean  of  the  signal  often  flucuates  while  the  transient 
signal  is  being  recorded  by  the  seismic  Instrument. 

Non-stationary  processes  have  time  changing  means 
and/or  variances.  Stationary  process  implies  that  the 
mean,  the  variance  and  the  autocovariances  of  the  process 
are  invariant  under  time  translations.  Thus  the  mean  and 
variance  are  constant,  and  the  autocovariances  depend  only 
on  the  lag  time. 

Non-stationary  time  series  have  been  modeled  by 
several  processes:  (1)  Harrison( 1964 )  used  an 

exponentially  weighted  moving  average,  EWMA,  to  forecast 
seasonal  short  term  sales.  (2)  A  modified  autogressive 
moving  average,  ARMA,  model  with  time  varying  coefficients 
of  the  form: 

Cjt*xt_j  +  dt  ,  r  «  max(p-l,q-m) 
has  been  used  by  P.  Whittle  (1965),  to  model  non- 


stationary  time  series.  Prediction  formulas  applied  in 
the  stationary  case  are  shown  by  Niemi(1983)  to  also  be 
valid  in  the  non-stationary  case.  The  effects  of  the  non- 
stationarity  on  the  estimation  of  the  parameters  of  the 
underlying  ARMA  model  are  not  shown  to  be  significant.  (3) 
The  autoregressive  integrated  moving  average,  ARIMA. 
method,  based  on  Gaussian  stochastic  processes  was 
developed  by  Box  and  Jenkins  (1970)  to  model  homogeneous 
non-stationary  time  series. 

Homogeneous  non-stationarity  implies  the  changing 
mean  can  be  described  by  a  low  order  polynomial  in  time. 
However,  the  coefficients  of  the  polynomial  are  not 
constant  but  vary  with  time.  The  observations  are 
described  by  random  stochastic  trends  (polynomials ) . 
Tintner  (1940),  Yaglora  (1955),  and  Box  and  Jenkins  (1976) 
argue  that  homogeneous  non-stationary  sequences  can  be 
transformed  into  stationary  sequences  by  taking  successive 
differences  of  the  series. 

In  practice  It  is  usually  the  first  or  second 
Integral  of  a  non-stationary  process  which  is  stationary. 
The  ARIMA  technique  integrates  the  non-stationary  time 
series  until  it  is  stationary  and  then  models  the 
resulting  time  series  as  an  ARMA  process.  Under  fairly 
general  conditions  the  prediction  interval  for  a  future 
observation  in  an  ARIMA  scheme  is  robust  with  respect  to 
symmetric  non-normality  of  the  error  distribution, 
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Non-parametric  methods  that  make  no  assumptions  about 
the  density  distribution  of  the  time  series  are  outlined 
by  Kassam  and  Thomas ( 1975 ) .  Cobben(l982)  outlines  a  non- 
parametric  detector  based  on  a  sign  test  that 
discriminates  between: 

H0:  x(t)  =  n(t)  //data  =  noise// 

Ha:  x(t)  =  s(t)  +  n(t)  //data  =  signal  +  noise// 

included  in  Kassam  and  Thomas  (1980).  The  sign  test  was 
extended  by  Cox  (1955)  to  detect  steps  and  ramps  in  the 
presence  of  additive  noise;  i.e.,  signals  with  sharp  onset 
and  gradual  onset  first  arrivals.  Cobben's  method  assumes 
a  stationary  signal  which  may  be  a  valid  assumption  for  a 
time  window  of  less  than  100  seconds,  (C.H.  Chen,  1978). 

Although  background  noise  may  possibly  be  considered 
stationary,  it  cannot  be  assummed  to  be  Gaussian.  Various 
known  and  unknown  processes:  microseisms  which  are  "non- 
Gaussian";  thermal  noise  due  to  current  across  resistors 
which  is  "Gaussian";  and  seismometer  noise  which  is 
"distributed  as  1/f"  contribute  to  form  a  non-Gaussian 
density  distribution.  Although,  Heuts  (1981),  has  shown 
that  an  ARIMA  prediction  model  (based  on  an  assumption  of 
Gaussian  and  stochastic  data)  is  robust  with  respect  to 
symmetric  non-normality  of  the  error  distribution, 
parametric  techniques  are  generally  incorrect  for 
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distinguishing  noise  from  the  first  arrival. 

Part  of  the  difficulty  in  automatically  classifying 
seismic  signals  and  locating  their  origin  is  a  direct 
result  of  Incorrect  first  arrival  times  picked  by  methods 
assuming  an  underlying  Gaussian  distribution  for  the 
signal.  Since  the  seismic  signals  have  an  unknown 
underlying  distribution  a  method  for  first  arrival 
detection  using  non-parametr ic  statistics  should  be 
developed  in  the  hopes  that  it  will  improve  the  accuracy 
of  the  first  arrival  time. 

Once  the  first  arrival  time  can  be  swiftly  and 
accurately  located  the  morphological  and  statistical 
methods  of  classifying  seismic  signals  can  be  used  with 
greater  accuracy  since  tht  correct  sequence  of  events  (or 
states)  that  characterize  the  seismic  signal  will  be 
readily  available.  The  first  arrival  on  the  seismic  trace 
is  identified  as  the  first  observation  which  is 
statistically  different  from  the  observation  before.  The 
point  chosen  is  dependent  on  the  signal  to  noise  (S/N) 
ratio  and  the  amplification  of  the  seismic  signal  (P.J. 
Hatherly ,  1982 )  . 

Some  signal  detection  schemes  used  to  record  events 
in  real  time  use  the  spectrum  found  using  the  Walsh 
(Goforth  and  Herrin,  1981)  or  Fourier  (Blandford,  1983) 
transforms  of  the  signal  to  distinguish  the  signal  from 
the  background  noise.  These  methods  while  effective  for 


real  time  detection  of  events  can  only  locate  the  first 
arrival  time  within  one  window  (a  block  of  observed  data 
values  transformed  to  provide  a  spectrum  for  comparision) . 
Figure  1  illustrates  the  fast  Walsh  transform. 

If  this  window  is  large  enough  to  detect  the  first 
arrival  (greater  than  two  times  the  longest  period  of  the 
expected  signal,  2.5  seconds  in.  our  case)  it  will  not  be 
small  enough  to  accurately  determine  first  arrival  times 
for  the  purpose  of  classifying  and  locating  the  origin  of 
events.  This  indicates  an  effective  automatic  first 
arrival  picker  should  be  implemented  in  the  time  domain. 

Methods  Attempted  in  this  Study 


Since  the  signal  is  nonstationary  with  an  unknown 
distribution  function  it  is  logical  to  develop  a  detector 
that  does  not  assume  stationar ity ;  i.e.,  constant  mean  and 
constant  variance,  or  a  "known"  distribution  function. 
There  is  a  broad  range  of  amplitudes  and  bandwidths  which 
characterize  local,  regional  and  teleseismic  earthquakes 
and  explosions.  Thus,  we  must  develop  features  to 
distinguish  a  broad  category  of  signals  from  the 
background  noise. 

Since  the  seismic  signal  we  record  has  an  unknown 
distribution,  it  is  important  to  determine  how  valid  some 
of  the  parametric  techniques;  i.e.,  techniques  that  assume 
a  know  distribution,  are  when  applied  to  a  seismic  signal. 
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Figure  1.  A  figurative  example  of  a  6.4  Second 
Detection  Window  and  its  Corresponding  Normalized  Walsh 
Spectrum  Plotted  on  Top  of  the  Normalized  Walsh  Spectrum 
for  the  Background  Noise. 


In  particular,  Niemi  (1983)  wrote  a  paper  indicating  an 
ARMA  scheme  was  robust  with  respect  to  symmetric  non¬ 
normality.  The  discrete  density  distributions  for  the 
signal  and  noise  were  plotted  for  100  seismic  events 
recorded  at  each  of  the  four  seismic  stations  in  far  west 
Texas  (Lajitas,  Marathon,  Shafter  and  Tres  Cuevas)  to  find 
out  if  the  density  distributions  of  the  seismic  signals 
and  noise  from  each  station  could  be  characterized  as 
symmetric  even  though  the  "true"  underlying  density 
distributions  were  unknown.  Histograms  of  the  amplitude 
of  the  seismic  signal  were  computed  to  approximate  the 
discrete  density  functions  for  the  noise  and  signal  plus 
noise  of  each  event  trace.  (See  figures  2  through  4)  In 
all  cases  the  discrete  density  functions  for  both  the 
signal  and  noise  were  symmetric.  This  assumption  allows 
us  to  model  the  non-stationary  seismic  signal  with  a 
simplier  model  than  the  ARIKA  model  for  which  methods  of 
computation  have  been  more  fully  explored. 

Two  methods  were  attempted  to  develop  an  automatic 
first  arrival  picker.  The  first  technique  uses  the  AR- 
spectral  estimation  of  the  signal  and  the  background  noise 
to  develop  a  synthetic  waveform  to  cross-correlate  with 
the  event  trace  and  pick  out  the  first  arrival  phase.  The 
second  technique  employs  a  non-parametr ic  test  within  a 
sliding  window  detector  to  identify  the  seismic  event 
signal  and  pick  the  first  arrival  phase. 
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Figure  3.  The  First  60  Seconds  for  a  Regional 
Event,  Recorded  at  Each  of  the  Four  Seismic  Stations, 
Lajitas,  Marathon,  Shafter  and  Tres  Cuevas,  are  Plotted 
from  Bottom  to  Top  Respectively.  The  Corresponding 
Discrete  Density  Functions  for  the  First  20  Seconds  of 
Background  Noise  and  the  Next  20  Seconds  of  Signal  Plus 
Noise  for  Each  Trace  are  Shovm  from  Left  to  Right. 
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Figure  4.  The  First  60  Seconds  for  a  Teleseismic 
Event,  Recorded  at  Each  of  the  Four  Seismic  Stations, 
Lajitas,  Marathon,  Shatter  and  Tres  Cuevas,  are  Plotted 
from  Bottom  to  Top  Respectively.  The  Corresponding 
Discrete  Density  Functions  for  the  First  20  Seconds  of 
Background  Noise  and  the  Next  20  Seconds  of  Signal  Plus 
Noise  for  Each  Trace  are  Shown  from  Left  to  Right. 
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DESCRIPTION  OF  THE  DATA  SET 

Instruments  Used  to  Record  the  Data 

Seismic  events  recorded  at  four  stations  in  southwest 
Texas,  Lajitas,  Maratnon,  Shafter  and  Tres  Cuevas,  were 
collected  during  the  Group  of  Scientific  Experts  Technical 
Test  (GSETT)  in  which  Southern  Methodist  University  was  a 
participant.  GSETT  was  conducted  from  October  15,  1984  to 
December  14,  1984  to  test  procedures  for  exchanging 

seismic  data,  which  include:  the  extraction  of  Level  I 

signal  parameters  (seismic  phase  identifiers,  arrival 
times,  signal  amplitudes  and  periods);  the  exchange  of 
these  parameters  primarily  via  the  Global 

Telecommunications  System  of  the  World  Meteorological 
Organization  (WMO/GTS);  and  the  collection  and  assesment 
of  these  data  at  the  Data  Center  in  Washington,  D.C. 

Table  1  gives  the  locations  and  descriptions  of  the 
instruments  used  to  record  the  seismic  events  during  the 
test.  The  instrument  responses  are  illustrated  in  figures 
5  and  6 . 

A  real  time  event  detector  utilizing  the  fast  Walsh 
transform  (Goforth  and  Herrin,  1981)  implemented  on  a  DEC 
RT/11  micro-computer  received  the  seismic  data  from  the 
Lajitas,  Marathon,  Shafter  and  Tres  Cuevas  stations  via 
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TABLE  1.  Locations  and  Descriptions 
of  Instruments  Used  to  Record 
the  Event  Data 


SEISMOMETER  DESCRIPTIONS 


NAME 

LRJITAS 

MARATHON 

SHATTER  TRES 

CUEVAS 

ID 

LA 

MA 

SH 

TR 

OWJEL  NUMBER 

2 

'3 

4 

l 

INSTRUMENT  TVPE 

23900 

23900 

23900 

18300 

LATITUDE  DEO  N 

23.333 

30.306 

23.924 

29.316 

LONGITUDE  DEG  H 

103.667 

103. 2SS 

104.371 

103.717 

ELEVATION  METERS 

1013 

1356 

1570 

1097 

VCO  HZ 

960 

1300 

1700 

2040 

LPF  CORNER.  HZ 

5 

5 

5 
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telephone  line.  (See  appendix  A  for  a  description  of  the 
fast  Walsh  transform  detector.)  The  detected  events  were 
recorded  on  magnetic  tape  for  later  graphic  display  and 
analysis  by  a  seismologist. 

The  seismic  analyst  was  responsible  for  identifying 
phase  (Pn,Pg,P  or  Lg)  and  picking  first  arrival  times, 
amplitudes  and  periods.  The  bias  of  the  analyst  directly 
effects  the  choice  of  local  events  v.s.  "noise". 
Discrimination  of  the  phase  of  the  first  arrival  (P  or  Lg ) 
is  dependent  on  which  phase  the  analyst's  experience 
indicates  the  waveform  matches.  Parameters  for  the  events 
chosen  by  the  analyst  as  "true"  events.  not  "glitches"  or 
"background  noise",  were  sent  to  the  U.S.  National  Data 
center  in  Washington  using  the  UNIX-net  software  of  the 
Eunice  operating  system. 

GSETT's  Final  Event  List  associated  all  the  event 
information  collected  from  throughout  the  world  and 
reported  the  origin  time,  location,  magnitude  and  number 
of  stations  used  to  define  each  event  which  was  detected 
during  the  two  month  time  period  of  the  test. 

The  residual  travel  time  (the  difference  between  the 
arrival  time  at  the  station  Indicated  by  the  best  least 
squares  fit  of  the  stations  associated  with  an  event  and 
the  first  arrival  time  picked  by  the  analyst)  was  listed 
for  each  station  used  to  define  an  event.  The  "true" 
estimates  calculated  by  least  squares,  asstiming  the 
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correct  location  of  the  event  was  found,  allow  us  to 
determine  if  the  automatic  detection  method  is  doing  as 
well,  worse,  or  better  than  the  human  analyst  in  ambiquous 
cases . 

Methods  Used  to  Detect  and  Record 
Events  in  the  Event  Files 

The  fast  Walsh  transform  detection  window  is  6.4 
seconds  long  with  a  50*  overlap.  (see  figure  1)  When  the 
spectra  of  the  detection  window  exceeds  the  background 
noise  threshold  the  detection  is  triggered.  Thirty 
seconds  of  data  preceding  the  detection  window  is  written 
at  the  beginning  of  the  event  file  so  that  if  the  Walsh 
detector  triggered  on  an  Lg  arrival  the  P  wave  first 
arrival  should  be  present  somewhere  in  the  first  30 
seconds  of  the  event  file.  Thirty  seconds  was  originally 
thought  to  be  adequate  to  make  sure  the  P  arrival  was 
included  in  the  event  trace  but  experience  in  working  with 
the  data  collected  during  the  GSETT  experiment  has  shown 
that  60  seconds  is  a  better  predetection  interval  for 
insuring  the  first  arrival  (Pn)  is  included  in  the  event 
file  when  the  detector  triggers  on  Pg. 


EVENT  SIGNAL  CLASSIFICATION 


Variation  in  the  Signals  to  be  Distinguished 
From  the  Background  Noise 


The  collection  of  seismic  events  we  are  attempting  to 
distinguish  from  the  background  noise,  (recorded  on  short 
period  seismic  instruments  with  a  bandwidth  of  0-20  Hz), 
encompasses  a  broad  range  of  periods  and  amplitudes. 
Local  events,  occurring  within  25  to  150  kilometers  of  the 
seismic  station,  are  characterized  by  amplitudes  which  may 
be  nearly  the  same  as  the  noise  and  by  periods  of  0.3  to 
0.4  seconds.  Regional  events,  occuring  on  the  same 
continent  as  the  seismic  station,  are  characterized  by 
varying  amplitudes  and  periods  of  0.5  to  0.7  seconds. 
Teleseismic  events,  events  that  travel  distances  greater 
than  about  2000  km  to  reach  the  seismic  recording  station, 
are  characterized  by  periods  of  0.8  to  1.2  seconds.  The 
high  frequency  component  of  the  background  noise  is 
similar  in  spectral  content  to  local  seismic  event  signals 
while  the  underlying  low  frequency  component  of  the 
background  noise  is  similar  in  spectral  content  to 
teleseismic  signals.  Pn  can  often  not  be  discriminated 
from  the  background  noise  by  the  analyst,  thus  Lg  which 
follows  Pn  is  often  picked  as  the  first  arrival  phase  for 
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Figure  7 .  Examples  of  the  Various  Types  of  Seismic 
Signals  Typically  Encountered.  (a)  33  1cm  -  Local  Event  - 
SPZ  recorded  in  Berkeley,  California.  (b)  140  km  -  Local 
Event  —  Closest  distance  Pn  and  Pg  can  be  observed 
separately  -  SPZ  recorded  in  Colorado.  (c)  6°  -  Regional 
Earthquake  -  SPZ  recorded  in  Winner,  South  Dakota.  (d) 
13°  -  Regional  to  Teleseismic  Earthquake  -  Gulf  of 

California  -  SPZ  recorded  in  Colorado.  (e)  91.5° 

Teleseismic  Earthquake  -  Offshore  Peru  -  SPZ  recorded  in 
Black  Rapids,  Alaska. 
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regional  events.  Figure  7  illustrates  typical  examples  of 
the  types  of  seismic  signals  encountered.  (R.B.  Simon, 
19S 1 ) 


Appropriateness  of  the  Data  Set  for  Testing 
the  Effectiveness  of  Any  Automatic 
Techniques  Devised 


The  data  set  collected  during  GSETT  represents  a 
collection  of  event  records  of  different  types  of  events 
with  continuous  background  noise  variations  from  four 


different  seismic  stations. 

This  data 

set 

has  variations 

in:  signal  type 

( local , 

regional 

and 

teleseismic )  ; 

background  noise 

(due  to 

changes 

in 

wind  speed, 

temperature,  humidity,  barometric  pressure  and  cultural 
activity) ;  instrument  type  and  local  structure  at  the 
seismic  station.  It  is  concluded  that  because  the  data 
set  includes  variations  in  all  the  parameters  normally 
expected  to  change  when  recording  seismic  events  the  data 
set  will  be  a  good  test  of  the  robustness  of  any  detection 
method  used  to  pick  first  arrival  times. 

Indentify  the  Characteristics  of  the  Seismic 
Signals  Which  Define  a  Class  Distinctive 
From  the  Background  Noise 

Any  classification  technique  actually  consists  of  a 
one-time  calculation  of  decision  boundaries,  followed  by  a 
comparison  of  each  sample's  feature  vector  and  the 
location  of  those  boundaries.  By  incorporating  these 
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facts  in  classification  algorithms,  their  efficiency  can 
be  improved  by  an  order  of  magnitude  or  more,  with  little 
or  no  reduction  in  accuracy.  ( Schowengerdt ,  R.A.,  1983) 

Easily  calculatable  features  available  in  the  time 
domain  for  discriminating  seismic  signals  and  background 
noise  are:  period  estimates,  (intervals  between  zero- 
crossings  and  slope  direction  changes);  amplitude;  and 
slope.  Figure  8  is  a  typical  example  of  the  background 
noise.  This  noise  data  was  recorded  during  a  background 
noise  study  February  8,  1985  at  the  Lajitas  test  site.  In 
the  next  three  sections  we  determine  the  effectiveness  of 
these  features  in  distinguishing  a  broad  class  of  seismic 
signals  from  the  background  noise  by  using:  (1)  a 
threshold  detector  based  on  an  estimate  of  the  mean  and 
variance  of  the  background  noise;  (2)  segmentation  of  the 
signal  based  on  the  similarity  of  features  to  find  a 
pattern;  and  (3)  division  of  the  features  based  on  their 
ranks  into  rank  quartiles  to  look  for  a  relative  pattern 
in  the  feature  distribution. 


A  Threshold  Detector  to  Identify  Features 


For  the  threshold  detector  the  mean  and  variance  were 
calculated  for  the  amplitude  and  slope  (average  difference 
between  the  ith  sample  and  the  samples  adjacent  to  it)  of 
the  first  100  observations  representating  background  noise 


assuming  a  normal  distribution.  The  threshold  for  signal 
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Figu-e  8.  Background  Noise  Recorded  on  11  Different  Instruments 
During  a  Noise  Survey  at  the  Lajicas  Test  Site,  February  8,  1985.  In¬ 
strument  1  -  GS  13  Z  (vault) ,  Instrument  2  -  GS13N  (vault) ,  Instru¬ 
ment  3  -  GS13  E  (vault) ,  Instrument  4  -  18300  (mine  site) ,  Instru¬ 
ment  5  -  S750  Z  (grouted  in) ,  Instrument  6  -  S750  Z  (vault) ,  In¬ 
strument  7  -  S750  A  (North  of  mud  hut  150m) ,  Instrument  8  -  S750A 
(East  of  mud  hut  150m),  Instrument  9  -  GS21  (buried  330'),  In¬ 
strument  10  -  GS21  (buried  50'),  Instrument  11  -  microphones  (wind). 
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detection  was  set  at  two  standard  deviations  from  the 
estimated  mean  of  the  background  noise.  Thus,  given  that 
background  noise  is  all  that  is  present  in  the  seismic 
signal,  there  is  a  95%  probability  that  the  amplitude  and 
slope  for  any  observation  will  be  within  two  standard 
deviations  of  its  estimated  mean.  Figures  9  through  11 
illustrate  the  results  of  this  test  on  three  types  of 
seismic  events  recorded  at  the  four  stations  Lajitas, 
Marathon,  Shatter  and  Tres  Cuevas.  The  state  table  for 
the  color  changes  occuring  along  the  event  traces  in 
figures  9  through  11  is  shown  in  Table  2.  The  mean  and 
variance  estimated  for  the  slope  and  amplitude  were  not 
very  effective  in  clearly  distinguishing  the  background 
noise  from  the  seismic  event  signal. 

Segmentation  Using  Affinity  Techniques 
to  Identify  Features 

Segmentation  of  the  trace  using  the  estimated  mean 
and  variance  of  the  slope  and  amplitude  of  the  background 
noise  was  done  using  the  nearest  neighbor  decision  rule 
and  an  affinity  algorithm  to  combine  the  samples  with  the 
most  similar  features  into  a  finite  number  of  segments. 
The  nearest  neighbor  decision  rule  decides  which  of  a 
segments  neighbors  is  most  similar.  Then  the  affinity 
algorithm  combines  those  segments  which  mutually  consider 
each  other  most  similar.  Figures  12  through  14  illustrate 
this  technique  on  the  three  types  of  seismic  signals, 
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local,  regional  and  teleseismic.  The  colors  represent 
different  segments.  No  clear  pattern  was  identified  in 
the  segmentation  to  make  this  a  useful  method  for 
distinguishing  signal  from  noise. 

Both  the  threshold  detection  and  segmentation  methods 
made  the  assumption  that  the  noise  was  Gaussian  and 
stationary  when  the  mean  and  the. variance  were  estimated. 
Since  this  is  untrue,  the  noise  could  not  be  effectively 
characterized  by  the  estimated  mean  and  variance  and 
overlap  occured  between  the  boundaries  defining  the  signal 
and  noise  classes.  The  subroutine,  "NEIGHBOR",  listed  in 
appendix  B  uses  the  nearest  neighbor  rule  and  affinity 
techniques  to  combine  the  segments. 

Using  Rank  Quartiles  to  Identify  Features 

The  next  attempt  to  characterize  the  signal  and  noise 
categories  was  to  rank  the  amplitudes  and  periods, 
estimated  by  the  distance  between  zero  crossings,  of  the 
event  trace.  Then  the  ranks  were  divided  into  quartiles 
and  color  coded  on  the  event  trace  so  any  patterns  in  the 
features  associated  with  the  signal  or  noise  might  be 
identified.  Table  3  illustrates  the  color  code  used  in 
figures  15  through  19  to  identify  which  rank  quartile  an 
observation  belongs. 

This  technique  showed  a  clear  pattern  for  ranked 
amplitudes  of  events  with  large  S/N  ratios  but  was 
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TABLE  3 


State  Table  for  Color  Changes  Used 
to  Graphically  Depict  the 
Rank  Quartiles 
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Figure  15 *a 

Figure  15.  The  Ranks  of  Amplitudes  and  Periods  are 
Divided  into  Rank  Quartiles.  The  Color  Changes  in  the 
Event  Traces,  Recorded  at  the  Stations  Lajitas,  Marathon, 
Shafter  and  Tres  Cuevas,  Bottom  to  Top  Respectively, 
Graphically  Depict  any  Patterns  that  Exist  in  the 
Features.  (a)  The  Color  Index  Corresponding  to  the  16 
Rank  Quartiles  the  Ranks  of  the  Periods  and  Amplitudes  of 
Each  Trace  are  Divided  into  as  Shown  in  Table  3.  {  / 
Shows  the  Rank  Quartile  Pattern  for  a  Regional  Event. 
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Figure  17. b 
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Figure  17. a 

Figure  17.  The  Ranks  of  Amplitudes  and  Periods  are 
Divided  into  Rank  Quartiles.  The  Color  Changes  in  the 
Event  Traces,  Recorded  at  the  Stations  Lajitas,  Marathon, 
Shafter  and  Tres  Cuevas,  Bottom  to  Top  Respectively, 
Graphically  Depict  any  Patterns  that  Exist  in  the 
Features.  (a)  Illustrates  the  Rank  Quartile  Pattern  for  a 
Regional  Event.  (b)  Illustrates  the  Rank  Quartile  Pattern 
for  a  Regional  Event. 
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Figure  18. a 

Figure  18.  The  Ranks  of  Amplitudes  and  Periods  are 
Divided  into  Rank  Quartiles.  The  Color  Changes  in  the 
Event  Traces,  Recorded  at  the  Stations  Lajitas,  Marathon, 
Shatter  and  Tres  Cuevas,  Bottom  to  Top  Respectively, 
Graphically  Depict  any  Patterns  that  Exist  in  the 
Features.  (a)  Illustrates  the  Rank  Quartile  Pattern  for  a 
Regional  Event.  (b)  Illustrates  the  Rank  Quartile  Pattern 
for  a  Regional  Event. 
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Figure  19. a 

Figure  19.  The  Ranks  of  Amplitudes  and  Periods  are 
Divided  into  Rank  Quartiles.  The  Color  Changes  in  the 
Event  Traces,  Recorded  at  the  Stations  Lajitas,  Marathon, 
Shatter  and  Tres  Cuevas,  Bottom  to  Top  Respectively, 
Graphically  Depict  any  Patterns  that  Exist  in  the 
Features.  (a)  Illustrates  the  Rank  Quartile  Pattern  for  a 
Local  *  Event .  (b)  Illustrates  the  Rank  Quartile  Pattern 

for  a  Regional  Event. 
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ambiguous  for  signals  buried  in  noise.  The  rank  pattern 
for  estimated  periods  showed  no  strong  pattern  for 
distinguishing  signal  from  noise.  However,  the  ranks  for 
the  amplitude  and  estimated  periods,  which  were  added 
together  and  divided  into  quartiles,  was  the  the  best 
indicator  of  a  clear  pattern  in  the  features  to  be  used  to 
discriminate  signal  and  noise.  This  is  reasonable  since 
the  criteria  the  human  analyst  uses  to  distinguish  signal 
from  noise  especially  in  ambiguous  cases  with  low  signal 
to  noise  ratios  is  based  on  changes  in  both  period  and 
amplitude  from  the  preceding  background.  Slope  is 
sensitive  to  changes  in  both  amplitude  and  period  which 
makes  it  the  most  favorable  feature  to  use  in  creating  an 
automatic  first  arrival  detection  algorithm. 

Conclusions 

Since  parametric  techniques  appear  to  allow  too  much 
overlap  in  the  classifications  to  clearly  distinguish 
signal  from  noise,  a  non-pararaetric  technique  using  the 
slope  of  the  event  trace  appears  to  describe  the  detection 
method  with  the  best  chance  of  discriminating  background 
noise  from  a  broad  class  of  seismic  event  signals. 

The  feature  with  the  least  amount  of  overlap  between 
background  noise  and  seismic  event  signals,  "modified" 
slope,  was  calculated  through  discrete  integration  of  the 
event  trace  slope  over  segments  with  the  same  slope 
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direction.  Figure  20  demonstrates  this  calculation. 
Various  period  estimations  of  the  background  noise  were 
too  similar  to  one  or  more  of  the  local,  regional  or 
teleseismic  signal  classes  for  effective  discrimination 
between  signal  and  noise.  Amplitude  was  a  good 
discriminator  when  the  S/N  ratio  was  high.  However,  the 
detectors  failed  to  trigger  on  the  first  arrival  for 
events  with  low  S/N  ratios.  Both  the  slope  and  "modified" 
slope  of  the  event  trace  are  characterized  by  a  constant 
mean  equal  to  zero  and  a  variance  which  changes  with  time. 
This  makes  the  "modified"  slope  of  the  event  trace  a  more 
attractive  feature  for  distinguishing  between  signal  and 
noise  than  the  unmodified  event  trace  whose  mean  and 
variance  both  change  with  time. 
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DETECTORS 


AR-Soectral  Estimation  Used  to  Develop 
a  Synthetic  Waveform 

Assumptions 

The  autoregressive  (AR)  model  is  correct  for  modeling 
minimum  phase  stationary  processes.  To  estimate  the 
spectral  density  function  we  must  assume  that  we  can 
obtain  an  "adequate  approximation"  by  using  a  finite  order 
AR  model  (Priestley,  1981).  The  development  of  a 
synthetic  waveform  to  pick  out  the  first  arrival  phase  in 
an  event  trace  from  the  autoregressive  spectral  estimate 
of  that  same  event  trace  requires  several  assumptions. 
First,  in  order  to  use  this  technique,  we  must  assume  the 
first  arrival  phase  of  the  seismic  signal  generated  by  an 
earthquake  or  explosion  is  minimum  phase.  Second  we 
assume  the  corner  frequency,  or  dominate  period,  of  the 
seismic  signal  has  the  same  period  as  the  first  arrival 
phase.  Third  we  assume  that  the  background  noise  does  not 
contain  an  obvious  corner  frequency.  Forth  we  assume  the 
seismic  signal  is  stationary  or  at  least  can  be  estimated 
by  a  process  that  assumes  the  driving  function  is  white 
noise.  Seismic  signals  contaminated  by  noise  are 
nonstationary  processes  with  unknown  distributions.  Hannu 


Niemi,  (1983),  studied  the  effects  of  nonstationary  noise 
on  ARMA  models  and  concluded  that  the  prediction  formulas 
applied  in  the  stationary  case  were  also  valid  in  the 
nonstationary  case  provided  the  underlying  density 
distribution  of  the  nonstationary  process  was  symmetric. 
Discrete  density  distributions  for  400  seismic  signals  and 
their  associated  background  noise  were  plotted  and  in  all 
cases  the  density  distributions  were  symmetric.  This 
allows  us  to  conclude  that  an  AR  model  is  appropriate  for 
estimating  the  spectra  of  the  event  trace  given  our 
initial  assumptions. 

Description  of  Method 

The  first  step  in  developing  a  synthetic  waveform  was 
to  compute  the  coefficients  for  an  eighth  order  AR  model 
using  the  Burg  algorithm.  See  appendix  B  for  the 
subroutines,  "Polr_ArBURG"  and  "AR_SPECTRA" ,  to  compute 
the  Burg  coefficients  and  the  AR-spectra.  Coefficients 
were  computed  for  the  first  800  samples  of  the  event  file 
representing  the  background  noise  and  the  next  800  samples 
of  the  event  file  thought  to  contain  the  first  arrival 
phase.  Then  the  AR(8)  spectra,  p(f),  was  computed  for 
both  the  signal  and  the  noise  and  plotted  (Priestley, 
1981 )  . 
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22  «s  variance  of  the  white  noise 
f  *  frequency  Hz 

To  determine  the  shape  of  the  synthetic  waveform's  AR 
power  spectral  density  we  employ  the  doctrine  of 
parsimony.  The  minimum  number  of  poles  needed  to  describe 
the  spectral  shape  with  the  information  we  have  available 
is  three.  The  poles  are  located  at  the  low  cut-off,  high 
cut-off  and  corner  frequencies  of  the  signal  bandwidth. 
The  corner  frequency  is  assumed  to  be  the  frequency  with 
the  signal  spectra's  peak  amplitude.  The  low  and  high 
frequency  cut-offs  define  the  signal  spectra  interval 
where  the  signal  spectra  is  consistently  above  the  noise 
spectra.  The  amplitude  of  the  low  frequency  cut-off  is 
determined  by  the  roll-off  below  the  corner  frequency  due 
to  the  sensivity  of  the  recording  instrument, 
approximately  6dB  per  octave  for  the  Lajitas  station 
Instrument.  The  amplitude  of  the  high  frequency  cut-off 
is  determined  by  the  roll-off  above  the  corner  frequency 
due  to  attenuation.  The  expected  roll-off  is  -3dB  to  - 
4d8  per  octave.  Figure  21  illustrates  the  proposed 
method . 


The  coefficients  of  the  synthetic  AR(p)  model  are 
derived  from  the  synthetic  AR  power  spectral  density, 
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£(f),  described  in  the  preceding  paragraph.  The  AR 
theoretical  power  spectral  density,  p(f),  is  defined  by 


where 


P(  t ) 


S(f ) 


2 


and 


S(f ) 


S(f )S*(f ) 


S(f)  = _ V2T _ 

[1  -  0ie-i2Tf - Qpe-i2Tfpj 

Now  by  making  use  of  the  fact  that  the  theoretical  power 
spectral  density,  p(f),  can  also  be  written  in  terras  of 
the  autocovariance,  R(k),  ^ 

p(f)  *  f  R(k)dk: 

— 

we  can  obtain  the  autocovariance,  R(k),  of  the  synthetic 
AR  model  by  taking  the  inverse  fast  fourier  transform 
(IFFT)  of  the^synthetic  AR  power  spectral  density,  p(f). 

p(f)  =  J  e“i27rfkR(k)dk  >-LFFT..  J'  x(k)x(k-t)dt  =  R(k) 

Then  the  autocovariance,  R(k) ,  of  the  synthetic  process, 
x(t),  is  used  by  the  Yule-Walker  algorithm  (Kay  and 
Marple,  1981),  (see  subroutine  "AKA IKE"  in  appendix  B),  to 
obtain  the  AR(p)  coefficients  for  an  AR  model  of  order  p 
(determined  by  Akaike's  information  criteria  ,AIC, 
Priestley,  1981).  This  synthetic  AR(p)  model  is  then 
convolved  with  a  spike  plus  white  noise  to  create  the 
synthetic  waveform.  The  synthetic  waveform  is  then  cross- 
correlated  with  the  event  trace  to  try  and  pick  out  the 
first  arrival  phase. 


■aw 


'  ''  ^ ^  V  WI.TB WWW  WV 


Effectiveness  of  Method  and  Conclusions 

The  method  hypothesized  above  started  to  break  down 
with  the  assumptions  that  the  first  arrival  phase  has  the 
same  period  as  the  corner  frequency  and  the  background 
noise  had  no  obvious  corner  frequency.  The  period 
estimated  by  the  analyst  for  the  first  arrival  phase  of 
the  Lajitas  event  traces  only  correlated  well  with  the 
corner  frequencies  of  the  corresponding  AR(8)  spectral 
estimates  in  25*  of  the  events  studied;  i.e.,  teleseismic 
events  with  a  high  signal  to  noise  (S/N)  ratio. 

Often  a  corner  appeared  in  the  background  noise 
making  it  difficult  to  determine  the  correct  spectral 
bandwidth  for  the  signal  automatically.  For  the  above 
reasons  this  method  for  identifying  the  first  arrival 
phase  for  a  broad  category  of  seismic  events  was  abandoned 
in  favor  of  a  more  robust  non-parametric  technique,  the 
rank  sum  detector.  Representative  examples  of  the  AR(8) 
spectral  estimations  for  several  of  the  events  studied  are 
shown  in  figures  22  through  26  with  the  analyst's 
estimations  of  the  dominate  periods  for  the  noise  and  the 
first  arrival  phase  marked  for  the  Lajitas  seismic  events. 
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Figure  22. a 


Figure  22.  AR( 8 } -Spectral  Estimates  for  the  Signal 
and  Noise  for  Event  1.  (a)  Event  traces  recorded  at  each 
of  the  four  stations,  Lajitas,  Marathon,  Shatter  and  Tres 
Cuevas,  for  Event  1.  (b)  Lajitas  AR( 8 ) -Spectral 
Estimates.  The  Analyst  Estimates  for  the  Dominate  Period 
of  the  Noise  and  the  First  Arrival  Phase  are  Indicated 
with  Dashed  Lines  Drawn  to  their  Respective  Spectral 
Estimates.  (c)  Marathon  AR ( 8 ) -Spectral  Estimates.  <d) 
Shafter  AR(8)  Spectral  Estimates.  (e)  Tres  Cuevas  AR(8) 
Spectral  Estimates. 
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Figure  23.  a 


Figure  23.  AR( 8 ) -Spectral  Estimates  for  the  Signal 
and  Noise  for  Event  5  are  Normalized  by  the  Variance  of 
the  Noise.  (a)  Event  Traces  Recorded  at  Each  of  the  Four 
Stations,  Lajitas,  Marathon,  Shafter  and  Tres  Cuevas,  for 
Event  5.  (b)  Lajitas  AR ( 8 ) -Spectral  Estimates.  The 
Analyst  Estimates  for  the  Dominate  Period  of  the  Noise  and 
the  First  Arrival  Phase  are  Indicated  with  Dashed  Lines 
Drawn  to  their  Respective  Spectral  Estimates.  (c) 
Marathon  AR ( 8 ) -Spectral  Estimates.  (d)  Shafter  AR(8) 
Spectral  Estimates.  (e)  Tres  Cuevas  AR(8)  Spectral 
Estimates . 
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Figure  24. a 


Figure  24.  AH ( 8 ) -Spectral  Estimates  for  the  Signa_ 
and  Noise  for  Event  6  are  Normalized  by  the  Variance  of 
the  Noise.  (a)  Event  Traces  Recorded  at  Each  of  the  Four 
Stations,  Lajitas,  Marathon,  Shatter  and  Tres  Cuevas,  for 
Event  6.  (b)  Lajitas  AR ( 8 ) -Spectral  Estimates.  The 
Analyst  Estimates  for  the  Dominate  Period  of  the  Noise  and 
the  First  Arrival  Phase  are  Indicated  with  Dashed  Lines 
Drawn  to  their  Respective  Spectral  Estimates.  (c) 
Marathon  AR ( 8 ) -Spectral  Estimates.  (d)  Shafter  AR(8) 
Spectral  Estimates.  (e)  Tres  Cuevas  AR(8)  Spectral 
Estimates . 


♦**♦  WISE  AR(  81-VELDCITT  SPECTRA 
_  SIGNH.  NOISE  ARl  81-VELDCITT  SPECTRA 


IMSWUS 


I  1  I ' 1 - 1  I  ■ 

“  ~  **  “  n&Bcr  “ 

Figure  2^.b 


1 


START  TINE:  1984  289  23:12:47.61  H42832313.R2I  REGIONAL 


**♦♦  WISE  ARt  8)-VELDCITT  SPECTRA 

_  SIGWR.  NOISE  ARC  81-VELDCITT  SPECTRA 

eeesAff'iEs 


FREOiScT 


Figure  2k. c 

59 


'  /  / 

N  S.  %  % 


>^4  ^  |6  N  ""fc  *VA  ^.VVV,.  ’»V '  *“ 


START  TI«:  1304  283  Z3: 12*44.4*  24283231 J. *18  RESICMf. 


***  noise  ar<  ai-vELxm  spectra 

_  SIGNAL  *4CIS£  AR!  81-VELOCITT  SPECTRA 

WBSAPPLES 


U.4 

FWfluexrr 


Figure  24.  d 


START  TI«t  1984  289  23112:33.82  T42S92313.W8  REGICm. 


— ♦  NOISE  AR (  8)-VEL0CITT  SPECTRA 

_  SIGNAL  .NOISE  AR!  81-VELDCITT  SPECTRA 

8BRSAPPLES 


t  ».m 

FKEJJJtHCT 


Figure  24. e 


START  TI«:  1964  290  0<19:17.  1  742900019.051  TBJESE1311C 


START  TIME:  1964  296  0:29:49.81  S4 2986021 . 823  TELESEIS4IC 


START  TI«:  1984  290  0:19:20.21  H42900819.053  TELESEISMIC 


START  TI«:  1964  290  0:19:17.  1  L42900019.0S1  TELJE5EI31IC 


A*x*» 


I - 1  I  I  I  I 

9.00  6.67  13.33  20.00  26.67  33.33 

SECONDS 


f  !  |  I 

40.00  46.67  S3. 33  60.00 


Figure  25. a 


Figure  25.  AR ( 8 ) -Spectral  Estimates  for  the  Signal 
and  Noise  for  Event  8  are  Normalized  by  the  Variance  of 
the  Noise.  (a)  Event  Traces  Recorded  at  Each  of  the  Four 
Stations,  Lajitas,  Marathon,  Shafter  and  Tres  Cuevas,  for 
Event  8.  (b)  Lajitas  AR{ 8 ) -Spectral  Estimates.  The 
Analyst  Estimates  for  the  Dominate  Period  of  the  Noise  and 
the  First  Arrival  Phase  are  Indicated  with  Dashed  Lines 
Drawn  to  their  Respective  Spectral  Estimates.  (c) 
Marathon  AR( 8 ) -Spectral  Estimates.  (d)  Shafter  AR(8) 
Spectral  Estimates,  { e )  Tres  Cuevas  AR(8)  Spectral 
Estimates . 
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Figure  26.  AR( 8 ) -Spectral  Estimates  for  the  Signal 
and  Noise  for  Event  11  are  Normalized  by  the  Variance  of 
the  Noise.  (a)  Event  Traces  Recorded  at  Each  of  the  Four 
Stations,  Lajitas,  Marathon,  Shatter  and  Tres  Cuevas,  for 
Event  11.  (b)  Lajitas  AR ( 8 ) -Spectral  Estimates.  The 
Analyst  Estimates  for  the  Dominate  Period  of  the  Noise  and 
the  First  Arrival  Phase  are  Indicated  with  Dashed  Lines 
Drawn  to  their  Respective  Spectral  Estimates.  (c) 
Marathon  AR ( 8 ) -Spectral  Estimates.  (d)  Shafter  AR(8) 
Spectral  Estimates.  (e)  Tres  Cuevas  AR(8)  Spectral 
Estimates . 
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Non-pa rame trie  Detectors 

Description  of  the  Detector  Used  to  Implement 
the  Non-parametric  Tests 

The  non-parametric  tests  used  to  distinguish  signal 
from  noise  in  the  detector  were  performed  on  the 
"modified"  slope  of  the  event  trace  observation  data. 
This  "modified"  slope  is  obtained  by  discrete  integration 
of  the  slope  of  the  event  trace  over  segments  with  the 
same  slope  direction.  Figure  20  illustrates  how  the 
slope  and  "modified"  slope  were  calculated  for  each  data 
point  in  the  event  traces.  The  subroutine,  "SLOPE",  in 
appendix  B  calculates  the  "modified"  slope  of  the  event 
trace . 

The  number  of  observations  selected  for  testing  the 
performance  of  each  non-parametric  technique  as  a  detector 
was  100.  The  sample  size  was  chosen  so  that  the  time 
period  of  each  observation  window  would  be  greater  than 
the  largest  period  expected  to  occur  in  the  event  trace. 
The  periods  of  the  various  components  of  the  event  trace 
signal  and  noise  range  from  0.1  to  1.2  seconds.  A  window 
of  100  observations,  taken  at  a  sample  rate  of  40  samples 
per  second,  correspondes  to  two  and  one-half  seconds.  The 
first  100  observations  of  each  event  trace  is  assummed  to 
be  a  representative  example  of  the  background  noise  for 
the  time  period  of  the  event  trace.  The  first  100 
observations  representing  background  noise  are  denoted  by 
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X(l),  1=1,100,  in  the  descriptions  of  the  non-parametr ic 
tests  following  this  description  of  the  detector.  The 
groups  of  100  sequential  observations  taken  at  later  times 
in  the  event  trace  and  compared  with  the  background  noise, 
X ( i ) ,  are  denoted  by  Y(i),  i=l,100.  The  null  hypothesis 
used  for  each  non-parametric  test,  HQ  t  ±s  the  assumption 
that  both  X(i)  and  Y(i)  are  taken  from  the  same 
population;  i.e.,  that  the  cummulatlve  distribution 
functions  for  X{i)  and  Y(i),  F(x)  and  P{y) ,  respectively 
are  the  same . 

Since  we  know  the  first  30  seconds  into  the  event 
trace  precedes  the  6.4  second  window  which  triggered  the 
fast  Walsh  detector,  we  make  the  assumption  that  the  first 
arrival  phase  of  the  seismic  event  (P  or  Lg)  is  within  the 
first  60  seconds  of  the  event  trace.  To  detect  the 
seismic  signal  a  window  of  100  observations  is  moved  in 
increments  of  10  samples,  0.25  seconds,  down  the  event 
trace.  Each  window  of  100  samples  is  compared  to  the 
first  100  samples  of  the  trace  representing  background 
noise  using  a  non-parametric  statistical  test.  Detection 
occurs  when  the  non-parametric  test  indicates  the  two 
windows  of  100  samples  are  not  from  the  same  population; 
i.e.,  the  underlying  distributions,  F(x)  and  F(y),  are  not 
the  same. 

Three  non-parametric  tests,  (1)  the  two  sample  sign 
test,  (2)  the  run  test,  and  (3)  the  rank  sum  test,  were 
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tried  in  the  detector  to  determine  the  most  effective  test 
for  discriminating  between  signal  and  noise.  The 
following  describes  each  test,  the  assumptions  each  test 
made,  the  validity  of  those  assumptions  with  respect  to 
the  data  set  and  the  effectiveness  of  each  test  in 
discriminating  between  the  background  noise  and  the  signal 
plus  noise.  The  detector  subroutine,  "DETECTOR",  is 
listed  in  appendix  B. 


Two  sample  sign  test 


Assumptions 


The  data  for  the  two  sample  sign  test  consists  of  two 
random  samples,  N  observations  from  the  control 
population;  i.e.,  background  noise,  and  N  observations 
from  the  treatment  population;  i.e.,  signal  plus  noise. 
In  our  case  we  choose  N  observations,  X( 1 ) , . . . ,X(N) ,  from 
the  beginning  of  the  event  trace  representative  of  the 
background  noise  and  N  observations,  Y ( 1 ) , . . . , Y(N) ,  from 
the  remaining  portion  of  the  event  trace.  Unfortunately, 
the  observations  taken  from  a  seismic  signal  are  not 
independent  random  variables.  Instead  there  is  a 
dependency  between  the  observations.  Walsh  (1949,1951) 
has  shown  that  the  sign  test  will  have  similar  results  if 
the  observations  in  the  two  samples  are  mildly  dependent. 

The  two  sample  sign  test  is  sensitive  to  changes  in 


the  location  and  spread  of  the  distribution.  If  a  signal 
is  present  the  spread  of  the  X(i)  observations  should  be 
significantly  different  from  the  spread  of  the  background 


noise,  Y(i)  observations.  The  sign  test  run  on 
seismic  data  assumed  that  no  ties,  (i.e.,  X(i)  =  Y| 
were  present. 

Procedure 


To  test 


1 . 


Z(i)  = 


<5U)  = 


H0: 

F(x) 

«  F(y) 

variables 

1 

If 

[  Y(  i  )  - 

X(i)] 

>  0 

0 

if 

[  Y  (  i  )  - 

X(i)  ] 

=  0 

1 

if 

[  Y(  i  )  - 

X(i)  3 

<  0 

1 

if 

Z(i)  > 

0 

0 

if 

Z(i)  < 

0 

Set 


N 


Sn  -  Z  ( 1  )  6(1)  . 

i — i 

The  statistic  Sn  is  the  number  of  positive  Z's. 

3.  For  a  one-sided  test  of  HO  versus  the  alternative, 

Ha :  F(x)  =  F(y) , 

at  the  alpha  level  of  significance, 

Reject  Hq  lf  Sn  <  [N  -  b( alpha , N , 1 /2 ) ] 
Accept  HQ  if  Sn  >  [N  -  b( alpha ,N, 1/2 )] , 
where  the  constant  b( alpha , N , 1/2 )  under  the 


the 
i)  )  . 


null 


Si 


hypothesis  satifies  P[Sn  >  b(alpha,N, 1/2) ]  =  alpha.  That 
is,  b(alpha,N, 1/2)  is  the  upper  alpha  percentile  point  of 
the  binomial  distribution  with  sample  size  N  and  p  =  1/2. 

For  large  sample  approximation  under  the  null 
hypothesis  define: 

Sn*  =  Sn  -  E  ( Sn)  «=  Sn  -  N/2 
SQRT[ Var ( Sn) ]  SQRT(N/4) 

When  Hq  is  true,  the  statistic  Sn  has  an  asymptotic  (N 
tending  to  infinity)  standard  normal  distribution, 
z(alpha).  The  normal  theory  approximation  test  of  HQ 
versus  the  alternative  Hq/  at  the  alpha  level  of 
significance  is 

Reject  Hq  if  sn*  >  z(alpha) 

Accept  Hq  if  Sn*  <  z(alpha). 

In  regard  to  ties,  X(i)  =  Y(j),  if  there  are  zero 
values  among  the  Z's,  discard  them  and  redefine  N  to  be 
the  number  of  nonzero  Z's.  Figure  27  is  an  example  of  the 
two  sample  sign  test  with  N  =  30,  k  =  10  for  a  confidence 
interval  (0.137-0.583),  and  alpha  «  0.01. 

Effectiveness 

A  large  number  of  ties  occur  when  the  X(i)  and  Y(i) 
observations  taken  from  seismic  signals  are  compared. 
These  ties  are  due  to  the  limited  dynamic  range  for  small 
values  of  seismic  observations  near  zero;  i.e.,  the  same 
values  occur  more  often  near  zero  since  there  is  a  more 
limited  range  of  values  to  choose  from.  Since  the  ties 
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Figure  28.  The  Performance  of  the  Two  Sample  Sign 
Test  When  Used  to  Detect  Seismic  Events.  The  Null 
Hypothesis  for  this  Test  Assumes  ail  the  Observations, 
X ( i )  and  Y(i),  Come  from  the  Same  Population,  the 
Background  Noise.  The  Blue  Color  Indicates  a  Rejection  o 
the  Null  Hypothesis,  i.e.  a  Signal  is  Present.  (a)  Even 
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Figure  29.  The  Performance  of  the  Two  Sample  Sign 
Test  When  Used  to  Detect  Seismic  Events.  The  Null 
Hypothesis  for  this  Test  Assumes  all  the  Observations, 
X(i)  and  Y(i),  Come  from  the  Same  Population,  the 
Background  Noise.  The  Blue  Color  Indicates  a  Rejection  of 
the  Null  Hypothesis,  i.e.  a  Signal  is  Present. 


were  not  taken  into  consideration  in  calculating  the  value 
of  the  sign  test  they  dilute  its  sensitivity  to 
differences  in  the  spread  of  the  distributions.  Ties  were 
not  accounted  for  because  of  the  increase  in  computation 
time  it  would  require;  i.e.,  recalculation  of  the  test  for 
the  reduced  number  of  untied  observations.  Figures  28  and 
29  illustrate  the  performance  of  the  two  sample  sign  test 
when  used  to  detect  signals.  The  sign  test  subroutine, 
"SIGNEXP" ,  is  listed  in  appendix  B. 

Run  test 

Assumptions 

The  run  test  selects  N  observations,  X(l) . X(N), 

of  representative  background  noise  taken  from  the 
beginning  of  the  event  trace  and  M  observations, 
Y ( 1 Y (M) ,  at  some  time  later  in  the  event  trace. 
Then  the  observations  are  ordered  in  ascending  order  and 
the  number  of  runs  (groups  of  X  or  Y  observations)  are 
counted.  The  run  test  in  this  application  assumes  the  two 
samples,  X  and  Y,  are  independent  random  variables.  Since 
there  is  a  dependence  between  the  observations  in  each 
sample  the  designated  level  of  significance  for  the  run 
test  will  not  be  preserved. 

If  X(i)  and  Y(j)  are  from  the  same  population  then  X 
and  Y  will  be  well  mixed  and  the  number  of  runs  will  be 
large.  However,  if  X(i)  and  Y(j)  are  from  widely  separate 
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populations  then  there  will  only  be  two  runs.  The  run 
test  is  sensitive  to  differences  in  both  shape  and 
location  of  the  distributions. 


Procedure 


To  test 


H0:  F(x)  =  F(y)  . 

1.  Order  the  observations,  X(i),  1=1, N  and  Y(i),  i=l,M, 
together  in  ascending  order  from  least  to  greatest. 

2.  Set  Z  equal  to  the  number  of  distinct  groups  of  Y's. 
The  statistic  Z  is  the  number  of  runs  in  the  ordered  array 
of  observations. 

3.  For  a  one-sided  test  of  HQ  versus  the  alternative 

Ha:  F(x)  =  F(y )  , 

at  the  alpha  level  of  significance, 

Reject  HQ  if  z  >  zO 
Accept  HQ  if  z  <  zO  , 

where  the  constant  zO  is  the  largest  integer  which 
satisfies  . 


P[Z  =  z]  =  alpha. 


That  is,  zO  is  the  lower  percentile  point  of  the 
distribution  of  Z  for  sample  sizes  N  and  M. 

The  following  specifies  P[Z  =  z]  under  the  null 
hypothesis  so  that  we  can  determine  the  integer  zO  for  a 
given  test  size. 
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Figure  30.  Example  of  the  Run  Test.  The  X ( i )  and 
Observations  are  Ordered  in  Ascending  Order.  Then 
the  Number  of  Runs  of  X  and  Y  Values,  Z,  are  Counted.  If 
the  Number  of  Runs  is  Less  than  the  Threshold,  zO,  the 
Null  Hypothesis  that  X  and  Y  Both  Come  from  the  Same 
Population  is  Rejected. 
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Fig-ore  31. a 

Figure  31.  The  Performance  of  the  Run  Test  When  Used 
to  Detect  Seismic  Events.  The  Null  Hypothesis  for  this 
Test  Assumes  all  the  Observations,  X(i)  and  Y(i),  Come 
from  the  Same  Population,  the  Background  Noise.  The  Blue 
Color  Indicates  a  Rejection  of  the  Mull  Hypothesis,  i.e.  a 
Signal  is  Present.  (a)  Event  19.  (b)  Event  14. 
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Figure  32.  The  Performance  of  the  Run  Test 
to  Detect  Seismic  Events.  The  Null  Hypothesis 
Test  Assumes  ail  the  Observations,  X(i)  and  Y( 
from  the  Same  Population,  the  Background  Noise. 
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the  effectiveness  of  the  run  test  in  distinguishing  the 
differences  in  shape  between  two  different  populations. 
Since  this  non-parametric  technique  does  not  have  a  good 
way  of  handling  the  number  of  ties  naturally  occurring  in 
seismic  signal  comparisons,  this  technique  was  rejected  in 
favor  of  the  rank  sum  test.  Figures  31  through  32 
illustrate  the  performance  of  .the  run  test  as  a  signal 
detector . 

Rank  sum  test 

Assumptions 

The  rank  sum  test  selects  N  sequential  observations, 

X( 1) _ ,X(N) ,  of  the  "modified"  slope  of  the 

representative  background  noise  from  the  beginning  of  the 
event  trace.  Then  M  sequential  observations  of  "modified" 
slope,  Y ( 1 Y ( M ) ,  are  chosen  from  some  time  later  in 
the  event  trace.  The  N  plus  M  observations  are  ordered  in 
ascending  order  and  assigned  a  rank  based  on  their 
position  in  the  ordered  sequence.  The  ranks  for  the  Y(i) 
observations,  Ry(i),  are  summed  and  the  absolute  value  of 
the  difference  between  that  rank  sum,  Ty,  and  the 
estimated  mean,  E(Ty),  is  compared  with  a  predetermined 
threshold  value,  k. 

The  rank  sum  test  assumes  the  X(i)  and  Y(i) 
observations  are  independent  random  variables.  Since 
background  noise  is  not  a  purely  random  process  and  the 


seismic  signal  can  be  thought  of  as  the  background  noise 
"treated"  with  the  addition  of  a  seismic  event,  neither 
the  X(i)  or  the  Y ( i )  observations  are  independent  random 
variables.  Serf ling  (1968)  investigating  the  robustness 
of  Wilcoxon  test,  upon  which  the  rank  sum  test  is  based, 
let  the  two  samples,  X  and  Y,  be  independent  of  each  other 
but  let  the  random  variables  within  a  sample  be  possible 
dependent.  It  is  found  that  the  robustness  of  the  test 
statistic  for  the  Wilcoxon  two-sample  procedure  under  the 
null  hypothesis  with  departures  from  the  standard 
assumption  of  random  samples  depends  upon  the  grade  of 
correlation  of  the  variables  X(i)  and  X(i+1).  In  other 
words,  similar  results  for  the  rank  sum  test  should  be 
obtained  when  only  mild  dependence  occurs  between  the 
observations  in  the  samples.  Since  we  do  not  know  the 
exact  degree  of  dependence  between  the  observations,  there 
is  a  probability  that  the  significance  level,  (alpha), 
assumed  for  the  rank  sum  test  does  not  reflect  the  true 
significance  level. 

Procedure 

To  test 

H0:  F(x)  =  F ( y )  . 

1.  Order  the  N  and  M  observations  from  least  to  greatest 
and  let  Ry(i)  denote  the  rank  of  Y(i)  in  this  ordering. 
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The  statistic  Ty  is  the  sum  of  the  ranks  assigned  to  the 
Y's  . 

3.  For  a  one-sided  test  of  HO  versus  the  alternative 

Ha:  F(x)  =  F ( y ) , 

at  the  alpha  level  of  significance. 

Reject  Hq  if  Ty  >  W(alpha,M,N) 

Accept  Hq  if  Ty  <  w(alpha,M,N)  , 

where  the  constant  w(alpha,M,N)  under  the  null  hypothesis 


satif ies 


P[Ty  >  w ( alpha , M , N ) ]  =  alpha. 


Values  of  w(alpha,M,N)  are  given  in  Table  C.l,  Appendix  C 
(Hollander  and  Wolfe,  1973). 

For  large  sample  approximation  under  the  null 
hypothesis  define 


Ty  =  Ty  -  E(Ty) 

SQRT[Var (Ty) ] 


When  Hq  is  true,  the  statistic  Ty  has  an  asymptotic 
(minimum  of  N  or  M  tending  to  infinity)  standard  normal 
distribution.  The  one-sided  normal  approximation  theory 


for  the  test  statistic  is 


Reject  H( 
Accept  H 


if  Ty  >  z (alpha ) 
if  Ty*  <  z( alpha ) . 


For  M  >  7  and  N  >  7  the  asymptotic  normal  approximation  is 
quite  accurate.  (Mood,  GraybilJ  and  Boes ,  1974) 

If  there  are  ties;  i.e.,  the  ith  observation  in 


ascending  order  is  tied  with  the  k  observations  following 
the  ith  observation  then  the  summation  of  the  indices 
divided  by  k+1  is  the  average  rank  assigned  to  each  of  the 
tied  values.  For  the  large  sample  approximation,  compute 
Ty  using  average  ranks,  and  replace  Var(Ty)  by 


where  L  is  the  number  of  tied  groups  of  ranks  and  f(i)  is 
the  size  of  the  ith  tied  group.  An  untied  observation  is 
considered  to  be  a  tied  group  of  size  1.  Hence  if  there 
are  no  tied  observations  the  right  hand  side  of  Var(Ty) 
reduces  to 

MN  (M +N+1) 

12 

Figure  33  illustrates  an  example  of  the  rank  sum  test 
for  N=30,  M=30 ,  and  alpha  *  0.001.  The  rank  sum  test 

subroutine,  "RANKTEST" ,  is  listed  in  appendix  B. 

Effectiveness 

Since  the  ties  are  handled  by  averaging  the  ranks 
this  test  is  sensitive  to  the  relative  changes  in  the 
shape  of  the  seismic  signal's  density  distribution. 
Figures  34  and  35  illustrate  the  application  of  the  rank 
sum  test  to  detect  a  seismic  event. 

Conclusions 


The  rank  sum  test  was  the  most  effective  of  the  non- 
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Figure  33.  Enanple  of  the  Rank  Sun  Test.  The  X(i) 
and  Y{  i )  Observations  are  Ordered  in  Ascending  Order  am 
Assigned  at  Rank  Indicating  their  Position  in  the 
Sequence.  Then  the  Ranks,  Ry ( i ) ,  for  the  Y(i)  Values  are 
Sunned,  Ty .  If  the  Absolute  Value  of  the  Rank  Sun,  Ty , 
minus  its  Enpected  Value,  E{Ty) ,  is  Less  than  the 
Threshold,  K.  we  Accept  the  Null  Hypothesis  that  X(i)  and 
Y(i)  Come  fret  the  Same  Papulation. 
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Figure  34.  The  Performance  of  the  Rank  Sum  Test  When 
Used  to  Detect  Seismic  Events.  The  Null  Hypothesis  for 
this  Test  Assumes  ail  the  Observations,  X(i)  and  Y(i), 
Come  from  the  Same  Population,  the  Background  Noise.  The 
Blue  Color  Indicates  a  Rejection  of  the  Null  Hypothesis, 
i.e.  a  Signal  is  Present.  (a)  Event  19.  (b)  Event  14. 
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Figure  35.  The  Performance  of  the  Rank  Sum  Test  Whe.i 
Used  to  Detect  Seismic  Events.  The  Null  Hypothesis  for 
this  Test  Assumes  ail  the  Observa1- ions  ,  X(i'  and  Y(i)  , 
Come  from  the  Same  Population,  tne  Background  Noise.  The 
Blue  Color  Indicates  a  Rejection  of  the  Null  Hypotnesis, 
i.e.  a  Signal  is  Present. 
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parametric  tests  in  discriminating  between  signal  and 
noise.  After  studying  figures  28,29,31,32,34  and  35  one 
can  see  that  the  rank  sum  test  does  the  best  Job  of 
discriminating  the  event  signal  from  background  noise  in  a 
trace  where  an  event  signal  is  known  to  be  present.  The 
success  of  the  rank  sum  test  is  due  to  its  method  of 
handling  the  large  number  of  ties  indigenous  to  seismic 
signals . 
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DETECTION  TECHNIQUE  AND  FIRST  ARRIVAL  PICKER 
EMPLOYED  TO  FIND  THE  FIRST  ARRIVAL  PHASE  OF 
EARTHQUAKES  AND  EXPLOSIONS  IN  SEISMIC 
SIGNALS  CONTAMINATED  BY  NOISE 

Description  of  Non- parametric 
Detector  and  Picker 

The  rank  sum  test  applied  to  the  "modified"  slope  of 
the  event  trace  observation  data  most  clearly 
distinguished  seismic  signals  from  the  background  noise  of 
the  three  non-parametric  techniques  tested.  Several 
Implementations  of  a  detector  based  on  the  rank  sum  test 
were  run  on  a  training  set  of  48  event  traces.  The  level 
of  significance  for  the  rank  sum  test  was  varied  to  find 
the  threshold  for  detection  which  would  allow  the  rank  sum 
test  to  clearly  discriminate  between  signal  and  noise. 
However,  the  "unmodified"  application  of  the  rank  sum  test 
to  discriminate  between  signal  and  noise  was  not 
"consistent"  in  picking  the  seismic  event  due  to  varying 
S/N  ratios.  To  allow  the  threshold  for  detection  to 
change  with  S/N  ratio  a  "modified"  version  of  the  rank  sum 
test  was  Implemented. 

The  minimum  and  maximum  rank  sums  for  a  each  trace 
were  determined  in  this  "modified"  implementation  of  the 
rank  sum  test.  If  the  range  of  the  rank  sums,  (i.e.,  sum 
of  the  ranks  in  each  window),  for  a  trace, 


maximum  rank  sum  -  minimum  rank  sum, 

Is  greater  than  3600,  which  was  observed  to  constitute  the 
"high"  S/N  case,  then  the  threshold  for  signal  detection 
is  computed  to  be  the  minimum  rank  sum  plus  2700.  If  the 
range  of  the  rank  sums  is  less  than  3600,  which  was 
observed  to  constitute  the  “low"  S/N  case,  then  the  signal 
detection  threshold  is  computed  to  be  the  minimum  rank  sum 
plus  75  percent  of  the  range  of  the  rank  sums .  Based  on  a 
detailed  analysis  of  the  rank  sums  test  when  applied  to 
the  event  traces  in  the  data  set,  this  allows 
approximately  99.9  percent  of  the  rank  sums  due  to 
background  noise  to  fall  below  the  threshold  in  the  "high" 
S/N  cases  and  99  percent  in  the  "low"  S/N  cases. 

The  automatic  first  arrival  pick  is  defined  as  the 
first  zero-crossing  preceding  the  first  observation  in  the 
first  detection  window  whose  modified  slope  Is  at  least 
1.05  times  greater  than  the  maximum  modified  slope  value 
of  the  representative  background  noise. 

The  event  detector  based  on  the  rank  sum  test  takes 
the  first  N=100  observations  in  the  event  trace  to  be  a 
representative  sample  of  the  background  noise  present  in 
that  event  trace.  The  rank  sum  test  is  computed  for 
observation  windows,  N  equal  100  samples  in  length,  moved 
in  increments  of  10  samples,  0.25  seconds,  down  the  event 
trace.  The  rank  sum  of  each  window  is  computed  and 
compared  with  the  rank  sum  of  the  background  noise  window. 


The  subroutine  to  detect  and  pick  the  first  arrival  phase, 
"RANK2700",  is  listed  in  appendix  B. 

Results 

The  detector  which  used  the  "modified"  version  of  the 
rank  sum  test  described  in  the  preceding  section  to 
discriminate  between  signal  and  noise  was  run  on  152  of 
the  event  files  recorded  during  the  GSETT  experiment. 
These  event  files  were  selected  to  give  the  broadest 
representation  of  the  wide  variety  of  seismic  signals  we 
would  expect  to  encounter  if  the  detector  were  run  over 
all  events  recorded  by  the  fast  Walsh  detector  on  a  daily 
basis . 

To  analyze  the  performance  of  the  "modified"  rank  sum 
detector  we  looked  at  128  of  the  event  traces.  Since  24 
of  the  152  event  traces  the  detector  was  run  on  were  also 
included  in  the  training  set  used  to  develop  the 
"modified"  version  of  the  rank  Siam  detector  they  were 
excluded  in  this  analysis  to  give  us  an  unbiased  look  at 
the  detector's  performance.  One  hundred  and  twenty-four 
automatic  picks  were  compared  to  their  respective 
analyst's  picks.  Sixty-five  were  within  0.4  seconds  of  the 
analyst's  pick.  After  analizing  the  automatic  picks 
differing  from  the  analyst's  picks  by  more  than  0.4 
seconds,  it  was  concluded  that  5  appeared  to  be  more 
correct  than  the  analyst's  original  picks.  Twenty-four  of 
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the  automatic  picks  were  ambiguous,  either  the  automatic 
pick  or  the  analyst  pick  could  be  correct.  Thirty-five  of 
the  automatic  picks  had  errors  larger  than  0.4  seconds 
relative  to  the  analyst  picks:  4  local;  28  regional;  and  3 
teleseisaic  events.  The  representative  background  noise 
for  22  of  the  regional  and  local  automatic  picks  with 
errors  larger  than  0.4  seconds  contained  P  arrivals.  This 
indicates  that  the  violation  of  the  assumption  that  the 
first  100  samples  is  background  noise  causes  a 
deterioration  in  the  accuracy  of  the  automatic  first 
arrival  pick.  The  apparently  erroneous  automatic  picks 
that  occur  for  teleseismic  events  are  caused  by  small 
emergent  precursors  to  the  P-arrival  which  the  analyst 
picked . 

Occasionally,  the  assumption  that  the  first  100 
observations  in  the  event  trace  represent  background  noise 
is  violated.  Either  because  the  fast  Walsh  detector  did 
not  detect  the  P  arrival  but  instead  triggered  on  the  Lg 
arrival  or  because  the  signatures  from  two  seismic  events 
occurring  near  the  same  time  overlap  each  other  on  the  the 
event  trace.  This  causes  either  the  P  arrivals  or  the 
coda  from  a  preceding  event  to  contaminate  the  first  100 
observations  assummed  to  contain  only  background  noise. 
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TABLE  4.  Errors  (Automatic  Minu9  Analyst  Picks  in 
Seconds ) ,  and  Ranges  of  the  Rank  Sums  for  the 
152  Events  Used  to  Test  RANK2700 


ALL  EVENTS 


EVENT  NAME 

RANGE 

EJ5RCR 

S4aaaut«.S37 

5057.0 

8.0 

142832313.303 

4802.5 

0.0 

H429USS1.0SI 

3742.0 

8.0 

742332333.319 

5323.0 

0.0 

T42S7381S.046 

5751.5 

8.0 

L429717S3.01S 

5272.0 

0.0 

342971753. 018 

6174.0 

8.0 

T 429020 19.351 

5564.5 

-0.1 

342308136. 331 

3581. a 

0.1 

M42S020S5.039 

4234.0 

0.1 

T42891418.034 

5343.5 

-0.  1 

H42312I 40. 040 

5143.0 

0.  1 

T42912213.8S7 

4014.5 

8.1 

1.42920526.308 

5145.5 

-0. 1 

342920528.301 

4184.0 

0. 1 

T42920S28. 008 

5536.0 

-0. 1 

L42921S33.036 

5469.0 

0.1 

M42321S33.023 

5871.0 

0.1 

342921533.326 

5261.5 

0.1 

L429218S9.031 

3338.0 

0.1 

T  42921859.  028 

4275.0 

0.1 

L42332129.009 

5495.5 

-0.  1 

T42932129. 006 

5737.0 

0.1 

L42SS1 418.037 

4920.5 

-0.1 

T 4235 1908. 057 

5539.5 

0.1 

L42S62037.016 

3720.5 

0.  1 

542970815.033 

5909.5 

0.1 

S42832313.316 

5454.0 

0.1 

M42892328. 007 

3311.5 

0.  1 

842300013.353 

4573.5 

0.1 

T 42302054.  044 

5904.0 

0.2 

142532303.019 

5848.0 

0.2 

842332303.235 

5889.5 

0.2 

342321430.030 

3994.0 

-0.2 

842351908. 050 

5725.0 

0.2 

S42911603.000 

4508.5 

0.2 

S 42300 02 1.023 

4569  5 

0.2 

842891 418. 0S3 

5475.5 

0.2 

1.42500013.051 

4862.0 

-0.2 

142902054.044 

5534.  a 

-3.2 

8423717S3.030 

5039.5 

0.2 

T428S2304.016 

5740.0 

0.2 

842901820.047 

4770.0 

0.3 

342311807.026 

4369.0 

8.3 

842921430.014 

4350.0 

0.3 

T429U  807.032 

4579.5 

-0.3 

L 42351308. 057 

5912.5 

8.3 

T42321429.018 

3944.0 

-0.3 

T42832327.ai2 

5168.0 

0.3 

T  42311851 . 207 

5053.5 

8.3 

T42911358.05S 

4692.  5 

8.3 

142321823.033 

4112.5 

8.3 

742891905.008 

3880.0 

8.3 

T42900136. 021 

4554.5 

0.3 

742302203.030 

3713.0 

0.3 

ALL  EVENTS 


EVENT  NAME 

RANGE 

EJtRCR 

T  42321 533.  036 

5751.5 

-0.4 

642362038.001 

4773.5 

0.4 

S42301821 . 025 

4187.0 

-0.4 

L4231 1827. 323 

45-2.5 

8.4 

L423121  40.015 

5730.0 

0.4 

K4231 1837. 323 

4144.5 

0.4 

L42312213. 357 

4443.5 

— .  4 

142321  423. 323 

2235.0 

0.5 

842352123.  331 

5070.0 

-0.5 

742832313.008 

sses.s 

0.6 

L4298224?. 213 

4145.3 

-0.6 

T  423 103 25.  029 

5643.0 

0.6 

542332129.318 

6325.5 

-0.5 

H42891815. 358 

3481.5 

0.6 

L*231 1851 .303 

4626.5 

-0.6 

842310325.376 

SS40.0 

0.7 

842302213. 873 

4052.0 

0.7 

L4291 1339. 217 

4390.0 

-0.7 

547367038.214 

4264.5 

B.a 

T4291 1632. 203 

519'. 5 

0.8 

542892327.218 

4036.0 

-0.8 

T 423022 46. 338 

3742.0 

0.9 

1142321830.208 

3958.5 

0.9 

842373815.043 

6135.5 

-0.9 

842311602.038 

3788.3 

0.9 

842971351.000 

2403.5 

-1.3 

U423738! 5.046 

6091.0 

-1.1 

342351309.000 

5251.5 

-1.1 

742921829.039 

4849.5 

-1.2 

L42831719.050 

4166.5 

1.2 

742891814.046 

4812.0 

1.2 

L 42371952. 327 

4583.0 

1.2 

S 42332303.  335 

5628. 5 

-1.2 

1.47310325. 

5353.0 

1.4 

742362377.016 

3655.0 

1.5 

842891719.057 

3826.0 

1.7 

142832327.016 

4615.5 

-1.8 

742891719.347 

4242.5 

1.9 

342321823.023 

5601.5 

-2.0 

84231221  4.092 

4067.0 

2.1 

342921858.056 

3923.5 

2.3 

S4231 1851. 011 

5646.0 

-2.4 

842831306.040 

3349.5 

2.7 

T  423121 40. 315 

5516.5 

-2.9 

L 42302203. 327 

2860.0 

3.0 

742971752.352 

5535.0 

3.3 

842320528.004 

5799.5 

-3.4 

742371351.000 

3419.5 

3.9 

842302246.029 

4954.5 

4  4 

842832704.028 

4193.5 

-5.1 

54233221 1 . 000 

2880.3 

5.1 

542302247.016 

2907.0 

5.4 

542911358. 358 

6332.0 

-6.3 

342892704.012 

5496.2 

-5.3 

L42991815.018 

417B.0 

6.5 

<c- 
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TABLE  4  (Continued) 


ALL  EVENTS 


EYEXT  NAME 

RANGE 

ERROR 

342631719. 025 

4605.5 

7.4 

L <2381621 . 016 

4405.5 

7.5 

T 4230 1620. 037 

3232.0 

-9.1 

342371350.040 

4355.5 

10.1 

342312133.030 

3740.0 

10.4 

342631615.016 

3470.5 

10.4 

H42300137. 013 

4208.0 

12.6 

H 423 11 353. 003 

4133.0 

13.4 

342631303.013 

3351.2 

_ ;  c  i 

H426  S2313. 020 

4760.5 

-15.3 

142631306.030 

3066.5 

-17.5 

L 42632204. 033 

5575. S 

-19.3 

342310325.016 

5376.0 

-22.3 

H42321&59.  054 

4645.5 

-23.6 

342312215.001 

5845.0 

-23.7 

L42900136. 053 

4227.5 

-25.0 

342302054.05! 

6aS4.2 

-26.3 

142311602.007 

4522.0 

-33.3 

NUMBER  Of  EVENTS  z  126 
HETIAN  =  a. 5 
*D*  =  3.2 

VARIANCE  =  44. 3 


LOCAL  EYEXTS 


EYEXT  MANE 

RANGE 

ERROR 

L42371753.018 

5272.0 

0.0 

542371753.018 

6174.0 

0.0 

14236 2037. 016 

3720.5 

0.1 

M  42302255. 038 

4234.0 

0.1 

H423121 40. 040 

5143.0 

0.1 

742892304.016 

5740.0 

0.2 

142302054.044 

5304.0 

0.2 

H4237 1753. 030 

5339.5 

0.2 

L42302254.044 

5304.0 

-8.2 

H  42362038. 00 1 

4733.5 

0.4 

L 42312140. 015 

5730.0 

0.4 

542562038.014 

4264.5 

0.8 

M42971S51.000 

2403.5 

-1.0 

L4297 1352. 007 

4583. 0 

1.2 

142362037.016 

3655.0 

1.5 

142312140.015 

5518.5 

-2.9 

142371752.052 

5505.0 

3.3 

142371951.000 

3419.5 

3.9 

H 4283230 4. 023 

4193.5 

-5.1 

S 42892304.  012 

5436.3 

-6.3 

S  4237  1  350. 040 

4355.5 

10.  1 

542912139.030 

3740.3 

10.4 

142892304.039 

5576. 5 

-19.3 

542382054.051 

6694.3 

-26.8 

NUMBER  Of  EVEXTS  = 

24 

HETlfiN  =  B.8 

HEBN  r  4.0 
VARIANCE  =  326.  3 
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TABLE  4  (Continued! 


REGiONAL  EVENTS 


REGIONH.  EVENTS 


EVENT  WWE 

Rf*GE 

E35KCR 

EVENT  NAME 

RANGE 

ERROR 

542331418.337 

5357.3 

0.3 

742831814.346 

4812.2 

1.2 

L42S32313.003 

4802.5 

8.3 

L4283171S.  353 

4166.  5 

1.2 

«42auasi.a5i 

3742.0 

3.8 

S42332303. 335 

5626.5 

-1.2 

T4233 2333. 219 

5323.3 

8.0 

742321823.333 

48*9.5 

-1.2 

542892313.016 

5454.3 

0.1 

L4231 3325. 332 

5353.3 

1.  * 

842892323.307 

3311.5 

0. 1 

842831713.057 

3825.3 

I  •  7 

S 42333 138.33! 

3581.0 

0. 1 

L42892327. 016 

4615.5 

-1.8 

T 4283! 413. 334 

5343.5 

-3. 1 

742831719.047 

4242.5 

1.3 

742312213.  aS~ 

4014.5 

8.1 

S42321823. 323 

5601.5 

-2.0 

L 42321533. 336 

5463.0 

0.1 

842312214.352 

♦067.3 

2.  '. 

M42321S33.323 

5871.0 

8.1 

542321853. 256 

3923.5 

2.3 

542321533.326 

5261.5 

0. 1 

542311851.011 

56*6. 3 

-2.  ♦ 

L42321853. 331 

3338.3 

0.1 

842831306.3*3 

33*9.5 

2.7 

T  4232 1 3SS.  328 

4275.3 

0. 1 

L4233222S. 327 

2863.3 

3.3 

L42332123.  333 

5435.5 

-0.  1 

842332246.323 

435*.  5 

«.  * 

T42332123.036 

5737.3 

3. 1 

542332211.333 

2883.3 

S.  1 

L42831 418. 337 

4323.5 

-0. 1 

S 423322 47. 316 

2337.3 

5.  4 

S42321 433. 333 

3334.3 

-3.2 

542311358.058 

6332.3 

-6.3 

842331418.353 

5*75. 3 

3.2 

L42£3!8!S.3'.8 

4)’8.0 

6.5 

L42332332.aid 

5848.3 

3.2 

542831713.025 

4825.5 

7.  * 

M42332333. 335 

5883.5 

3.2 

L42381821 .216 

4435.5 

7.  S 

342311633.333 

4538.5 

0.2 

742331823.037 

3232.3 

-3.  1 

T42321 423.31 8 

3344.3 

-0.3 

542831815.018 

3473.5 

10.* 

T42332233. 333 

3713.3 

0.3 

842308137.213 

4238.3 

12.6 

742332327.312 

5158.3 

0.3 

842311353.038 

4133.3 

13.  * 

S42S1 1837. 326 

4363.3 

3.3 

542831303.313 

3361.3 

-15. 1 

142321823.333 

4112.5 

3.3 

8428S2313. 323 

4760.5 

-15.3 

T42S1 1337. 332 

4573.5 

-0.3 

L42631336. 333 

3068.5 

-17.5 

742831335. 338 

saea.3 

0.3 

542910325.316 

5376.3 

-22.3 

74231 1851.337 

5053.5 

8.3 

842321853.35* 

4645.  5 

-23. 6 

T4231 1358. 356 

4632.5 

3.3 

542912215.331 

58*5.3 

-23.’ 

742303136.321 

4554.5 

3.3 

L429031 36. 353 

*227.5 

-25.2 

842921 438. 31 4 

4350.3 

0.3 

L4291 1632. 237 

*522.3 

-33. 3 

842331820.347 

4773.3 

0.3 

142311837.323 

4572.5 

3.  4 

HUMBER  Of  EVENTS  = 

88 

742321533.336 

5751.5 

-3.  4 

8ECJSN  =  0.6 

142912213.057 

4443.5 

-0.  4 

8EFW  =  3.5 

842911807.829 

4144.5 

3.4 

VFK1PNCE  -  135. 6 

S42901821 . 325 

4187.3 

-0.4 

L42321 423.323 

2235.3 

3.5 

842332123.331 

5073.3 

-3.5 

742310325. 323 

5643.3 

3.6 

342332123.318 

6385.5 

-0.6 

742832313.338 

5583.5 

3.6 

142311851.383 

4626.5 

-3.6 

8*2831815.338 

3481.5 

0.6 

142302247.313 

4145.0 

-0.6 

842313325.326 

5548.0 

3.7 

842302210.322 

4052.3 

0.7 

1.42311353.017 

4330.3 

-0.7 

742311602.003 

5137. S 

3.8 

512832727. 01 8 

4036.0 

-3.8 

742332246.038 

3742.3 

3.3 

842311602.338 

3788.3 

3.3 

842321830.308 

3358.5 

3.3 

ifflS 
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TABLE  4  (Continued] 


TELESEISM 1C 

EVENTS 

EVENT  NfWE 

RANGE 

ERROR 

T 12970915.3*6 

5751.5 

0.0 

1*2300019.051 

5564.5 

-0. 1 

L  *2920526. 308 

51*5.5 

-0.1 

S 12920528. 301 

1181.3 

0.1 

T12320528.306 

5536.0 

-0.1 

T1295 1908. 057 

5533. 5 

0.1 

612373615.039 

5909.5 

0.1 

11*2930019. 053 

1573.5 

0.1 

L *2900013. 051 

1862.0 

-0.2 

612900021.023 

1569.5 

0.2 

11*2951 908. 350 

5725.0 

0.2 

L  *295 1906. 357 

5912.5 

0.3 

(1*2970815.0*3 

6135.5 

-0.9 

L1297 0615. 3*6 

6091.0 

-1.1 

5*2951909.303 

5251. S 

-1.1 

11*2920528.20* 

5799.5 

-3.* 

NUMBER  Of  EVENTS 
HECIAN  =  0. 1 

MEAN  =  0.S 

VARIANCE  =  aS9.S 


=  16 
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noise  poorly.  One  way  to  improve  the  effectiveness  of  the 
detection  method  is  to  extend  the  time  preceding  the  6.4 
second  window  triggering  the  fast  Walsh  detector  to  60 
seconds.  Past  experience  in  picking  the  first  arrival  has 
shown  that  while  the  Pn  arrival  sometimes  precedes  the  Pg 
arrival  triggering  the  Walsh  detector  by  more  than  30 
seconds  all  of  the  Pn  arrivals  are  included  in  the  first 
60  seconds  preceding  the  window  triggering  the  detector. 
Then  only  overlapping  seismic  events  will  cause  the 
performance  of  the  detector  to  deteriorate.  Of  the  400 
event  traces  in  the  data  set  only  four  were  overlapping 
events.  Table  4  lists  the  errors  relative  to  the 


analyst's  picks  in  the  automatic  first  arrival  picks  and 


the  range  of  the  rank  sums  for  all  the  events  and  for  e^ch 
individual  type  of  event;  i.e.,  local,  regional  *  .id 
teleseismic.  The  illustrations  depicting  the  performam « 
of  the  automatic  detector  and  first  arrival  picker  on  15'.. 
event  traces  are  included  in  appendix  D. 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 


The  rank  sun  detector  Is  extremely  robust  and 
excellent  for  distinguishing  a  broad  category  of  seismic 
signals  from  noise.  However,  the  time  preceding  the 
window  triggering  detection  in  the  fast  Walsh  detector 
should  be  increased  to  60  seconds.  To  illustrate  this 
conclusion  let  us  examine  an  extreme  case  where  we  would 
like  the  detector  to  pick  either  Pn  or  Lg  as  the  first 
arrival  for  later  event  association  and  li  cation.  Suppose 
the  fast  Walsh  detector  fails  to  trigger  on  the  Pn  arrival 
from  a  regional  event  13°  from  the  recording  station.  The 
maximum  time  interval  between  Pn  and  Pg  at  13°  is  57 
seconds.  So  if  the  fast  Walsh  detector  triggered  on  Pg, 
the  Pn  arrival  would  still  be  included  in  the  event  trace 
file  if  the  predetection  time  was  increased  to  60  seconds. 
If  the  fast  Walsh  detector  triggered  on  Lg  instead  of  Pg, 
then  we  conclude  the  Pn  arrival  must  be  too  far  below  the 
background  noise  level  for  even  a  human  analyst  to  discern 
otherwise  the  fast  Walsh  detector  would  have  triggered  on 
Pg  instead  of  Lg .  In  this  event  the  human  analyst  would 
pick  Pg  as  the  first  arrival  when  we  would  rather  have  him 
pick  Lg  for  later  association  purposes.  Since  the  time 


the  recording  station  is  3  minutes,  the  automatic  detector  would  select  Lg  as 
the  first  arrival.  The  real  time  fast  Walsh  detector  writes  seismic  events  oc¬ 
curring  closer  together  than  60  seconds  into  a  continouous  event  file.  There 
should  be  no  significant  problem  with  overlapping  events  in  extending  the  pre¬ 
detection  time  to  60  seconds  for  each  event  file.  Then  99%  of  the  automatic 
picks  should  be  correct  within  a  median  error  of  plus  or  minus  0.2  seconds 
(see  Table  5) . 

The  next  step  is  to  develop  a  criteria  to  break  up  the  class  of  signals 
detected  into  local,  regional  and  teleseismic  events  and  then  into  earthquakes 
and  explosions.  Decision  tree  or  stratified  layer  classification  is  designed  to 
take  advantage  of  such  situations  (one  or  two  features;  i.e.  slope,  period  or 
amplitude)  to  improve  efficiency  and,  if  possible  accuracy.  The  decision  tree 
classifier  progresses  through  a  series  of  stages  or  layers;  at  each  layer  certain 
classes  (local,  regional  and  teleseismic)  are  separated  in  the  simplest  manner 
possible.  It  is  flexible  and  permits  different  features  and  classifications  (boun¬ 
daries)  to  be  used  to  separate  different  classes.  (Schowengerdt,  R.A.,  1983). 
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TASL2  5  Errors  (Automatic  Minus  Analyst  Picks  in 
Secondi ).  23  Raises  of  <=h.  Rank  Sums  or  Events 
Shlch  do  not  Violate  the  Initial  Assumptions 


«_l  £vev73  except  -moss  wm 
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HJ_  fVESTS  EXCEPT  THOSE  MHICM 
VIOLATE  THE  IHITIPL  HSSJWTICNS 
Of  THE  WiICCTOa  CE7ECT2« 


EVENT  WWE 

J  *1331418.327 
L12832234.23  3 
Li28322’.2.233 
8*2533127.213 

l  iisaur..  2‘.  3 

T42531923.237 
U2313225.332 
L *25 11 532. 237 
J1i23r.S32.233 
Ti23r.532.233 
L123r.3Sl.3a3 
812911851.351 
Ki23r.3S3.3aa 
T42912I *3. 3lS 
8*2321923. 233 
Si222’.3S3.3So 
M 123321 23. 321 
Si2322'.2S.  313 
S 12322232. 32E 
T 12322233.  313 

512351239. 200 
512352233.31 1 
T 12352337.  216 
T 12373915.  2l6 
L12371752. 213 
512371752.313 
T *2371752. 352 

812371351.200 
T *2371 351.200 
S12S21S22.J26 
K12232223. 237 
L42321353.J31 
H 123333 1 3. 352 
T123213S3. 323 
L 423221 23. 2B9 
T12333013.3S1 
142331113.337 
T  42332123. 206 
312323125.221 
T423S1 113.331 
342332213.215 
T123S1S33.3S7 
L123S2337. 21 5 
KU31 2'.  43. 3*3 
K 12532253. 333 
S *23723 l 5. 233 
T 123122:2. 257 
L 123735:3.  233 
542323523.231 
T 42322523. 208 
L42321533. 335 
K42321S32.323 
T 1233235 1.3*1 


*ff4GE  Essca 


5357.3 

a.  a 

5575.5 

a.  a 

1802.5 

a.  a 

4  *.Z2 .  3 

a. a 

4*25.5 

3.3 

2252.3 

0.2 

5353.3 

3.2 

*522.3 

3.3 

3733.3 

3.3 

5197.3 

9.3 

*626.5 

a.  a 

37*2.3 

a.  a 

*199.3 

a.  a 

5318.3 

a. a 

2953.5 

a.a 

3922.5 

a. a 

5373.3 

3.3 

6i2S.5 

a.a 

5523.3 

a.a 

5122. a 

a.a 

» _ 5 

a.a 

*25*.  3 

a.a 

3633.3 

a.a 

5751.3 

a.a 

5272.3 

a.a 

6171.3 

a.a 

ssas.a 

a.a 

2*32.5 

a.a 

3*13.5 

3.3 

5251. S 

a.: 

2211.5 

3.1 

3358.3 

a.  i 

4573.3 

a.: 

1273.3 

a.: 

5155.3 

•3. 1 

5551.5 

-3. 1 

1923.5 

-a.i 

5737.3 

a.  i 

3591.3 

a.: 

52*9.5 

-a. : 

545*.  3 

3.1 

ecoo  t 

3.1 

3722.5 

a. : 

51 49. 3 

a.  i 

4231.3 

a.i 

SS23.5 

a.  i 

*011.5 

a.i 

51*5.5 

1 

1181.3 

a.t 

5325.2 

-a.i 

5*63.3 

a.i 

5871. a 

a.: 

SS01.3 

0.2 

HPHE 

5*2321*23.333 
T*2a92234.2lS 
S  *253032 1.222 
K42231 *13.253 
L *2222223. 313 
K 42532232. 225 
K12571753, 323 
51251 1533. 223 
L42332254.241 
K423S13aa.2S3 

L 125033 1 3. 3S I 
512311307.325 
T 12311307 . 332 
T12392327.312 
T4220312S.321 
K4232H33.3!  1 
142911951.337 

T 12321 *23. 313 

T42331S3S. 203 

712311353.355 

K4253I823.317 

U232!  823.339 

L123S'.338.357 

T  *2332239 . 323 

K1231 1307. 223 

L42912212.357 

T12321523.J36 

K423 62338. 231 

51233132’. 325 

L123121 *3. 215 

L *23 11307. 223 

L12532217.313 

T 128323 13.  208 

T12310325.229 

8*2891815.238 

>1*2310325.225 

L123U3S3.317 

8*2532213.322 

S *2352327. 313 

7*25322*6.333 

812572815.213 

U2373815.24S 

•*  ^37'  272 . 2 -S 

L *237 1352. 237 

7*2391314.216 

2*2891719.252 

8*2991713.357 

L12392327. 315 

7*2331719.3*7 

8*2831506.3*3 

L 125 32229. 327 
8*2323523.23* 
812392234.323 
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2391. 3 

■•••3  e  2 

57+a.  2 

a.a 

*565.5 

3.2 

5*75.5 

a. 2 

£8*9.2 

•'.2 

C88S. 5 

a.: 

5323. 5 

a.  2 

*523.5 

a.  2 

5531.3 

-3.2 

5725.3 

a.a 

1862. 3 

-a.a 

1369.3 

0.3 

1573.5 

-3.2 

5168.3 

a.  2 

155*.  5 

a.  2 

4353.3 

3.3 

S3S3.S 

a.  3 

29*1.2 

-a. 2 

3898.3 

a.  2 

*632.5 

a.  2 

*773.3 

a.  3 

*112.5 

a.  3 

5912.5 

a.  2 

3712.3 

a.  2 

414*. 5 

a.* 

*4*3.5 

-a.  * 

S7S1.5 

-a.  * 

4723. S 

a.  * 

4187.3 

-a.* 

572a. a 

a.* 

*572.5 

a.* 

41*5.3 

-0.  5 

5532.5 

a.  6 

56*3.3 

a.  6 

3*81.5 

0.5 

55*3.3 

9.7 

1393. a 

-3.7 

1352.3 

3.7 

*396.3 

-a.a 

37*2.3 

3.3 

6125.5 

-3.3 

6091.3 

-l.l 

*8*3.5 

-i .  a 

*562.2 

i.a 

*812.3 

1.2 

*166.5 

1.2 

3826.3 

1.7 

*615.5 

-1.8 

*242.3 

1.9 

32*3.3 

2.7 

2863.3 

3.3 

5799. S 

-3.  * 

*133.3 

-5.1 

100 
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TABLE  5  (Continued) 
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ALL  EVENTS  EXCEPT  THOSE  WHICH 
VIOLATE  THE  INITIAL  ASSUMPTIONS 
OF  THE  RAWC2700  DETECTOR 


LOCO.  EVENTS  EXCEPT  THOSE  WICH 
VIOLATE  THE  INITIAL  ASSUMPTIONS 
Of  the  nswz7 aa  QETECTOR 


EVENT  WANE  RANGE  ERROR 

EVENT  NAME 

RANGE 

ERROR 

342437304.012  5436. 0  -6.3 

L428 32304. 039 

5576.5 

0.0 

S42891719.025  4805.5  7.4 

T423121 40. 815 

5518.5 

0.0 

S 42962038. 81 4 

4264.5 

0.0 

UMBER  OF  EVENTS  =  108 

T42362337. 016 

3655.0 

0.0 

8ETIAN  s  0.2 

L423717S3.018 

5272.0 

0.0 

MEAN  =  0.5 

542371753.013 

6174.0 

0.0 

VARIANCE  =  1.3 

'  T 42371752. 052 

5505.0 

0.0 

842371551.000 

2403.5 

0.0 

T 42371 9S 1.000 

3419.5 

0.0 

842902055. 038 

4234.0 

0.1 

L42362037.016 

3722.5 

0.1 

842912140.040 

5149.0 

0.1 

L423020S4.044 

5504.0 

-0.2 

T 42892304. 016 

5740.0 

0.2 

T42302054.  044 

5904.3 

0.2 

842971753. 230 

5339.5 

0.2 

L423121 40.015 

5730.0 

0.4 

842562038.001 

4733.5 

0.4 

L 4237 1352. 007 

4583.0 

1.2 

842892304.028 

4193.5 

-5.1 

542832304.012 

5436.0 

-6,3 

NUMBER  OF  EVENTS  = 
MEDIAN  *  0. 1 

BEAN  -  0.7 

VARIANCE  =  11.5 

21 
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TABLE  5  (Continued) 


REGIONAL  EVENTS  EXCEPT  THOSE  MH1CH 
VIOLATE  THE  INITIAL  ASSUMPTIONS 
Of  THE  RANK 2700  DETECTOR 


REGIONAL  EVENTS  EXCEPT  THOSE  NHtOI 
VIOLATE  THE  INITIAL  ASSUMPTIONS 
Of  THE  SANK 2722  DETECTOR 


EVENT  NAME 

3*2031*18.037 

L  *2832313- 023 

11*2300137.013 

L*2S01 821.01 6 

T 4230 1820. 237 

L42S1032S.332 

L42S11602.007 

H*23UB02.038 

T *231 1602. 023 

L423118S1.003 

11*2311851. 0S1 

H423U353.008 

H4232 1330. 008 

342321853.056 

8*2332123.231 

S42332123.018 

5*2332303.035 

T42932303. 013 

L42S21S33. 036 

0*2321533.023 

3*2321533. 026 

S42300136. 031 

1*2321853. 031 

L42891418.037 

T  *232 '.833. 22  8 

L *2332123 . 009 

3*2832313.016 

8*2832328.007 

T42332123. 006 

7*2831418. 034 

T42312213.0S7 

8*2831*18.053 

1*2332303.013 

8*2332303.035 

5*2321*30.030 

3*2311603.000 

T42831305.088 

T *291 1851. 007 

0*2901820.0*7 

1*2311858.056 

T42832327.212 

0*2321*30.01* 

T *2302203. 030 

1*2321123.018 

3*2311607.026 

T 423 11 807. 232 

1*2300136.  Z21 

L *232 1823. 039 

5*2301821.025 

T42321S33.036 

L42312213.0S7 

L423U807.029 

0*2311807.023 


RANGE  ERROR 


5857.0 

0.0 

*802. S 

9.0 

*238.0 

0.0 

4*05.5 

0.0 

3292.0 

0.0 

5958.3 

0.8 

*522.0 

0.8 

3788.0 

0.0 

5197. 5 

0.0 

*628. S 

0.8 

37*2.0 

0.0 

*133.3 

0.0 

3958.5 

0.0 

3923.5 

0.0 

5073.0 

0.0 

6385. S 

0.0 

5628. S 

0.0 

5323.3 

0.0 

5463.0 

0.1 

5871.3 

0.1 

5261.5 

0.1 

3S81.0 

0.1 

3338.3 

0.1 

4920.5 

-0.1 

4275.0 

0.1 

5*95.5 

-0.7 

5*54.0 

0.1 

3311.5 

0.1 

5797.0 

0.1 

53*9.5 

-0.1 

*014.5 

0.1 

5*75.5 

0.2 

58*8.0 

0.2 

5883.5 

0.2 

3994.0 

-0.2 

*508.5 

0.2 

3880.0 

0.3 

5059.5 

0.3 

4770.0 

0.3 

*692.5 

0.3 

5166.0 

0.3 

*350.0 

0.3 

3713.3 

0.3 

33*4.0 

-0.3 

*363.0 

0.3 

*579.5 

-0.3 

*554.5 

0.3 

*112.5 

0.3 

*187.0 

-0.4 

5751.5 

-0.  * 

4*43.5 

-0.4 

*572.5 

0.4 

*1*4.5 

0.4 

EVENT  NAME 

T 42310325. 023 

T*28 92313. 008 

L42302247.013 

8*2831815.038 

0*2302270. 322 

8*2310325.026 

L429U9S9.017 

3*2832327.018 

7*2302246.038 

T4292 1823. 038 

Li2831719.050 

1*2831814.0*6 

042831719.057 

U2832327.016 

1*2891719.0*7 

842831306.0*0 

L *2302209. 027 

S42891719.025 

NUMBER  OF  EVENTS 
8EDIAN  *  0.2 

MEAN  =  0.S 

VARIANCE  =  4.5 


RANGE 

ERROR 

56*3.0 

0.6 

5583.5 

0.6 

*1*5.0 

-3.6 

3*81.5 

0.6 

*052. a 

0.7 

55*0.0 

0.7 

*330.8 

-0.7 

4096.0 

-0.8 

37*2.0 

0.9 

4849.5 

-1.2 

4166.5 

1.2 

*612.0 

1.2 

3826.0 

1.7 

4615.5 

-1.8 

42*2-5 

1.9 

33*9.5 

2.7 

2660.0 

3.0 

*805.5 

7.4 

71 

.8 
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TABLE  5  (Continued) 
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TILESEISMIC  EVENTS  EXCEPT  THOSE  WHICH 
VIOLATE  THE  INITIAL  ASSUMPTIONS 


EVENT  NAME 

342351339- 333 
T42370815.046 
1.42928528.308 
$42322528.331 
T42928S28 . 308 
M42S800I9. 353 
T423SI908.3S7 
S 42320815. 333 
T 42900019. 3S1 
S42333321 . 323 
1442351308.350 
L4233 3319. 051 
L4235 1338. 357 
M4297381S.043 
L4297081S.346 
H42323S28. 334 


RANGE 

ERROR 

5251.5 

8.3 

5751.5 

8.3 

5145.5 

-0.1 

4184.0 

8.1 

5535.3 

-3.1 

4579.5 

8.1 

5539.5 

0.1 

5333.5 

0.1 

5564.5 

-0.1 

4563.5 

3.2 

5725.3 

0.2 

4862.3 

-3.2 

5912.5 

3.3 

6135.5 

-0.9 

6331.3 

-1.1 

5739.5 

-3.4 

NUMBER  Of  EVENTS  -  16 
MEDIAN  =  8.1 

MEAN  =  8.4 

VARIANCE  =  22.3 
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APPENDIX  A 


FAST  WALSH  TRANSFORM  DETECTOR 


The  Walsh  functions  are  an  ordered  set  of  rectangular 
waveforms  whose  amplitudes  take  the  values  +1  or  -1.  They 
are  arranged  in  order  of  increasing  number  of  zero- 
crossings  per  time  interval,  so  called  sequency  order. 

Sequency  is  defined  as  one  half  the  number  of  zero- 
crossings  per  interval .  The  sequency  of  Walsh  transforms 
is  analogous  to  the  frequency  of  sines  and  cosines  used  in 
Fourier  analysis. 

Any  time  series  can  be  expressed  as  the  weighted  sum 
of  a  series  of  Walsh  functions.  The  Walsh  transform 
provides  the  coefficients  for  the  summation  in  the  same 
way  as  done  by  a  Fourier  transform.  The  advantage  of  the 
fast  Walsh  transform  is  that  it  can  be  computed  much 
faster  than  the  Fourier  transform  because  it  involves  only 
integer  addition  and  shifting.  The  macro  subroutine  for 
the  fast  Walsh  transform  used  in  the  real  time  event 
detector  is  listed  in  appendix  C. 

One  complication  of  the  fast  Walsh  transform  is  that 
it  does  not  produce  a  spectrum  in  sequency  order. 
Instead,  it  produces  a  Paley  ordered  spectrum  (sometimes 
called  natural  order)  in  which  the  subscripts  have  been 
bit-reversed.  The  seismic  event  detection  program  uses 
lookup  table  A.l  to  convert  a  sequency  subscript  from  the 
passband  into  a  bit-reversed  Paley  subscript  in  the 
transformed  spectrum.  ( GSETT  Report,  Paul  Golden,  1985) 


The  time  series  data  in  the  detection  window  is 
transformed  using  the  fast  Walsh  transform.  The  absolute 
values  of  the  weighted  coefficients  are  summed  across  the 
pass  band.  The  sum  is  then  normalized  by  shifting  to  the 
value  it  would  have  had  if  the  weights  had  all  been  1. 
This  normalized  sum  is  the  magnitude  which  is  compared  to 
the  threshold  for  the  channel.  A  magnitude  below  the 
threshold  triggers  the  detection. 

When  detection  occurs  the  signal  is  written  to  an 
event  file.  This  event  file  is  later  graphically 
displayed  to  enable  a  seismic  analyst  to  discriminate 
between  an  "event"  signal  and  "glitches"  or  "noise"  and  to 
pick,  the  event  parameters.  (Goforth  and  Herrin,  1981) 
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appendix  B 

Subroutines 


i  4  A  ».  r.  y.  jT.  i 


L «  Sjl  ^  •>.  V  -  *f  a.  Km  n*  m.  a  »  R  -  A  .  r  «  *s  .  pi  «  r  r  _  »  _  r  -  V^_  p^  *  «  pi  _  mV  j«  -  rf\_  _ 
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LIST  OF  SUBROUTINES 


AR_SPEC7RA  -  COMPUTES  BRIPI  SPECTRAL  ESTIMATES 

RKAIKE  -  COMPUTES  *  YULE-WfCKER  COEFFICIENTS 

COV  -  COMPUTES  AUTOCORRELATIONS  FOR  A  TIME  SERIES 

DETECTOR  -  DETECTOR  USED  TO  TEST  N0NPARAHE7R I C  TESTS 

IFWM64  -  INTEGER  FAST  W«_SH  TRANSFORM  IN  MACRO  FOR  64  POINT  WINDOW 

NEIGMBCR  -  COMBINES  SEOMEHTS  WITH  NEAREST  NEIGHBOR  DECISION  RULE 

POP  -  COMPUTES  DENSITY  DISTRIBUTION  MlSTDR GRAMS 

Polr_A-BURG  -  COMPUTES  AR  BURG  COEFFICIENTS 

RANKFREO  -  COMPUTES  R»*S  AND  FREQUENCIES  FOR  N  OBSERVATIONS 

RANKTEST  -  COMPUTES  THE  RANK  SUM  St'“  TST  FOR  N*M  OBSERVATIONS 

RANK2733  -  OETECTOR  BASED  ON  A  MODIFIED  RANK  SUM  TEST 

RUNEXP  -  COMPUTES  THE  RUN  TEST  FOR  N«M  O&SKVATIONS 

SIGN  EXP  -  COMPUTES  THE  SIGN  TEST  FOR  N  OBSERVATIONS 

SLOPE  -  COMPUTES  MODIFIED  SLOPE  FOR  A  TIME  SERIES 
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c- 

c- 

c 


subroutine  sLSPECTRAtFMAx.FMiH.coEF.HCQEF. 

^pudkuu  PF,PFLOO.NF,VAR«ARMAX> 

TNl.  •bhnwlln.  oo-^U.  U.  «  ~ll“^ 

tflMd  bM  K.8.  Prl-Uy  In  Sp~lrol  Ti**  5*r,~ 

(paparboekl #  p*  681#  1981* 


I  FRESHEN  C7  MAXI** 

I  FRESHEN  CT  MINIMUM 
If*  COEFFICIENTS 
1WL1TU0E  OF  AR-SPECTRA  fit  F 
1  AMPLITUDE  IN  d8  OF  AR-SPECTRA  AT  F 
(VARIANCE  OF  NHITE  NOISd.  OF  .  R»-~-- 
(HUMBER  OF  AR  COEFFICIENTS 


REAL  FMAX 
REAL  FX1N 
REAL  COEFtNCOEFI 
REAL  PFI2S11 
REAL  PFU30125U 
RES.  VAR 
INTEGER  NCOEF 
COMPLEX**  OFlff. CONST, SUN 
COMPLEX**  •rwyiJ/fB.B'S. 263185387181/ 

^^,28318538718/  rews  w  «*-SPECTRA 


:_  raPUTE  SPECTRAL  ESTIMATE  FOR  EACH  FRE8UEJCT 


l COMPUTE  FREQUENCT  INCREMENT 
lIK-lllh  FRESHEN  Cl 


F  =  FMIN 

SRRTE  =  FLOAT  INF -11  /1FMAX-FH1N1 
DO  t  -  l.NF 

F  =  (FLOAT (K-1)/SRATE)  *  FHIN 
SUM  =  OPJttl.  8.8.81 
IF  (NCOEF. GE.11  TFEN 

00  1*1  »NC2EF  ICOW-EX  FREOUENCT 

END  00 
END  IF 

const  -  SUM»CCNXISUMI 
PFtri  *  l. /RES. (CONST) 

PF IK)  *  PF tKl *2. *VAR 
pfLCGtKl  =  28.*AL0G18(PFIK1 ! 

(yiKAX  -  MBXIPFUQGIKI.ARMAXI 
END  00  IK 
RETURN 
END 


IVAR  =  VARIANCE  OF  NOISE  DATA 
I  AMPLITUDE  IN  &  OF  SPECTRA 
(MAXIMUM  SPECTRAL  AMPLITUDE 
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m 
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SUBROUTINE  *flIKE(X.N.M,RCf.l«..COEF) 

C 

c 

C—  Modified  fro*  Qroy-Nood»ard  RRMfl  •p*otral  grtalyala 

C—  paoKOQ*  for  Lha  IBM  PC  by  Kothloon  ft.  flidoo  1-30-86. 

C—  Thl«  •uorouUn*  aeopulM  IKa  Tulr-WolWor  ooaff Iclania 
C—  a*  dooerlbod  by  Kay  and  Mar-pi*.  (1382). 

C 

REEL  Kill  ISERIES 

REPL  FPE138)  IHORK  fWRflf  FOR  F1NBL  PREDICTION  ERROR  IFPE) 


REPL  RCflll  I RUTOCORR ELATION 


RE8L  ROB.  381  I  WORK  RRRPT 
REPL  COE?  1 1 )  ITULE-nSJCER  COEFFICIENTS 
INTEGER  N  (LENGTH  OF  X  SERIES 
INTEGER  M  IMftXJHUM  ORDER  OF  FPE 
INTEGER  >f_  I  ORDER  SELECTED  BT  FPE  • 

C 

HRG=188 
lCCEf=t 
3UH=8.0 
DO  SI  1=1. N 
SUM  =  SUM  *  XII) 

51  CONTINUE 
XBW=5UN/N 
XSQ^.0 

00  S2  1=1. N 

XU)  =  XII)  -  XB« 

xss  =  xsa  *  xiiiaxin 

52  CONTINUE 
PV  =  XSO/N 

C  CPLCULPTICN  OF  RR  COEfS.  USING  TULE-WRLXER  EQUATIONS 
RU.l)  =  RCFI2) 

FPE  111  =  It .0-tftll. llaflll.il )  la(N*2)/IN-2) 

Ml  =  M-l 
00  2  1=1. Ml 
SN  =  8.8 
SO  s  g.8 
00  3  Jsl.I 

3N  =  SN  ♦  flll.JlaflCf  U*2-J) 

30  =  SO  ♦  RU.JlaflCf IJtll 

1  CONTINUE 

B  11*1. 1*11  =  IRCF 11*2)  -  5N!/ 1 1.0- SO) 

FPElUl)  *  FPE(IU(N*I*2)»(N-I-I)all.fr-R1I*1.I»11»R(I*1.I*1)  )  / 
♦  I  (N-I-21alN*I*l!  1 
00  4  J=l.I 

fl(M.J)  =  R(I.J)  -  flU»l.I*l)»fifX.I-J»l) 

4  CONTINUE 

2  CONTINUE 

JFIICOEF.  ES.  8)  GO  TO  IS 
16  CONTINUE 
15  CONTINUE 
1=1 
J=2 

7  IF(FPEU)  .LE.  FPEIJ1)  DO  TO  S 
I  =  J 
J  =  J  *  1 
DO  TO  6 

5  J  =  J  *  l 

4  IFU  ,GT.  Ml  DO  TO  6 
DO  TO  7 

6  *.=! 

DO  65  1=1, ML 
(5  COEfUlsRIML.I) 

RETURN 

END 
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SU8R0UTINE  COVIX.L.XB,  W.Q.N) 

C 

C —  Th l a  wiBroutlrwa  calculator  N  autoeorrrlatlon*.  tb» 

C—  mbbI*  won,  <rid  ihm  varlano*  *qual  to  ihr  loro  log  aulo- 

C—  varlano*  for  lh*  Input  t!*«  oorloa.  nil).  1=1, L. 

C—  Tcfc*n  fro*  lh*  Grog  Woodoord  ARMS  apaelral  Mil  sail  on 
C —  pack  ago.  1365. 

C 

c—  x=series  -iwnrr 

C—  L^ENGDl  OF  SERIES  -  INPUT 

C—  XB=ME PH  -  OUTPUT 

C—  W=V«IANC*  -  OUTPUT 

C—  C=flUTC~RRE_ATICNS  -  OUTPUT 

C—  N=NUM2ES  OF  AUTOCORRELATIONS  TO  CALCULATE  -INPUT 

C—  NOTE :  0  ( 1 1  =AUTOCCRREJ?TIONS  ST  LAG  1-1 

C 

REAL  XU) 

RES.  Oil) 

c 

C-  COMPUTE  TIC  SAMPLE  TOR  XII).  1=1.1 

c 

XB=0.B 
DO  25  1=1. L 
25  XB=Xfl*X!I! 

XB=XB/FU»T(LI 

C 


C—  COMPUTE  DC  SUTOCOVSR I ANCES  ST  LAGS  =  8.  IN-1) 

C 

DO  30  1=1. N 
0(11=0.0 
NK=L-I*t 
DO  40  JSl.MK 

40  D(I)  =  tX(Jl-X0)«(XlI*J-l)-XB)«CU) 

30  CONTINUE 
C 

C—  SAMPLE  VARIANCE  EQUALS  DC  AUT0C0YRR I RNCE  AT  LAG  =  0. 

c 

W=0111 

c 

C—  NORMALIZE  DC  SUTCOOVRR I ANCE5  TO  08TBIN  DC  SUTOCBRRELBT t ON 
C 

00  50  1=1. N 
50  Din=OtU/W 

return 

END 
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SUBROUTINE  DETECTOR  INSAMPLES.N. MEANER .INC* 

*  Tflf  spa  FN .PATTERN. 1ST. SN. RANK! 

-  SUBROUTINE  TO  PERfORH  TWO  SAMPLE  SION  TEST 

—  HOODT.CRATBIU.  4  DOES 

REAL  TEST 

REAL  PATTERN!!)  ll.=»*.2.** 

REAL  TRCtl)  (TRACE  CONTAINING  SIGNAL 
RES4.  SECLENU1  (SAMPLE  LENGTH  Of  PATTERN 

fccai  UCOU 

gE,qL  SN  IWW  SIGNAL  TO  NOISE  RATIO  FOR  TRACE 
•  yjT?  yFj?  UDTUlP*  c  - 

INTEGER  INC  I  INCREMENT  WINDOW  HOVES  DOWN  TRACE 
TWTFCT2  1ST  V I STP^RT  IN  HlNQQN 

INTEGER  N  I  NUMBER  Of  SAMPLES  IN  NOISE  PART  OF  TEST 

INTEGER  K  IACCEPT  OR  REJECT  CRITERIA 

REa.  Tisaa.3!  1WORK  ARRAT  for  ordering  samples 

RE*.  RANK  111  ISUMATION  Of  RANKS  Of  SIGNAL 
LOGICAL  ACCEPT/. FALSE./ 

C _  COMPUTE  RANK  SUN  TEST  OVER  TRACE 

C 

M2  s  2*N 
SN  *  MEAN 
ENMKX  =  8.8 
DO  Jst.N 
L  «  M 

PATTERN  (NSAMPLES-LI  *1.8 
SEOLtN INSAMPLE S-Ll  *  t-8 
PA  TERN  1 1 )  *  t-8 
SEGLE.Nl I)  *  1.8 
711.11  *  ABStTRCtJ»2-l)l 
Tll*N.ll  *  ABSlTRCl J»21 1 
EKPAX.  =  MAXIT ( 1. 1)  .ENMAX) 

ENHAX  =  MAX  IT  I  I«N.  11  /ENMAX) 

END  DO 

B*«X  =  EJtlAX*  1 . 35 
IR  =  8 

DO  1=1, INSAMPLE5-N! . INC 
IR  =  IR  *  I 
00  Jsl.N 
jj  =  J  ♦  (I-U 
TU*N.l)  =  BBSITRCtJJO  1 

END  DO  1J*WN  _  _ 

CrtJ.  RAMTTESTlN.Tll.U  ,T1N*1»1)  ,K.TY. ACCEPT) 
SN  s  MAXITT.SN1 

rahkiiri  =  tt 

IF  (ACCEPT!  THEN 
00  L=i > II*IINC-1)  1 
PATTERN l D  =  2.8  1RE.ECT  H8 
SEOLENiLl  =  1.8 
END  00 
ELSE 

DO  L*I .  1 1* ( INC-11  1 
PATTERNIL1  =  1.8  IACCEPT  H0 
SE0LEN1L!  *  1.8 
END  00 
END  IF 
O©  DO  II 

SN  =  SN/RANKIU  IMAX  SIGNAL  RAF*  /NOISE  RANK 

RETURN 

END 
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.THU 


IFWM64 


i  evrTN  £4  INTEGER  FR5 T  WALSH  TRANSFORM. 
t»L£MENTS  FDRTRW  VERSION  IN  1FVF64.FOR  . 

S  NflMU  IN  COMMENTS  ARE 

ROM  THE  FORTRAN  VERSION, 
av  nWIGKT  OBJRINC  ON  22-N0V-62. 

Lk!*i>sep-<.4  to  NondL  I-U^  ~-f,~ 

IKTE0ER«2  X184I 
CPLL  IFNH641X! 

PRODUCES  TRANSFORM 

SIT-REVERSED  DTfCIC  INATURrul  ORDER. 

Q f  elements 

512=64  ,M  =  2.6 

PNR=6 

L=SIZ/2 


.PSECT 

.ENTRT 

REGISTER  USAGE; 
IS 
R1 
R2 
R3 
R4 
RS 
R6 
R7 


IFWN64 

IFVIM64. 


1 

L»2 

R 

POINTS  TO  X  IP 1 
POINTS  TO  XtPPL) 
J 

X»2 
XP 


.MfRO/RS.R*.®'*6'*7' 

ITERATION  CONTROL 
L  AS  WORD  INDEX 


K  AS  NORO  INOEX 


10$‘ 


20*; 

3 


35$' 
40$ 1 


45$  • 
50$  t 


MOVZSL 

HOVL 

H0V23L 

NOVL 

AOOL3 

M0V2BL 

M0V2BL 

MOVW 

PD0W3 

BVC 

BLSS 

MOVW 

era 

MOVW 

SUBW3 

BVC 

BLSS 

MOVW 

era 

MOVW 
ADDL2 


aPWR,R0 

•<L»2>.R1 

•UR2 

4lflPl.RO 

R3.RI/R4 

al.RS 

•2.R6 

IR3I/R7 

R7.lR4l.lR3I* 

40$ 

35$ 

«.x60ei.-2iRO! 

40$ 

..X7FFF,-2(R3! 
IR41  .R7.IR41* 
50$ 

45$ 

m*  xaaei.-oiR*1 

5jj5 

a/X7FFF/-2lR4) 
•2,R6 


CNPL 

R6.R1 

BUO 

30$ 

A0CL2 

RI.R3 

flDOU 

RI.R4 

iNa 

RS 

CNPL 

R5.R2 

BLES 

20$ 

ASHL 

a-l.Rl.Rl 

ASH. 

al,R2.R2 

S08GTR 

R0. 10$ 

RET 

.END 

<L»2 


/ITERATION  COUNT 


/R 


/AOCRtXIP)).  P*1 
iRODRIXIPPLIl •  PPL=L»1 
/j=: 

,-&Zl  (WORD  INDEX! 

,xp=xipi 

,XtPl=XP*X(PPLI.  P=P*1 
/no  ovarf lo*.  OX 

».lgn  aNong-d  to  nagatlv. 

/ Iqrgaat  nagotiv**!  lao  1AB51 

/uaa  la-gaat  poalllva 
,XIPPL!=XP-X(PPLJ.  PPL=P?L 
/no  overflow.  OX 
talgr  enanvad  to  nagollvo 

/ Iqrgaat  nagotl''**!  l*°  1Rft~  ‘  • 

/uaa  Iqrgaat  poaltlva 
/2X=2X*2 

zIF1X.LE.L1 

/P=P*L 

iPPL=PPL*L 

/J=J*1 

/IFIXLE.R1 

iL^yo 

»R=*»2 


115 


\  .  N  *4- 


— '  ™>wwwiinreopgp WHII  U.'VJH  W  mur 


c 

c— 

c— 

c— 

c— 

c— 

c— 

c— 

c— 

c 

c— 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 
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c 

c— 

c 


c 

c— 
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SUBROUTINE  NEIGHBOR  !MA&EN.SEO.NSEGMENTS.NHIK) 


Tha  hwmI  rwtqrvber  daoialon  rula  ftrat  ooapcraa  tha 
faatira  vector-  far  acch  tippinl  of  an  array  of  obaarvot  Iona. 

Each  MgMnt  la  toRAod  »lth  a  3  or  l  to  Indloata  Ita  affinity 
to  althar  tha  pravlaua  or  na*t  aagaant  raaaacttvaly.  Than 
lha  aagaanta  ahtah  hova  an  affinity  for  aocn  othar  «-a  ooaolnab 
and  thalr  faaturaa  ora  ovaroyad.  Thta  sentlnuaa  until  tha 
nuadar  of  aaqaanta  bacoaaa  laaa  than  or  ooual  to  tha  alnlaua 
nuadar  of  aayaanta  daalrad. 

Coda  vrltton  by  Kathlaan  A,  ft  I  oar ,  Morcn  14.1984. 

SEGU.J1  -  IWMJT  VECTOR  FOR  ITH  SEGMENT,  J*i.8 
SEGII.il  a  AMPLITUDE 
3EGU.2J  *  ESTIMATED  PERIOD 
SEGU.31  a  NUMBER  Of  OBSERVATIONS  IN  ITH  SEGMENT 
SE0U.4I  a  ABSOLUTE  DIFFERENCE  BETWEEN 
1-1  AND  I  SEGMENT  AMPLITUDES 

SEGlt.SI  a  ABSOLUTE  DIFFERENCE  BETWEEN 
1-1  AND  I  SEGMENT  PERIODS 

SEGU.61  a  ABSOLUTE  DIFFERENCE  BET€EN 
1*1  AND  I  SEGMENT  AMPLITUDES 

SE0U.71  a  ABSOLUTE  DIFFERENCE  BETWEEN 
I«1  ANC  I  SEGMENT  PERIODS 
SEGM.Bl  =  8  INDICATES  WFINITY  HITH  NEXT  SEGMENT 

1  INDICATES  AFFINITY  WITH  PREVIOUS  SEGMENT 
NSEGMENTSMENT5  a  NUMBER  Of  SEGMENTS  IN  TRACE 
NMIN  a  MINIMUM  NUMBER  OF  SEGMENTS  IN  TRACE 

RES-  SEGIMAXLEN.Sl 

INTEGER  NSEGMENTS  I  TOTAL  NUMBER  OF  SEGMENTS 
INTEGER  NMIN  (MINIMUM  NUMBER  OF  SEGMENTS 

DETERMINE  AFFINITY  OF  EACH  SEGMENT 

00  HHILEINSEGKNTS.  GT.NMIN! 

SEGI1.8I  a  8.  1 1ST  5AMPU  AFFINITY  H/NEXT  SEGMENT 

SEG1NSEWENTS.81  a  1.8  I  LAST  SAMPLE  AFFINITY  W/PREVIOUS 

DO  Ia2. NSEGMENTS- 1 

SEGtI.4)  =  ABStSEOtl-l.ll-SEOU.t!)  IPREV  AMPLT  OIF 

SEGtl.Sl  a  ABSlSEGU-l.21-SEGII.211  IPREV  PER  DIF 

SEGlt.SI  a  AflSISEGUal.ll-SECIl.il!  (NEXT  AMPLY  DIF 

SEOU.71  a  ABS(S£Oll»1.2)-SEOII.2) 1  (NEXT  PER  DIF 

IF  t51GU,4)  .OT.SEGU  ,61 1  THEN 
IF  (SEGII.Sl.GT.SE0a.71l  THEN 
SEGI I .8!  S  8.8  IPERICD  4  AMPLITUDE  AFFINITY  WfHEXT  SEG 
ELSE 

SEGtl.Sl  =  8.8  I  AMPLITUDE  AFFINITY  H/NEXT  SEG 

END  IF 
ELSE 

IF  tSEGU.51.LE.SEGtt.7U  THEN 
SEG!  1.81  =  1.0  I  PERIOD  t  AMPLITUDE  AFFINITY  K/PREV  SEG 
ELSE 

SEGtl.Sl  a  1.0  I  AMPLITUDE  WFINITT  W/PREV  SEG 

END  IF 
END  IF 
END  DO 

COMBINE  SEGMENTS  HITH  AN  AFFINITY  FOR  EACH  OTHER 

Raj  UNIT  in.  MIMBER  OF  NEW  COMBINED  SEGMENTS 
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JF  15EG12.S'  .LT.8.S!  T>€N 
K  =  K  .  1 

END  If 

qq  i  *  2,KSEO<e^S 
l  s  t  *  1 

(SEGtI-M! 


SEG1K.1I  =  - 

SEG1K.2)  *  tSEGH-1.21 
SEGIK.31  =  tSEGU-t.3) 


,  SEGII.V"^- 
,  SEGII.2! 1/2. 
.  SE51I.3)) 


ELSE 

K  =  K  ♦  1 

SEG1K.U  =  SE21I.H 
SEG1K.2)  «  SEC! 1.2) 

SEGIK.31  =  SEGU.31 
END  If 

®i5  tSEGl  1*1.81  .LT.B.S1  THEN 
K  =  K*l 

SEG1K.11  =  SEGll.ll 
SEG1K.21  =  5EG11.2) 

SEGIK.31  =  SEG11.31 
END  If 
END  If 

END  00  -  »up\j 

IF  tSEGU-l»8)  .LT.8.S1  .MEN 

*  =  i«ril-l.U  ♦  SEGlMll^- 


FOR  NEN  SEG 

(fWLITUCE 

(PERIOD 

(LENGTH 


(AMPLITUDE 

(PERIOD 

(LENGTH 


(AMPLITUDE 

(PERIOD 

(LENGTH 


(C0H3INE  FEATURES  FD»  st0 


SEG1K.21  =  - 

SEGIK.31  =  (SEGlI-1'2) 

END  If 

B^SDo1»I^LEtNSEaCNTS- CT'MHINI 

RETURN 

END 


SEG1 1.2) 1T2. 
►  SEG1 I .31 1 


(AMPLITUDE 

(PERIOD 

(LENGTH 


G 

v! 

si 

< 

% 

k. 

» 

0 

s 

> 

& 


4 
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y^V.V.VW  A 


PROGRfW  POP 


C 

C—  TH I «  progrea  ealculot**  th*  d*r\»lly  distribution 

C—  h  I  •  cogr » •  for  Um  book  ground  not**.  flroi  IB 00  point*/ 

C—  and  tho  ilyial  *  nal«o>  oooord  1908  point*/  of  th*  **l**le 

C—  *v*nt  trse*. 

C 

PSMfTEI  IMAXLBW900) 

PARAMETER  (MAXNTRC=41 
PfKAMETER  IIBUFSIZ=256) 

PARAMETER  ININTERV*1000! 

CHARACTER»80  BLANKS 
CHAAACTER»1S  FLNAME  IHAXNTRC! 

CHARAC7E8»12  TRCYPE 
CHfWAC7ER«B0  TCHHE 

LOGICAL  SAMESCALE/. FALSE./  IADJUST  ALL  TRACES  TO  SAME  SCALE  IF  .TRUE. 
LOGICAL  CLASS 

REAL  CLASS1NT  ICLASS  INTERVAL  FOR  PRQ9. OENS. CALCULATION 

REAL  TRACE (MAXLEN.MAXNTRC)  ITRACES  FOR  SEISMIC  EVENTS 

REAL  P0N0ISEININTERV,MAXNTRC1  IPROBAfllLITT  DENSITY  DIST.  NOISE 

REAL  PDSIGNALININTERV,MAXNTRC1  IPROBAfllLITT  DENSITY  DIST.  SIGNAL*NOISE 

REAL  SAMPLERATE/  40 .  /  (SAMPLES  PER  SECOND 

REAL  SCALETRC IHAXNTRC 1/1.0. 1 .0, 1. 8, 1 . 8/  ISCALE  TRACES  FOR  PLOTTING 

REAL  FSVAUJE/0.0/  IFIRST  SAMPLE  VALUE  IN  SECONDS 

REAL  XIWLOT/24./  IX  LENGTH  OF  PLOT  IN  INCHES 

REAL  YINPLCT/20./  IT  LENGTH  Of  PLOT  IN  INCHES 

REAL  AMIN 

REAL  »WX 

REAL  GMIN/100./ 

REAL  GMAX/-100. / 

REAL  FHIN/0.0/  IFREDUENCY  MINIMUM 

REAl.  FHAX/0. 0/  IFREDUENCY  MAXIMUM 

REAL  AHT/0. 21/  (SIZE  OF  CHARACTERS  IN  INCHES 

REAL  XIN  IX-AXIS  LENGTH 

REAL  YIN  IY-AXIS  LENGTH 

REAL  EVTXIN/24./  IEVENT  X  AXIS  IN  INCHES 

REAL  EYTT IN/12./  IEVENT  Y  AXIS  IN  INCHES 

REAL  PRC8XIN/S./  I PROBABILITY  DENS  X  AXIS  IN  INCHES 

REAL  PROSY  IN/S./  IPROBAfllLITT  OENS  Y  AXIS  IN  INCHES 

REAL  EYTMRX/-I0.0/  IEVENT  MAXIMUM  PLOTTED 

REAL  EVTMIN/I0.8/  IEVENT  MINIMUM  PLOTTED 

REAL  FREOMIN/0.0/  IRELATIVE  FREDUENCY  MINIMUM  PLOTTED 

REAL  FREDMAX/I.B/  I  RELATIVE  FREDUENCY  WIXIMUM  PLOTTED 

REAL  IXCFF/1.0/ 

REAL  ITOFF/1.0/ 

REAL  FAC/1.0/  I  FACTOR  TO  SCALE  PLOT 

REAL  XORIGIN/0. 0/  I  INITIAL  PLOT  ORIGIN 

REAL  YORI5IN/0. 0/  I  INITIAL  PLOT  ORIGIN 

REAL  AFACTOR/0. 0/  IY-AXIS  TRACES  OVERLAPP  IF.LE.0. 

INTEGER  ITIMEIfll  I  ARRAY  TO  COMPUTE  CORRECT  TIME  FOR  1ST  SAMPLE 
INTEGER  NPCFSMPS/1000/  INUMflER  OF  SAMPLES  IN  PROBABILITY  OISTR. 

INTEGER  HSWLES/2490/  IHJMSER  OF  SWPLES  IN  EACH  TRACE 

INTEGER  NCLRS5INT/100/  INUMBER  OF  CLASS  INTERVALS 

INTEGER  NPLOTS/I/  UNITIALIZE  THE  NUMBER  OF  PLOTS  PLOTTED  TO  1 

INTEGER  STRT5AMP IMRXNTRC1 /I. 1.1.1/  IFIRST  SAMPLE  TO  BE  PLOTTED 

INTEGER  XLABEL 1 20 1  I  LABEL  FOR  X-AXIS 

INTEGER  TLABELI20I  I  LABS.  FOR  Y-AXIS 

INTEGER  PLABEL1201  I  LAB  EL  FOR  PLOT 

INTEGER  NCHPRSX/0/ 

INTEGER  NCHAR5Y/0/ 

INTEGER  NCHPRSP/0/ 

INTEGER  NTRACES/4/  I  HUMBER  OF  TRACES 

lnl*g*r  Ibuff  lbuff*r  *lx*  for  »tandard  ooloo»o 

lnt«g*r  no  lnu*6*r  ct^iaard  colcc*e 

lnt*g*r  I gu  I  logical  unit  nu»o*c  for  PLDTSIN  naw*ll*t 
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*■»  'V*  *  *tri  Kj*  ilt  X. 


rr.  v  v : v  ww  ?y?  * 


-  CONTINUED  - 

PR0GRW1  PDF 


o 


C 

C 


c 

c— 

c 


c 

c 

c— 

c 


Integer  11/22/  llogleal  tnlt  re. mbmr  for  plotting  device  Incut 

Integer  lo/21/  llogleal  jilt  nt»ber  for  platting  dovloe  output 

Integer  le/1880/lreeolutlan  of  plotter  Inc/Inch 
Integer  |e/0/  (plotter  eldth  In  Inehee 
Integer  lt/0/  IcMart  tupe  0  a  1/2  froe  plot  edge 

1  I  s  I  1  f roe  plot  edge 

Integer  |r/2/  1 2==eu  1 1 1  pen  plotter  elth  l*  offeet 

Integer  lx/0/  I0=run  In  leeedtate  eode 

INTEGER  NEWEN/9S3/  1333  FDR  DL2T12,-3  FDR  CfiLCOMP 
INTEGEReZ  BUFFER  I IBUFSIC!  I ELEMENTS  3-16  OF  FIRST  RECORD  CONTAIN 
I  ITR, IMO. I0AY, IHR, IHIN.ISEC.IHSEC.JDST  OF  FIRST 
I  SAMPLE  OF  EVENT  TRACE 
INTEGER  ITR (HAXNTRC! 

INTEGER  IMOIHRXNTRCl 
INTEGER  I0AY1HAXNTRC'. 

INTEGER  IHRtMAXNTRCl 
INTEGER  IMIN(WUNTRC) 

INTEGER  ISECIMAXNTRC! 

INTEGER  IHSEC 1MAXN7RC1 

INTEGER  JOAY (HAXNTRCJ 

INTEGER  I  DEC/2/  INUMSER  OF  DECIMAL  PLACES 

BYTE  BUC2eIBUFS!Zl  I  BUFFER  EQUIVALENT  IN  BYTES 

EQUIVALENCE  IBuFFERIll  .BUIllI 

COMMON  /CHARSI2/BHT , IDEC.nFSCTOR 

NAHELIST/LISTIN/FLNAHE,XINPL0T,YINPLDT,NS»1PLES,NTRPCES, 
STRTSAHP.FSVALUE.SAMPLERA7E.  1 1, lo.  I; 
le.lt. Ir. lx.NEWPEN.EVTXIN.EVTf  IN. PROBXIN, 

PRCBY I N .  FRE5MAX , FRECM I N . EVTMAX . EVTM I N , FAC , 
SAMES DALE . NCLA5S I  NT . I DEC . NPDFSMPS , RHT , 
SCALTTRC.AFACTOR 


READ  IN  NAMELIST 

HRITEt6,el  'TYPE  IN  FILENAME  CONTAINING  PDF  NA|€LIST' 

REAOI5.304)  FLNAME 1 1 1 

OPEN  (UN!T=22,Nf*£rfLNAME  U ) .  TYPE= '  OLD ' ) 

ISTAT  =  0 

READ  122. NHL =LJSTIN. IQSTRTalSTRT) 

00  NHILEIISTST.EQ.B) 

READ  IN  TRACES 

GHIN  *  100.0 
GHAX  a  -100.0 
FMIN  a  100.0 
FMAX  a  -100,0 
I DEV  a  78 
DO  I=1,NTRACES 
OPEN  I  UNIT  =  IDEV, 

*  NWIE  =  FLNAME  II) , 

*  BLOCKS  ICE  =  2eIBUFSI2. 

*  RECL  =  IBUFSIZ/2. 

*  FORM  -  ‘UNFORMATTED'. 

*  ACCESS  a  -DIRECT', 

*  TYPE  =  'OLD') 

ISTAT  =  0 

1BU  =  1 

READ II DEV'  IBLK.IOSTAT  a  ISTAT1  BU 
IF  IISTAT.EQ.01  THEN 
ITR  1 1 )  =  BUFFER  13) 
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-  CONTINUED  - 

PROGRAM  POP 


IHOdl  a  BUFFER  UBI 
lofitm  =  BUfTsmn 

IHR(I)  »  BUFFER  1 121 
IMINUIa  SUFFER ( 131 
ISECtn  a  BUFFERIH1 
IhSECIIl  =  BUFFER  I IS) 

JDAYtl)  =  BUFFERU6) 

J  a  a 
KK  =  a 

IBLK  =  IBLK  ♦  I 

REKHIDEV'IBLK.IOSTAT  =  ISTHT!  BLK 
DO  HHILEIISTBT.EQ. B-ANO.  J.LT.HAXLEN1 
DO  Lai.IBUFSIZ 
UJ*L 

IFIK.GE.  STRTSAMPU1.AND.KK.LT.MAXLEN1  THEN 
KK  a  KK  *  1 

TRACEIKK.n  =  BUFFER  ILlaSCALFTRC  1 1) 

IF  1KK.LE.1F0FSMPS2)  THEN 
chin  =  abs:trace:kk.ih 

GHAX  a  MAXIGNAX.GMINI 
END  IF 
END  IF 
end  do 

J  a  J  *  IBUFSIZ 
IBLK  =  IBLK  ♦  1 

READ (I  DEV' IBLK. 10STAT  =  ISTAT1  BLK 
DC  00  1 1WILE I ISTRT. EC. 0. AND.KK.LT.  HAXLBI) 

END  IF 

CLOSE  1 UN IT= I DEVI 
END  DO 
C 

C—  DEFINE  LIMITS 
C 

GMIN  a  -GWX 

IF  IEVTHAX.CT.EVTMIN1  THEN 
GHIN  =  EYTMIN 
GHAX  a  EVTHAX 
END  IF 

IF  iSAHESCALEl  THEN 
EVTHAX  a  GMAX 
EVTHIN  a  GHIN 
END  IF 
C 

C—  INITIALIZE  PROBABILITY  DENSITY  CLASS  INTERVALS  TO  ZERO 
C 

00  I=1/NCLAS5!NT 
DO  Jsi.NTRACES 

ponoiseii.ji  =  a. a 

P0SIGNAL1I.J1  =  a. a 
END  00 
END  DO 
C 

c—  CALCULATE  CLASS  INTERVAL  USING  GREATEST  MINIMUM  AND 
C—  GREATEST  MAXIMUM  VELOCITY  AMPLITUDES 
C—  FOR  THE  FOUR  EVENT  TRACES  PLOTTED 
C 

CLASStNT  a  I  GHAX  -  GHIN) /FLOAT  INCLflSSINT) 

FCVALLE  a  GHIN  *  CLASS  I  NT/2. 

c—  CALCULATE  PROBABILITY  DENSITY  DISTRIBUTION  FOR  NOISE  SAMPLE  1-*FDFSMPS 
C 


CONTINUED  - 
PROGRAM  PDF 


00  Jsi.NTRACES 
DO  J=1.>«5F5MPS 
K=t 

OASS  =  .FALSE. 

00  WHILE t. NOT. CLASS.  A*fl. K.LE. NCLASSINT) 

IF  ITRACEII.J1.GT.IIUCLASS1NT.CMIN))  THEN 

t  -  K  *  1 
ELSE 

pdndiseik.J)  =  ponojseik ,j)  ♦  i.a 
CU^E  =  .TRUE. 

END  IF 

end  00  IHHILEt. NOT. CLASS) 

END  DC 

“fKSSS*  PONCISEtl  .J)  /FLOAT INPCFSMPS!  I  RELATIVE  mOUDCf 
F >«X  =  MflXlFmX.PONOISEtl.JU 
END  DO 
END  00 

C-  CRTJLATE  PRGBA8LITY  OENSITT  DISTRIBUTION  FDR  5IC*AL*NC’.SE 
C —  NP0F5MP5-1  TO  2»NPCf5MPS 
C 

NP0FSHPS2  =  2»NPDFEHPS 
00  J=1 .NTRACE5 

DO  INPCFSMPS. 1.NP0FSHPS2 
K=1 

CLASS  =  .FALSE. 

DO  WHILE !. NOT. CLASS. AND. K.LE. NCLASSINT) 

IF  (TRACEtl.Ul  .OT.  (K*CLASSINT-n3NINn  THEN 

K  =  K  *  1 
ELSE 

PDSIGNflLtK.Jl  =  POStGNALtK.J)  *  1.8 
CLASS  =  . fRUc. 

END  IF 

END  DO  i WHILE  I. NOT. CLASS) 

END  DO 

M^SSTpOSI^tU^FU«T,^S,  (RELATIVE  FRE5UENCT 
FW  =  WlFmX.PDSIGNALtI.Jl) 

END  00 
END  DO 

IF  (FREflmX.LT.FREOHIN)  THEN 
FREBHAX  s  FHAX 
FREBMIN  =  FMIN 
END  IF 

C-  SET  STATUS  TO  8  AND  THEN  TRT  TO  READ  ANOTHER  WANELIST 
C 

ISTAT  =  8 

REW3!22.NH-3_15TIN.  IQSTAT=ISTAT) 

END  DO  IWHlLEtISTAT.EB.81 
30*  FORMAT (Al6) 

END 
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LTj.in.niv-k.1  v1n’TWJ»-KrK'rj  TKrVKV>TWVVWr1,-yVS 


subroutine  pajuv-ajRoix.N.roef.iP) 


Modified  fro*  Grag-Noodsard  ARMA  •poctral  anolg»lB 
pocWag*  for  tho  IBH  PC  bg  o.hagwara  1-23-96 

Thl»  oudrouilno  soooulM  t>*  3urg  oooff lolonU  oo 
dooorlbod  bg  Kog  and  Harp  I*.  113621. 


Pqra»ot»r  (Ho*Ordor=S01 
ParaBotor  (Ma*PU=10000) 

InUgor  N 

tool  X  IN] 

R»a;  RIMaxOrdBr) 

Rto!  C  IHajturaor— -1 ) 

R»a !  S  IHaxC.-aor— 1  ] 

tool  Coof 1  IP] 

Roai  BTlHaxOrdor-1  #21 


Inuobor  of  Inoui  pi# 
I t 1 oo  oorloo  InpuU 


DIMENSION  31 (■ axpU) .32 laoxpiol .PI EoaxpU) .P2 CbobpU) 

C 

IftN.gt.MoKfto)  lh«n 
WrlUI6.ll  N 

1  Por*ot  I  lh0,  'Mjador  of  Input  point*  iooood#  boxIbub'1 

HrltotU.l)  N 

Stop  'Po(r_ftr3URG  -  orror  1‘ 

End  If 
C 

C—  CREATE  a  ZERO  HEfW 
C 

suN=«.a 

00  1=1. N 

SUK=SUh*X(!l 

Eno  Oo 

XBf«=sun/N 


C—  INITIALIZE  EQUATIONS 


1ENC=n-1 
00  Jsl.IENO 


Ptij)  =x(j»u-xafw 

qkji=xiji-xb« 

o 

Tw=a 
IEND=**-1 
DO  JS2.IEN0 

TXP=TW>  ♦  ( X (J)  -XBfW )  u2 
o 

3(1*9)  =  l XU  ]  -XBW!  b«2*(X  IN)  -XBAR1  *«2*2»TMP 

THP=0 

IEND=N-1 

00  Jsl.IENO 

Ttf>=TH,*Pl(J]«onj] 


Rtl*0) =TMP 

C! 1*9) =2* IRE  1*0) /SI  1*0) 1 
IEN0=»(-2 
00  Jsl.IENO 

P2  ( J1  =P  1 1 J*  1 )  -C  U  *0  ]  >0 1  IJ*  1 ) 
02IJ)=ai  IJl-Ctl*9l«Pl  (J1 


T=U*CU*0)bb2]bS11^)-4»C  !1*0)bR  (1*91 

S(t*ll=T-iat  (N-0-ll-CU*0)«Pl  (N-0-l))B*2-!Plll)-CU*«)»ai  II!)  «—2 
00  NISI. IP 

IEND=*-2 


r:v:v:'*»:v:v>v 


y.V.v^p^.’p'/.v.'.O.'.v.v  vVvV1 


CONTINUES  - 

SUBROUTINE  P0LRJV8URG1X,N.C0EF. IPl 


00  IK=1.IEND 

PltHCl^IIKl 
01  (lx:  =02  UK) 

End  do 

TMP=a 

IEND=N-N1-1 

DO  K=t.IENO 

THP=THP«Pl  Klofll  IK1 

Ena  do 

Rt:*NU=Tw> 

TNP=2 

IENC=N-N1-1 

oo  k=;,:enc 

TMP=TMP*(PUK!»«2*ai  IK1««2) 

Ena  Oa 
S  I1*NU  =TMP 

C(l«NI)=2»(Ra*Nl)/S(l«Nt)! 

IDC=n-*U-2 
00  J=1 , IENO 

P2lJ!rf’.  OI!-Cll*Nll«fll  1J»U 

Ene  do 

T=(l»Ct’.*NU»«2)»S[l*Nl)-*-C!lnNn  ■filial) 
T1=T- IQl  lN-Nl-U-CIl«NU»PHN-Nl-in*«2 
5a^t»n=Ti-(P!(i)-ca*NU*>:ii)>««2 


- NOW  CSLCXnTE  THE  3' 5- - —  — 

DO  KI=1.IP 

flTtru2)=t 

n 

DC  I0RDER=I < IP 

00  U=UI0R0ER 

BT(KI.l)s«TIKI.2) 

End  do 

Nl=IORCER-I 

IBOU*l 
DO  K=WIEND 

IT=N1*1 

IF(IT.EO.K)WIRST=0 
IFUT.NE.  K!  RFIRST=flT  UHC.l) 

ITl=»U*l-< 

IFim.ES.  IT)RSEXN2=2 
IFIJT1.NE.  IT)fiSECONC=flTU*Nl»HC.U 
BTll-HC.21=flFIRST-CU»NU«ASEC0N0 


00  IK=l,IP 

COEF  i  IK)=-flT  (IK*1»2) 

End  do 

RETJRN 

END 


SJSRCJTINE  RANKFREOIN.OATB.RflNK.FREQ.RAMCHAX) 

Thl«  •uta-oulln*  ooaput»«  lh«  rank  »vd  fr«QU«rcv 
of  oceu-ro^M  for  Men  doio  Mint  In  or  orraj  of  N 
doia  valu#«. 

RES.  DATA  INI  ICCNTA1NS  DATA 
REAL  RWK!N)  I  CONTAINS  RANK 
RES.  FREQIN)  ICCNTAINS  FREQUENCY 
REAL  RAMOIAX 

initially  give  all  a  data  points  a  rank  of  1.0 

and  a  FREQUENCY  OF  8. 0 

DO  1=1 -H 

RANK!’!  =  1.0  UNITIALIZE  RANK 

ARE;; I)  =0.0  UNITIALIZE  FREQUENCY 

END  DC 

CALCULATE  RANK  OF  DATA  VALUE 
AND 

CALCULATE  FREQUENCY  OF  DATA  VALUE 

DO  1*1. N 
DO  j=;.m 

IF  (DATA  ( I) .  GT.OATAl.il 1  THEN 
RANK ( I I  =  RANK ( I J  ♦  1.0 
END  IF 

IF  (DATA  1 1).  ED.  DATA  (J1 1  T>€N 
FREQ ( 1 1  =  FREQ (II  *  1.8 
END  IF 
END  DC 
END  30 

KLUST  RANK  OF  DATA  ACCORDING  TO  FREGUENCT 
DO  1=1. N 

Rf**;;)  =  RANK ( II  ♦  (FREQtD-1.01/2.8 

RAMKHAX  =  MAXI  RANK!!!  .RANKHRX) 

END  DO 

return  ranks  and  frequencies 
return 

END 


I  INCREASE  RANK  Of  DATA 
I  INCREASE  FRE3LENCT  OF  DATA 
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K'V  umv  vm*  ->  -vK7U(7'>rrj^^jr^yxV'XyvJ^j’v  '»VL"V  '."VwVvv.  ^\^V.  r?  R»  v'*  * « 


C 

C- 

c- 

c 


SUBROUTINE  RANJCTEST  IN. XDATA.YOATA.K. TEST.  ACCEPT) 

.  SUBROUTINE  TO  PERFORM  RANK  SUM  TEST 

,  moody. gratb ill  ;  boes 

RES.  XOflTSINl 
REAL  YDATA(N) 

INTEGER  N  (NUMBER  Of  SAMPLES  IN  RUN  TEST 
INTEGER  K  I ACCEPT  OR  REJECT  CRITERIA 
REAL  21200-3)  1H0RK  ARRAT  FOR  ORDERING  SAnP-E. 
LOGICAL  ACCEPT 

> _  COMPUTE  RANJC  SUM  TEST  PCR  NINCONS  Of  S . ZL  N 

N2  =  2»N 
DO  1=1. N 

2(1. n  =  AaStXDATAir.  > 

2II+N.1)  *  AB5!TCA7S::n 

fun  QQ 

CALL  RWKFRE0tN2.2(l.l)  >2(1.2)  .2(1.3) , RANKMAX ) 
TY  s  0.0 
00  1=1. N 
J  =  I  ♦  N 
TY  =  TY  ♦  2U.2) 

END  DO 


C 

C 

C 

C 


Cl  =  FLOAT!  N2  ♦  1) 

C2  =  FLOAT  (N=N1 

EXV  =  FLOAT  (N)=Cl/"2. 

Vf*  =  C2»CI/12. 

TEST  =  A0SITT  -  EXV) 

IF  (TEST. GE. FLOAT tK)  1  THEN 
ACCEPT  =  .FALSE. 

ELSE 

ACCEPT  =  .TRUE. 

ENO  IF 
RETURN 
ENO 


(EXPECTED  VALLE  OF  Ty 
(VARIANCE  OF  Ty 

I REJECT  TIE  NULL  HYPOTHEC 1 5 
(ACCEPT  THE  NULL  HYPOTHESIS 


$ 


T.Vl 
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•  j 


SUBROUTINE;  RflNK27B0(NSAMPL£S.AMPTRC,TRC.N,MEW.  INC.TEST.SEGLEN, 
♦  PATTERN,  IST.RANGE.RHRX.RMIN.NRANK.RANK1 

C 

C—  Th  I  •  eubroutlne  ooeeutee  the  rank  mM  teet  a*  deeorlbed 

C —  In  the  Introauotlan  to  the  Theory  of  Statletlae.  Mood.  Orayblll 

C—  ana  Soee,  HcGra»-N 1 1 1 ,  3rd  Edition*  19?4.  The  obeervattane, 

C —  xll),  1=1, N.  ere  taken  from  the  repreeentailve  baav^roLna  nolee 

C—  oeeiMO  to  bo  prooont  In  the  flrot  2.5  oooondo  of  the  event 
C—  trace.  The  eeoond  group  of  obeervat lone*  gill*  1=1, N,  1e 
C“  ehoeen  free  a  elndoe  of  N  eaeelee  eoved  Inoreeentol  )y  along 
C —  the  event  trooe  for  eaoh  exaction  of  the  rank  eue  teet. 

C —  The  ranx  eue  teet  etatletle  le  ooeputod  for  eooh  etndoe 

C —  of  obeervatlane,  gill.  1=1. N.  and  the  reoreeentatlve  nolee 

C—  obeervatlone.  xltl.  1=1. N.  Detection  oocure  when  the  rank 
C —  eue  of  the  yfll  ranke  exceeae  the  threenoid  for  prooaole 

C —  nolee.  The  autoeatlc  flret  «"rlval  plok  te  the  flret 

C—  xero  oroeelng  preoealrg  the  flret  Modified  elope  value  thot 
C —  exceeas  the  eenleue  eoalfled  elope  value  of  the  repreeentatlve 

C —  nolee  eaeipiee.  mill.  1=1. N,  occurring  in  the  flret  deteotlon 

C—  w  I  naow . 

C 

REAL  TEST 
REflL  RMAX 
REAL  RMIN 

REflL  PATTERN!'.)  I  NULL  KTPCTHESIS  =  l.B  //  ALTERNATIVE  =  2.8 

RE~_  TgCHl  IARRAY  CONTAINING  FEATURE  OF  SIGNAL  TO  RANK 

REFL  W?*Rc:U  I  ARRAY  CONTAINING  WPLITUDES  OF  SIGNAL  POINTS 
REPL  SEGLENU1  I  LENGTH  Cf  PATTERN  SEGMENT 

REA.  WAh  I  EXPECTED  VALUE.  OF  RANK  SUM  TEST 

REAL  RANGE  I  RANK  SUM  RANGE  OF  VALUES 

INTEGER  NSAMPLES  INUMfiER  OF  OBSERVATIONS  IN  EVENT  TRACE 

INTEGER  NRANK  INUMSS3  OF  RAMC  SUM  VALUES  COIf’L'TED  FOR  TRACE 
INTEGER  INC  I  INCREMENT  WINDOW  MOVES  DOWN  TRACE 

INTEGER  1ST  IFIR5T  WRIVA  PICK  INDEX 

INTEGER  N  IHLRBER  OF  SAMPLES  REFERENCE  BACKGROUND  NOISE  WINDOW 
REAL  TI3M.31  l WORK  ARRAY  TO  RANK  NOISE  l  SIGNAL  H1MJ0W5 

REAL  RANK  I  *  I  I  SUMPTION  TY  OF  RAHCS  OF  SIGNAL 

LOGICAL  FOUND/.  FALSE./ 

LOGICAL  SECOND/. FAuSE./ 

LOGICAL  FIRST/. TRUE./ 

C 

C—  COMPUTE  RANK  SUM  TEST  OVER  TRACE 

c 

N2  =  2eN 
RMAX  =  MEW 
RMIN  =  MEW 

EHMAX  =0.8  I EVENT  NOISE  MAXIMUM  MODIFIED  SLOPE  MAGNITUDE 

C 

C—  STORE  REFERENCE  BACKGROUND  NOISE  IN  WORK  VECTOR 
C—  FIND  EVENT  NOISE  MAXIMUM 
C 

DO  1=1. N 
L  =  1-1 

PATTERN  I NSFM’LE5-LI  =1.0 
SEGLENINSAfPLES-L)  =  l.B 
PATTERN  1 1 !  =  1.8 
SEGLENC)  =  1,8 
TII.l)  =  ABS (TRC  t Ie2~T  1 ) 

TIIvN.ll  =  ABS ITRCI i=2! ) 

ENMAX  =  MRXtTII.il .ENMAX) 

ENMAX  =  MAXtTIIvN.ll  .ENMAX1 
END  30 

ENMAX  =  ENMBXet.81  I  THRESHOLD  FOR  FIRST  ARRIVAL  TICKER 

C 

C—  COMPUTE  RANK  SUM  FOR  A  WINDOW  OF  LENGWT  N  COMPARED 
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TOT 


-  contimjed  - 

SUBROUTINE  RAMI2738  (NSA^LES.  AHPTRC.TRC.N. MEAN.  INC. TEST.SEGLQL 


WITH  THE  REFERENCE  EBCXGROUfC  NOISE  WINDOW  OF  LENCHT  N 
IR  =  8 

SLIDE  THE  WINDOW  DO»J  THE  TRACE  IN  INCREMENTS  =  INC 


DO  I=l.(NS«m£S-NMNC 
IR  =  IR  -  1 
DO  J=1.N 
JJ  =  J  *  (I-U 
T(J-N.l)  =  ABSITRCIJJI) 

END  DO  IJsl.N 

CALL  RAhKFREQ(N2.T(l,l! . TI1.21  .Tll.3)  .RflNKMAX) 

rr  =  a. a 

DO  J=:.N 

rr  =  tt  ♦  tij*n.2i 

END  DO 


RMfiX  =  HflXirr.RMAX) 
RMIN  =  MINlTr.RMlNI 

rankiiri  =  rr 

END  DO  II 
NRAnk  =  IR 
RANGE  =  RMflX  -  RMIN 


l  RANK  SUM  MAXIMUM  VPUJE 
I  RANK  SUM  MINIMUM  value 
I  RANK  SUM  OF  TVE  IRlh  WINDOW 


ISIGNS.  tt(D  NOISE  RANK  SUM  INTERVAL 
NOW  DETERMINE  THeESHOLD  =  2788  *  RMIN.CR.8.75«(RMAX-RMIn!  ♦  RMIN 


IF  URMAX-RMIN1.LT. 3688.)  THEN  IU)W  SIO«L  TO  NOISE  RANK  SUMS 
TEST  =  IRMAX-RM!N1»8.7S  »  RMIN 
ELSE 

TEST  =  2700.  *  RMIN 
EMI  IF 


C 

C—  FIND  SIGNAL 
C 


ILDOKING  FOR  FIRST  DETECTION  WINDOW 


FIRST  =  .TRUE. 

00  I=1.NRANK 
K  =  II-llalNC  *  1 

IF  (RANKU1.GT.TEST1  THEN  (PROBABLE  SIGNAL  PRESENT  IN  WINDOW 
00  L=K.IK*(INC-1!> 

PATTERN (LI  =2.0  IREJECT  HB 

SEOLENIL)  =  1.8 
END  DC 
L  =  K 

00  WHILE  (FIRST. AND. IL.LT. (K+1N-1! 1 ) ) 

1ST  =  L 

IF  (ASS  (7RC IL I ) .  ITT.  ENNAX)  77CN 
MM  =  L  -  1 

DO  WHILE  (FIRST. AND. (MM. GE. 1)1 
1ST  =  MM 

IF  (AMPTRCIMM*l!»(Vf,TRC(MM).LT.0.81  THEN 
FIRST  =  .FALSE.  IFOUND  IST=2ER0  GROSSING 

ELSE 

ICOWTINUE  LOOKING  FOR  ZERO  CROSSING 
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MM  =  MM  - 
END  IF 

END  DO  1 WHILE 
ELSE 

L  =  L  *  1 
END  IF 

END  00  I WIILE 
IF  (FIRST!  THEN 
MM  =  L  -  1 


(OBSERVATIONS  IN  WINDOW  NOT  .GT.  NOISE 
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CONTIMED  - 

SUBROUTINE  RSMC27B01NSRHPLES.SMPTRC.TRC.N.HEPN. INC. TEST. SEGLQl. 


00  WHILE  (FIRST.SKJ.[HH.BE.l)l  IP1CK  1ST  ZERO  CROSSING 
1ST  =  HH  I  BEFORE  EMI  OF  WINDOW 

IF  (BHPTRCtHH*l)rtMPTRCIMHl.LT.B.B)  THEN 
FIRST  =  .FSLSE.  (FOUND  IST=ZERC  CROSSING 

ELSE 

I  CONTINUE  LOOKING  FOR  ZERO  CROSSING 


HM  =  MH  -  1 
END  IF 

END  00  I  WHILE 
END  IF 
ELSE 

00  Lac.  (K*(  INC-111 
PATTERN (L)  =  1.0 
SE3LEN1L)  =  t.B 
END  DO 
END  IF 
END  00  II 
RETURN 
END 


ISCCEFT  HB 


m 


M 


SUBROUTINE  RUNEXPtN.XDATB.YmTB.IZa.NRUNS.EXV.VBR) 


Thl*  subroutine  oaloulat**  U>*  run  Iwl  statistic 
for  tso  arrays  of  data#  sill.  1=1. N  and  yll).  1=1. N. 

Ths  thr*sncld.  l2.  tor  ths  nt-mc  m~  of  run*  bslow  rn  I  oh  «» 
reject  ths  null  hypothesis.  Ft*)  =  Fly),  I •  provided  oo 
Input  to  ths  subroutine.  It  is  oaloulatod  outoldo  ths 
•ubroutln*  using  th*  norsal  distribution  and  o  spsclflsd 
alpha  level.  Th*  algorithm  I*  bo**d  on  th*  daacrlptlon  of 
th*  rvn  t**t  glvsn  In  th*  Introduction  to  th*  Theory  of 
Statlatfcs.  Hood.  Gragblll  and  Bo*s.  MoGca*-H 1 1 1 ,  3rd  Edition. 
197t. 


REAL  XDATAIM 
REAL  TORTS  IN) 
RESL  EXV 
RES.  VSR 
INTEGER  N 
INTEGER  1 22 
RESL  21222.2) 


I  ©PETTED  VALUE  OF  NUMBER  OF  RUNS 

i vsr i shoe  of  number  of  runs 

I  NUMBER  OF  SAMPLES  IN  RUN  TEST 
ISCCEPT  OR  REJECT  CRITERIA 
I  WORK  ARRAY  FOR  ORDERING  SfWLSS 


C 

C—  COMPUTE  RUN  TEST  OVER  TRACE 
C 


N2  =  2sN 

□0  1=1. N 

211. 1)  =  RBS IXDATSI I) ) 

211.21  =  1.2 

2II-N.il  =  ABS (YOSTAt I) ) 
2II-N.21  =  -1.2 
END  00 


C 

C—  ORDER  COMBINED.  X(I)  AND  Til).  OATH  INTO  ASCENDING  ORDER 
C 


00  1=1, NJ- l 
00  J=I»1.N2 

IF  tZIJ.D.LT.ZlI.in  THEN 
temp  =  za.ii 

ZII. l)  =  ZIJ.l) 

ZIJ. l)  =  TEMP 
TEMP  =  ZII.2) 

211,21  =  2IJ.2) 

ZIJ.2)  =  TEMP 

EJO  IF 
END  00 
EMI  00 


C 

C—  COUNT  Th£  WJfflER  OF  RUNS 
C 


ZSIGN  =  211.2) 

WUNS  =  1 
00  I=2,N2 

IF  t  (ZSIGNsZU  ,2) )  .LT.2.2)  THEN 
ZSIGN  =  ZII.2) 

NRUNS  =  NRUNS  ♦  t 
END  IF 
E>C  00  II 


COMPUTE  EXPECTED  VALUE  AND  VARIWICE  OF  THE  NUMBER  OF  RUNS 
ASSUMING  THE  HULL  NTPC  THESIS  25  TRUE 


Cl  =  FLOAT I2sNsNl 
C2  =  FLOAT (N2) 

EXV  =  C1/C2  *  1. 

VAR  =  ICl*(CI-C2n/(C2**2*IC2-l.n 
RETURN 
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I 


■r, 


i 

*  w 

'4 


jfl 


CH 


•fit 


% 


y 


v 

& 


*4 


<~jsrcut::»’C  SLyEosAMPLEJ.TRc.si 

c 

C —  Tht •  •ufarouiln*  ooaputoo  th«  Intogrotod  slop*  of 

C—  a  Mlnlo  •  i  iyia  for  dloerolo  oogoonLo  of  tho  •  I  goal 
C—  *  I  th  oonilnueuo  oiopo  direction. 

C —  Kathloon  A.  Aldon,  March  16-  1986. 

C 

REAL  SIU  I  MODIFIED  SLOPE  OF  OBSERVATION 

REAL  TRCIl)  I  TRACE  CONTAINING  SEISMIC  SIGNAL  CBSERVATIONS 

INTEOER  NSAMPLES  INUMBER  Of  OBSERVATIONS  IN  TRACE 

C 

C—  ASSUME  SAMPLE  CAN  OF  TRACE  ItPUT  IS  ZERO 

C 

C 

C—  COMPUTE  SEISMIC  TRACE  MODIFIES  SLOPE 
C 

3(11  =  TRC (21  -  TRC1I)  ISLDPE  OF  1ST  SAMPLE 

00  I  =2 -NSAMPLES  -  1 

Sill  *  ( (TRCi 1*11 -TRC (II 1 ♦ t TRC  til —TRC 1 1  —  1 1 1 1/2.  Illh  SLOPE 

C 

IF  (S(II»S(I-11.GT.0.0)  THEN  ISLOPE  DIRECTION  IS  CONTIGUOUS 
Sill  =  SIX)  ♦  Sll-U  IINTEGRATE  SLOPE 

END  IF 
END  00 

S I  NSAMPLES!  =  TRC  (NSAMPLES!  -  TRC  (NSAMPLES- 1 1  ISLOPE  LAST  SAMPLE 
RETURN 


|Wfw  irairvotwffWTOliruwvirt  vf  ',>  >v  vrsJJWW  TOfTOTOWTOTOTOW^ 
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c— 

c— 

c— 

c— 

c— 

c— 

c— 

c— 

c— 

c 


SUBROUTINE  SIGNE*PCM.XDBTft.TfDBTB.Il.SN^v.TCri 

Thl-  wibrouil^  ealeulfli-  **•  M 

for  ih.  obnlul.  volo.  of  in  <rrc*.  of  doio,  .III. 

*111.  l=l.N.  Th.  ihrnnold.  k.  nieh  —I  b.  In.  U>on  ^ 

U«i  .latl.ile  io  r.jni  ih.  null  hgpoihnl..  Fl»  * 

provide  o.  Input  l*  U-  nbroultn..  11  *•  eol~‘ol!l~^ 
ih.  .ubnouiln.  using  ih.  nor.ol  distribution  and  a  .pnlfln 
oloho  tn.l.  Th.  oi  gor  I  ih.  I.  bond  on  ih.  dnerlpUon  of 
ih.  .Ign  ini  glv«n  In  th.  lhtroAjcilon  to  th. 

Statistic  Hood.  Groyblll  and  Bo...  MoGrw-M  1 1 1 .  3rd  Edition, 

197*. 

REAL  XOfl'SIN) 

INTCCD^K^^1  I  THRESHOLD  FOR  THE  TEST  STATISTIC 

l  I  NUMBER  Of  SAMPLES  IN  SI 04  TEST 

REAL.  SN  I  SUMPTION  OF  THE  POSITIVE  SI045  OF  (TUI  -  XlIH 

m  1 EXPECTED  VALUE  OF  SUMMATION  OF  SIGNS 

^c7  I  ABSOLUTE  VALUE  OF  ISN  -  EXV) 


C 

c— 

c 


COMPUTE  SIGN  TEST  OVER  TRACE 

SN  a  0.0 

00  1=1. N 

IF  IASS IXDSTAtl) l.LT.ABSITOATAI I) 1 )  THEN 
SN  =  SN  »  1.0 
END  IF 
END  DO 

COMPUTE  EXPECTED  VALUE  OF  SUMMATION 
AND  THE  TEST  STATISTIC 

EXV  =  FUDATINl/2. 

TEST  =  ABS (SN  -  EXV) 

RETURN 

END 
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TABLE  C.l  (Continued) 


*- 3 


X 

3 

a 

K» 

m  ■  13 

m  ■  14 

m  ■  IS 

m  ■  16 

m  *  17 

m  •  18 

m  «  19 

m  «  20 

24 

.527 

25 

.473 

26 

.420 

.500 

27 

.367 

.450 

.524 

28 

.316 

.400 

.476 

29 

.268 

.352 

.429 

.500 

30 

.224 

JOS 

.384 

.456 

.521 

31 

.182 

.261 

.338 

.412 

.479 

32 

.147 

.220 

.296 

.369 

.438 

.500 

33 

.116 

.182 

.254 

.327 

.396 

.461 

.519 

34 

.090 

.148 

.216 

.287 

.356 

.421 

.481 

35 

.068 

.120 

.181 

.249 

.317 

.382 

.444 

.500 

36 

.051 

.095 

.150 

.213 

.280 

.345 

.407 

.464 

.517 

37 

.035 

.073 

.122 

.180 

.244 

.308 

.370 

.4  29 

.483 

38 

.024 

.055 

.099 

.151 

.211 

.273 

.335 

.394 

.449 

39 

.015 

.041 

.078 

.125 

•  ISO 

.239 

.300 

.359 

.415 

40 

.009 

.029 

.060 

.102 

.152 

.208 

.267 

.325 

.382 

41 

.004 

.020 

.046 

.082 

.127 

.179 

.235 

.293 

.349 

42 

.002 

.012 

.034 

.065 

.105 

.153 

.206 

.262 

.317 

43 

.007 

.024 

.050 

.086 

.129 

.178 

.232 

.286 

44 

.004 

.016 

.038 

.069 

.108 

.153 

.204 

.257 

45 

.002 

.010 

.028 

.055 

.089 

.131 

.178 

■  .229 

46 

.006 

.020 

.042 

.073 

.111 

.154 

.202 

47 

.003 

.013 

.032 

.059 

.092 

.132 

.177 

48 

.001 

.009 

.024 

.046 

.077 

.113 

.155 

49 

.005 

.017 

.036 

.062 

.095 

.134 

SO 

.002 

.01 1 

.027 

.050 

.080 

.115 

51 

.001 

.007 

.020 

.040 

.066 

.098 

52 

.004 

.014 

.031 

.054 

.083 

53 

.002 

.010 

.023 

.04  4 

.069 

54 

.001 

.006 

.017 

.034 

.058 

55 

.004 

.012 

.027 

.04  7 

56 

.002 

.008 

.020 

.038 

57 

.001 

.005 

.015 

.  .030 

58 

.003 

.010 

.023 

59 

.002 

.007 

.018 

60 

.001 

.005 

.013 

61 

.003 

.009 

62 

.001 

.006 

63 

.001 

.004 

64 

.002 

65 

.001 

66 

.001 
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TABLE  C.l  (Continued) 


m  •  12 

m  ■  13 

m  »  14 

m  *  15 

m  •  16 

m  *  |  7 

m  *  18 

m  *  19 

m  ■  2  n 

34 

.524 

35 

.476 

36 

.431 

.522 

37 

.385 

.478 

38 

.34  2 

.435 

.521 

39 

.299 

.392 

.479 

40 

.260 

.352 

.439 

.519 

41 

.223 

.312 

.399 

.481 

43 

.190 

.274 

.360 

.44  3  ■ 

.518 

43 

.158 

.239 

.323 

.405 

.482 

44 

.131 

.206 

.287 

.368 

.446 

.517 

45 

.106 

.175 

.253 

.332 

.410 

.483 

46 

.085 

.148 

.221 

.298 

.375 

.449 

.516 

47 

.066 

.123 

.191 

.265 

.34  1 

.415 

.484 

48 

.052 

.101 

.164 

.235 

.308 

.381 

.451 

.516 

49 

.039 

.082 

.139 

.205 

.277 

.349 

.4  19 

.484 

50 

.029 

.065 

.116 

.179 

.24  7 

.318 

.387 

.453 

515 

51 

.021 

.051 

.096 

.154 

.219 

.287 

.356 

.4  2  2 

4X5 

52 

.015 

.039 

.079 

.131 

.192 

.258 

.326 

.392 

4  55 

53 

.010 

.030 

.063 

.110 

.168 

.231 

.297 

.363 

42h 

54 

.007 

.022 

.051 

.092 

.145 

.205 

.269 

.334 

.347 

55 

.004 

.016 

.04  0 

.076 

.124 

.181 

.242 

.306 

,3<.  K 

56 

.002 

.011 

.031 

.062 

.106 

.158 

.217 

.279 

.34  t 

57 

.001 

.008 

.023 

.050 

.089 

.138 

.193 

.253 

.314 

58 

.001 

.005 

.017 

.04  0 

.074 

.119 

.171 

.228 

288 

59 

.003 

.012 

.031 

.061 

.101 

.150 

.205 

.262 

60 

.002 

.009 

.024 

.050 

.036 

.131 

.183 

.239 

61 

.001 

.006 

.018 

.041) 

■072 

.113 

.162 

2  lh 

62 

.000 

.004 

.014 

.03  2 

.060 

.098 

.143 

1  94 

63 

.002 

.010 

.025 

.049 

.083 

.125 

.174 

64 

.001 

.007 

.019 

.040 

.070 

.109 

.155 

65 

.001 

,nu5 

.015 

.032 

.059 

.094 

137 

66 

.000 

.003 

.01 1 

.026 

.049 

.081 

.120 

67 

.002 

.008 

.020 

.040 

.069 

,  .10? 

68 

.001 

.006 

.016 

.033 

.05  8* 

.04  1 

69 

.001 

.004 

.012 

.027 

.04  9 

.074 

70 

.000 

.002 

.009 

.021 

.04  1 

067 

71 

.001 

.006 

.017 

.033 

.1)57 

72 

.001 

.005 

.013 

.027 

.04  X 

73 

.00(1 

.003 

.010 

.022 

.04  1 

74 

.000 

.002 

.007 

.018 

.034 

75 

V  .001 

.005 

.014 

.028 

76 

.001 

.004 

.01  l 

.0  23 

77 

.000 

.002 

.008 

.018 

’.  v.v.v. , 
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TABLE  C.l  (Continued) 


It  *  4 


*  m  M2 

m«  13  m  M4  m  *  15  mM6  m«n 

m  M8 

m  -  19 

3 

• 

u 

O 

78 

.000 

.002 

.006 

.015 

19 

.001 

.004 

.01 1 

to 

.001 

.0D3 

.009 

81 

.000 

.002 

.007 

82 

.000 

.001 

.005 

13 

.001 

.004 

.000 

.003 

85 

.000 

.002 

86 

.000 

.001 

87 

■ooi 

88 

.000 

89 

.OU0 

90 

.000 

n  -  5 


m  ■  5 

m  •  6 

m  *  7 

m  •  8 

m  »  9 

m  M0 

28 

.500 

29 

.421 

30 

.34  5 

.535 

31 

.274 

.465 

32 

.210 

.396 

33 

■  1SS 

.331 

.500 

34 

.111 

.268 

.438 

35 

.075 

.214 

.378 

.528 

36 

.04  8 

.165 

.319 

.472 

37 

.028 

.123 

.265 

.416 

38 

.016 

.089 

.216 

.362 

.500 

39 

.008 

.063 

.172 

.311 

.449 

40 

.004 

.04  1 

.134 

.262 

.399 

.523 

41 

.026 

.101 

.218 

.350 

.477 

42 

.015 

.074 

.117 

.303 

.430 

43 

.009 

.053 

.142 

.259 

.384 

44 

.004 

.037 

.111 

.219 

.339 

45 

.002 

.024 

.085 

.182 

.297 

46 

.015 

.064 

.149 

.257 

47 

.009 

.047 

.120 

.220 

48 

.005 

.03  3 

.095 

.185 

49 

.003 

.023 

.07  3 

.155 

50 

.001 

.015 

.056 

J27 

51 

.009 

.(W! 

.103 

52 

.005 

.030 

J82 

138 


TABLE  C.l  (Continued) 


«  ■  5 


»  m  »5 

m  •  6  m  *  7  m  ■  8 

m  *  9 

m  ■  10 

S3 

.003 

.021 

.065 

54 

.002 

.014 

.050 

55 

.001 

.009 

.038 

56 

.006 

.028 

57 

.003 

.020 

58 

.002 

.014 

59 

.001 

.010 

60 

.000 

.006 

61 

.004 

6: 

.002 

63 

.001 

64 

.001 

65 

.000 

n  •  6 


X 

m  *  6 

m  «  7 

m  *  8 

m  »  9 

m  •  10 

39 

.531 

40 

.469 

41 

.409 

42 

.350 

.527 

43 

.294 

.473 

44 

.242 

.418 

45 

.197 

.365 

.525 

46 

.155 

.314 

.475 

47 

.120 

.267 

.426 

48 

.090 

.223 

.377 

1523 

49 

.066 

.183 

.331 

.477 

50 

.047 

.147 

.286 

.432 

51 

.032 

.117 

.245 

.388 

.521 

52 

.021 

.090 

.207 

.344 

.4  79 

53 

.013 

.069 

.172 

.303 

.437 

54 

.008 

.051 

.141 

.264 

.396 

55 

.004 

.037 

.114 

.228 

.356 

56 

.002 

.026 

.091 

.194 

.318 

57 

.001 

.017 

.071 

.164 

.281 

58 

.011 

.054 

.136 

.246 

59 

.007 

.04  1 

.112 

.214 

60 

.004 

.030 

.091 

.184 

61 

.002 

.021 

.072 

.157 

62 

.001 

.015 

.057 

.132 

63 

.001 

.010 

.044 

.110 

139 
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TABLE  C.l  (Continued) 


m  •  9 

n  •  9 

n 

•  10 

Jt 

m  *  9 

m  «  10 

X 

in  *  9  m  *  10 

X 

m  »  1 

86 

J00 

122 

.000  .004 

121 

.12 

87 

.466 

123 

.000  £03 

122 

.10 

88 

.432 

124 

.000  .002 

123 

.09 

89 

.398 

125 

.000  .001 

124 

-08 

90 

.365 

.516 

126 

.000  .001 

125 

.07 

91 

.333 

.484 

127 

.001 

126 

.06 

92 

.302 

.452 

128 

.000 

127 

.05 

93 

.273 

.421 

129 

.000 

128 

.04 

94 

.245 

.390 

130 

.OuO 

129 

.03 

95 

.218 

.360 

131 

.000 

130 

.03 

96 

.193 

.330 

132 

.000 

131 

.02 

97 

.170 

.302 

133 

.000 

132 

.02 

98 

.149 

.274 

134 

.000 

133 

.01 

99 

.129 

IJ 

Aw 

OO 

135 

.000 

134 

.01 

100 

.111 

.223 

135 

.01 

101 

.095 

.200 

136 

.00 

102 

.081 

.178 

ft  *  10 

137 

.00 

103 

.068 

.158 

138 

.00 

104 

.057 

.139 

x  m  *  10 

139 

.00 

105 

.04  7 

.121 

140 

.00 

106 

.039 

.106 

105  .515 

141 

.00 

107 

.031 

.091 

106  .485 

142 

•  .00 

108 

-.025 

.078 

107  .456 

143 

.00 

109 

.020 

.067 

108  .427 

144 

.00 

110 

.016 

.056 

109  .398 

145 

.00 

111 

.012 

.04  7 

110  .370 

146 

.00 

112 

.009 

.039 

111  .342 

147 

.00 

113 

.007 

.033 

112  .315 

148 

.00 

114 

.005 

.027 

113  .289 

149 

.00 

115 

.004 

.022 

114  .264 

150 

.00 

116 

.003 

.017 

115  .241 

151 

.00 

117 

.002 

.014 

116  .218 

152 

.00 

118 

.001 

.011 

117  .197 

153 

.00 

119 

.001 

.009 

118  .176 

154 

'.00 

120 

.001 

.007 

119  .157 

155 

.00 

121 

.000 

.005 

o 

o 

r* 

Adapted  from  Table  B  of  A  Nonperamcrric  Introduction  to  Statistics,  by  C.  H.  Kraft  and 
C.  ran  Eedcn,  Macmillan,  New  York,  1968,  with  the  permisson  of  the  authors  and  the 
publirher.  Copyright  ©  1968,  by  the  Macmillan  Company. 


APPENDIX  D 


Illustration  of  the  Results  of  the  Automatic  Detector 
ana  Picker  When  Run  on  the  Event  Traces  Recorded 
by  the  Four  Seismic  Stations ,  Lailtas ,  Marathon , 
Shaf ter  and  Tres  Cuevas  for  38  Seismic  events 


""""  . . . . 


TABLE  D.l  Event  Pile  Names,  Errors  {Automatic  Minus 
Analyst  Picks  in  Seconds),  and  the  Range  of  Rank  Sums 
for  the  3fi  Events  Used  to  Test  the  RANK27Q0  Detector 


^2!_E  NflME 

ERRC8  IN  flUTOKfVT IZ  PICK 

RBNGE  Of  RflW 

dent  i 

L 42891 419. 337 

-0.1 

4922.5 

M4289W18.053 

0.2 

5475. 5 

342891 418. 817 

8.0 

5057.0 

742891418.934 

-0.1 

5343. 5 

DENT  2 

L42891719.8E9 

1.2 

4168.5 

>•42891719.957 

1.7 

3826.2 

342891719.925 

7.4 

4805.5 

T42891719.94? 

1.9 

4242.5 

EVENT  3 

uzasiais.aia 

6.5 

4178.8 

N428S 1815. 938 

a.s 

3481.5 

342891815.018 

10.4 

3472.5 

742891814.048 

1.2 

4812.0 

EVENT  4 

L 42891 90S. 939 

-17.5 

3068.5 

X428919B6. 048 

2.7 

334S.5 

342891993.918 

-IS.l 

3361.2 

T42831S05.0B8 

0.3 

3880.0 

DENT  5 

L42892304.233 

-19.3 

5576.5 

*  4283239 4. 928 

-5.1 

4133.5 

342892394.012 

-6.3 

5496.3 

742892394.015 

0.2 

5740.2 

EVENT  6 

L42S92313.003 

0.0 

4822.5 

H42892313.022 

-IS. 3 

4760.5 

342892313.016 

0.1 

5454.2 

742892313.008 

8.6 

5583. 5 

DENT  7 

U2892327.016 

-1.8 

4615.5 

842892328.  007 

1.1 

3311.5 

342892327.018 

-0.8 

4096.2 

T 42892327. 012 

0.3 

5168.2 

DENT  8 

L 429080 19. 851 

-0.2 

4862.2 

M 42308013. 053 

0.1 

4579.5 

342300021.023 

0.2 

4569. 5 

742900019,051 

-0.1 

5564.5 

DENT  9 

L 42300136. 853 

-25.0 

4227. S 

H42300137, 019 

12.6 

4228.2 

542900138.031 

0.1 

3581.0 

742300138.021 

8.3 

4554.5 

DENT  10 

L42301821.016 

7.S 

4405.5 

H 4230 1820. 047 

0.3 

4770.0 

342901821.025 

-0.4 

4187.0 

742301828.037 

-9.1 

3232.2 

DENT  11 

L  4239285  4. 044 

-0.2 

5524.8 

M42S02055.038 

0.1 

4294.0 

342302054.051 

-26.8 

6894.0 

742902054.044 

0.2 

5984.2 

DENT  12 

L 42902289. 027 

3.0 

2862.2 

H42302210. 022 

0.7 

4252.8 

342902211.000 

5.1 

2862.0 

742302209.038 

0.3 

3713.8 

r» 


»>.v 


< 
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FILE  N WE 


TABLE  D.l  t Continued) 

ERROR  in  OJTDHBTIC  PICK  RWOE  OF  W*«  SUHS 


EVENT  13 

L42902247.013 

-0.6 

H42302246.823 

4.  4 

342302247.016 

5.  4 

T42302246.038 

0.3 

EVENT  14 

1.42910325,932 

1.4 

M42910325.026 

0.7 

342318325.016 

-22.9 

T423182 25.023 

0. 6 

■  EVENT  15 

L42313S22.048 

-1.5 

M4291SS22.048 

0.2 

342313922.335 

0. 1 

T42910S22.040 

0.1 

EVENT  16 

-33.3 

L4231 1682. 007 

H429  U  602.038 
S4291 1603. 002 
T42311602.083 
EVENT  17 
L429  U807.029 
H4291 1807. 023 
342911887.026 
T423U  807.032 
■  EVENT  18 
L42911819.039 
H4231 1819. 035 
542311813.030 
T4231 1819.033 
EVENT  13 
14231 1851.003 
M42311851.051 
542311851.011 
742311851.007 
EVENT  20 
L42311S53.017 
M4231 1353. 008 
542311358.058 
T423  U35S.056 
EVENT  21 
L42S12140. 015 
M42312140.040 
542312133.030 
T42312140.815 
EVENT  22 
L42312213.057 
M4231221 4. 052 
542312215.001 
T 423 1221 3. 057 
EYENT  23 
L. 42320528. 008 
M 42320528. 004 
542320528.081 
T 423 20528. 008 
EVENT  24 
1.42521423. 020 
*42321430.01* 
542321430.030 
T42321423.018 
EVENT  25 
L42321533. 036 
M 42321533. 023 

542321533.026 

T42321S33.03S 


0.3 

0.2 

0.8 

0.4 

0.4 

0.3 

-0.3 

-5.5 

-23.9 

0.1 

0.2 

-0.6 

0.0 

-2.4 

0.3 

-0.7 

13.4 
-6.3 

0.3 

0.4 

0.1 

10.4 
-2.3 

-0.4 

2.1 

-23.7 

0.1 

-0.1 

-3.4 

0. 1 

-1.2 

0.5 

0.3 

-0.2 

-0.3 

0.1 

0.1 

0.1 

-0.4 


145 


4145.0 
435  4.  5 
2907. 0 
3742.0 

5358.0 

5540.0 

5376.0 

5643.0 

5437.0 
4327. B 

5260.5 

5462.  a 

4522. 0 
3788.0 

4588.5 

5137. 5 

4572.5 

4144.5 

4363.0 

4579.5 

3833.5 

3035.5 

3335.5 

4139.5 

4626.5 
31742.0 
5646.0 

5059.5 

4330.0 

4199.0 

6332.0 

4632.5 

5730.0 

5143.0 

3748.0 

5518.5 

4443.5 
4067.0 
5845.0 

4014.5 

5145.5 

5799.5 
4184.0 
5536.0 

2235.3 
4353. 0 

3994.3 
3344.0 

5469.0 

5871.0 

5261.5 

5751. 5 


Vi 

ft 
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TABLE  D.l  (Continued) 


nu  wwe  ejwo*  1N  amuBTic  pick  »*<*  of  rpnk  a*s 


Denotes 


EVENT  26 
L  4232 '.823. 333 
H42321S33. 308 
542321623.323 
T42S2 1323. 333 
EVENT  27 
L42321653.33! 
*142321853.054 
5*2321858.356 
T 42321853. 028 
•  EVENT  23 
L42322-.57.314 
H42S22157.314 
S42322157.314 
T42922157.31" 

■  EVENT  23 
L423223B4.342 
*1*232238*.  325 
542322384.333 
T42S22384.335 

■  EVENT  30 
L423381 48. 326 
M423381 *8. 229 
5423381*8.323 
T42338148.823 

a  EVENT  31 
L42331453.014 
M42331443.047 

542331450.018 
T42331 *53. 81  * 

EVENT  32 
L42332123.00S 
*142332123.031 
542332123.018 
T42S32123. Z86 
EVENT  33 
L42332203. 313 
*1*2332203.035 
5*2332383.835 
T *2332303. 013 
EVENT  34 
L42351308. 057 
H4235 1388. 052 
542351303.002 
T 4235 1388. 057 
EVENT  35 
L42362237 . 016 
*142362238. 001 
542362038.014 
T42362237.016 
EVENT  36 
L42370615. 046 
*142378815.043 
S 423 788 l 5. 833 
T 423788 IS. 0*6 
EVENT  37 
L42ST17S3.018 
F479717F1.330 
S42371753.018 
T42371752.052 
EVENT  38 
L42371352. 007 
II 4237 135 1.000 
342371350.043 
T423719S1.000 

as  the  Events 


U9ed 


1.2 

-1.8 

18.1 

3.3 

as  Part 


4583.8 
2*03.5 
4355. S 
3*13.5 

of  the  Training  Set 


'Vm.-IMUTKJIV. rv  V V.V\WWiSVS\ 


«K  SMS  OF  IHTEOWTED  SLOPE  NMH  K  *  11073.  HUE  =  10050.1  KME  =  3826.* 


Discars  ikttowtiqn  or  surf  seocwts  with  tw  s«  slope  direction 


START  TIHEi  1984  2SS  17i  19<Z4. K42891719.0571  REOICWL  SCALES  STt  1.0 


m 


m/V 


6.67  13.33  20.80  26.67  33.33  40.2 

SECONDS 


53. 33  60.00 


RAMC  SLKS  OF  INTEGRATED  SLOPE  N  =  1«  K  e  11186.  NJE  s  18050.0  RANGE  *  4166.5 


QIStFTTE  INTEGRATION  OF  SLOPE  SEWEXTS  WITH  TrtE  SFWE  SLOPE  DIRECTION 


fr|r(T(WHl 


START  TtICi  1964  288  17tl9il6.  1  142891719.850  REGI04RL  SOLED  BTi  1.0 


0. 60  6.57  13.33  20. 


SECONDS' 


T"  40.00  46.67  53.33  60.00 


Event  2  Regional  (P-Lg)  =  20  Seconds 


:W.vj 


■v  ^  ,-v’.  .vt  .v.  v  >crx.w 


OJSCXTTE  IMTEOWTICH  OF  3L«  3BGWXT3  WITH  TJC  3*C  3L«  DIJECTJCW 


STWT  TIME!  ISM  289  IBs  4sS4.52  T42B91985. JB8  HOIOH*.  SOUS  BTs  I. 


•i  > 


ff:.m 


0.20  e.S7  13.33  29.20  26. S7  33.33  40.20  46.67 

SECONDS 


6i  W*  **S  OF  IKTEOOTTEE  3UTE  M  =  IM  K  =  JJ675.  *JE  e  IM50.0  WNGE  =  336J.I 


DISCRETE  IMTEOWTIOH  OF  SLOTE  SEtSOTTS  HITH  TW  3F*C  3LDTE  01*0710* 


TTflrr  TJJCt  1964  289  19t  2)46  .  2  5428919*3. 118  HOIONfiL  SCTUD  BTs  l.« 


26. S7  53.23 

SECONDS 


srerr  tuci  ibm 


23 i  Si 42. SI  T42B82SS4.S16 


LOOL  SOLES  5Yi  1.5 


STOTT  TDCi  XS64  208  23i12i47.81  IU2BB23tS 


157 


*  *  * 


OJSOtETE  INTEOBTIW  OF  SLOPE  SEOOfTS  WJTN.TW 


DISCRETE  JHTEOWTIt*  OF  SLOPE  S£Q€XT$  KITH 


SECONDS 


oiscxrrt  inteowtion  c r  slope  seqcnts  mitm  t>€ 


SLOPE  OIKECTION 


DM  SUNS  OF  INTEGRATE!}  SLOPE  N  -  166  K 


START  TIME-  1964  296  22>46>43.  1  S429K217.S13 


REGIONR.  SCALE3  BTi  l.t 


V-'.*  «.»  K»  VJV*  V.W'4.'  ','T  ^.V^ryy.'  v;  %^y^  -J-. 


r«vT^r-'vv\'v,.v,’ 


lVi*j 


kVJj\ 


■>  ».»  ».» ;i 


INK  SUNS  OF  INTEGRATED  SLOPE  N 


DISCRETE  INTEGRATION  OF  SLOPE  SEOCNTS  WITH  TIC.SA*  SLIPS  DIRECTION 


ifl*  SLIPS  DIRECTION 


START  TIICi  1984  291  9.22i  S.8B  T 42918922.848  TELESEI9UC  SCALED  BT:  8.S 


■1  .u*Au 1 


SECONDS 


RAMC  SUNS  OF  INTEGRATED  SLOPE  N  =  190  t  -  11534.  KJE  =  18859.1  RANGE  r  S868.S 


DISCRETE  INTEGRATION  OF  SLOPE  SEGICNTS  NITN  TIC  SWC  SLOPE  DIRECTION 


s 


SECCND5  * 

Kurile  Islands  50. 3N  153. N  153.  5E  -  10/17/84  Origin  Time  :  0  9  :  1 1  :  0  4  .  F  -s 
Origin  Time:  09:  11  :04.8  -  Depth  1  34km+3  M.  :  5.0  ,s 

Residual  Error:  -  0.2  D  X 


1%  Pi  *■»  1 


IH«  SUNS  OF  INTEGRATED  SLOPE  K:|H  K  =  IIS89.  HUE  =  1B8S8. 0  RANGE  =  3786.8 


discrete  integration  cf  slope  sewents  with  the  smE  slope  direction 


STRUT  TIHEi  1984  291  16  t  2:  5.88  H429  U  582.838  REGIONAL  SCALED  BT :  1.8 


*'*v  * 

: i  v.t  f-  n 


j£vrr'j,|.* 


8.88  6.67  13.33  28.88  26.67  33.33  48.88  46.67  53.33  68.88 

SECONDS 


RANK  SUNS  OF  INTEGRRTEQ  SLOPE  N  s  108  K  r  12831.  HUE  =:  18858.8  (WOE  =  4522. 


DISCRETE  INTEGRATION  OF  SLOPE  SEGHENT5  WITH  Tt*  SRHE  SLOPE  DIRECTION 


START  TIHEi  1984  291  16:  1<S3.88  14291 1682. 8ST7  REGIONAL  SCALED  BT:  1.8 


8.88  6.67  13.33  28.88  26.67  33.33  48.88  46.67  S3. 33  68.88 

SECONDS 


Event  16  -  Regional  -  Probable  Lg  Arrival 


S 


1 


Kl 


DISCRETE  INTEGRATION  0 f  SUITE  SEGMENTS  WITH  THE 


SUITE  DIRECTION 


■me  9UQ  or  IKTEERRTD  SUFt  H  ■  Ui  K  ■  1**27.  HUE  ■ 


S 


a 

*1 


tf  , 


i 

w. 


I  j  |  r  j  j  I  I  I  1 

«.  02  6.67  13.33  22.80  26.67  33.33  48.88  46.67  53.33  68.00 

SECONDS 


„> 

$ 

$ 


§8 

.V, 


| - , - | - 1 - 1 - 1 - 1 - 1 - T 


(.00  6. 67 

13.33  22.00 

26 . 67  33.33 

SECONDS 

12.  ZZ  46.  67 

53.33 

Event  18 

-  Regional 

-  Probable  Lq 

Arr iva 1 

179 


«-  s."  «- 


~ 1 

60.00 


s?* 

v 

o 

V1 

vrl 


¥ 

‘3 


a 


mir»ir«vn irwirw Bimini ithmt"  vr* i«>»  >rxu-n* yiurun'.k *  u mt» wm *jt  nji rji ■jfjr 


SMB  OP  MTEEMTD  UPC  MU  D  11515.  MX  ■  1OT.9  MCE  ■  41SB.5 


j-i  DISCRETE  INTEGRATI&"  Of  SLOPE  SEGMENTS  MJTH  TX  SANE  SLOPE  DIRECTION 


up 


j-i  START  TlnEi  1994  291  IB;  19;  6.S1  142911519. BS9  I  REGIONAL  SCALED  Bit  1.5 


I - 1 - 1 - 1 - 1 - 1 - I - I - I - 1 

I.BB  6.67  13.73  29.50  26.  S7  33.73  40.00  46.67  53.33  60.00 

SECONDS 


W*«  SUKS  Of  INTEGRATED  SLOPE  N  =  100  t  s  12*74.  HUE  «  10*50.0  RANGE  •  3999.5 


I-,  DISCRETE  INTEGRATION  Of  SLOPE  SEGMENTS  NITh  THB  BANE  SLOPE  DIRECTION 


I-,  START  TIME:  1994  291  18<  1S< 57.  2  S429U819.B301  REGIONAL  SCALED  BT i  1.9 


I. BS  6.67 


13. 33  20.00  26.67  33.33  40.00  46.67  53.33  60.00 

SECONDS 


i|* 


SUMS  OF  INTtSMTED  SLOPE  M1M  (i  11934.  MX  *  IflSM  ■  M.J 


discrete  iHTEamuoN  of  slope  seocnts  kith  the  sh  slope  direction 


STOTT  Tift i  19*4  291  19iS9’33.S2  T42SU8S1.097  REGION*.  SOLED  BTt  S.S 


r*i  ■- jTi* 


6.67  13.33  20.00  26. 67  33.33  tB.BB  46.67  53.33  60  00 

SECONDS 


me  sum  or  imtawrEP  wan  n  ■  ua  t  ■  mat.  hue  »  lam.a  «m  ■  m.i 


DISCRETE  IWTECRflTICW  OF  SLOPE  SEGMENTS  WITH  TML  SMC  SLOPE  DIRECTION 


STOTT  T»C>  IK4  291  19<S0'39.B2  S4291 


HI 

lBSl.tlJ 


REGION*.  SOLED  BTt  1.9 


3.^  13.1 


22.  IS  26.67  33.33  40.00 

SECONDS 


46.67  53.33  69.00 


DISCRETE  INTEGRATION  Of  SLOPE  SEGMENTS  WITH  THEflSAHE  SLOPE  DIRECTION 


Dissert  iKTEcwrriw  of  slots  acoors  with 


H  V »]■’ L'i 1 1 L ^ '-  LVWW 


*h 


BMC  Site  OF  tHTEOBTIED  SUTE  N  =  1H  K  =  11783.  HUE  =  18BS0.9  RANGE  =  3958. S 
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■  .08  6.67  13.33  28.90  26.67  33.33  40.80  46.67  53.33  60.90 

SECONDS 
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BMC  SJHS  OF  INTEGRATED  SLOPE  N  =  190  K  =  12813.  KJE  =  100S0.0  RANGE  =  4112.5 


0ISCRETE  INTEGRATION  OF  SLOPE  SCOOTS  KITH  DC  SBC  SLOPE  DIRECTION 
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- 1 - 1 - 
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Event  26  -  Regional  (P  -  Lg)  =  48  seconds 
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agmmvwwiATO  w.'mv.K'sjw  m  minmmvm  jum',  wvrmivxvwm 


i - 1 - 1 - 1 - 1 - 1 - 1 - 1  i  i 

I.H  6.67  13.33  20.00  26.67  33.33  40.00  46.67 

SECONDS 


t-i 


HK  sue  OF  IHTEONTED  SLOPE  N  -  100  K  *  12098.  HUE  -  100S0.0  RflNCE  *  3998. 0 
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r - 1 — 

0.00  6.67 

Event  27 


— , - 1 - 1 - 1  i  i 

13.33  20.00  26.67  33.33  40.00  46.67 

SECONDS 

Regional  (Pn  -  Lg)  =  95  Seconds 
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g-,  mt  ue  v  untonTD  slope  h  ■  in  c  =  1212s.  ns  s  imss.i  ma  =  477s.  0 


j-,  discrete  iNTcawTiDH  of  sure  seoerra  kith  t»  srk  sure  direction 


I  I  III  lip  l| 

j-.  smrr  times  19*4  232  isisaiss.«r>  142921859.028  regional  soiled  ht.  r.s 


•.00  6.67  13.33  20.00  26.67  33.33  40.00  46.67  53.33 

SECONDS 


B-.  RK  SUMS  OF  INTEGRATED  SLDPC  N  =  100  K  =  12721.  HUE  *  10050.0  WNK  :  3923.5 


DISCRETE  INTEGRATION  OF  SLOPE  5EOIENTS  WITH  THE  SANE  SLOPE  DIRECTION 


START  TllCi  1904  292  10iSat23.01  S429218S8.056  REGIONAL  SOILED  BYi  1.0 
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i.s  J  •  | 


0.00  6.67  13.33  20.00  26.67  33.33  40.00  46.67 
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S3.  S3  60.00 
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RANC  sue  OF  INTEGRATES  SLOPE  N  =  IBS  K  =  12222.  HUE  a  10050.0  RANGE  a  S018.S 


DISCRETE  INTEGRATION  Of  3L0TE  SCGMNT3  HIT*  THE  S«C  SLOPE  DIRECTION 


H'lp  411  I  'I  »Y  \ 
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...  ,  .Is* 
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Event  29  -  Regional  (Pn  -  Lg)  =  37  seconds 
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ST  APT  TIME'  1964  292  Z3i  4i  2.81  T 42922384. 8SS  REGIONAL  SOLED  BTi  8.5 


—  seS& 
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8.87  13.33  28.88  28.67  33.33  48.88  46.87  S3. 33  68.88 

SECONDS 
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26.67  33.53 

SECONDS 


Event  30  -  Samoa  Islands  -  15. OS  17 1 . 1W-10 / 19 /8 4 
Origin  Time:  01:28:16.3  Depth:  1  km  +75  M^:  4.9 
Residual  Error:  -0.6  kJ 
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TllCl  IM  2*3  23:  3:31.  I  U29323I3.019  REG!  Oft  SOLES  ST:  S.5 
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DISCRETE  INTEGRATION  OF  SLOPE  SEGMENTS  WITH  T*  SAME  SLOPE  bMdlW 


START  TU€:  1904  293  23  :  3:40.51  S42S32303. 035  REGIONAL  SCALED  0Y:  1.0 


It 


r lip  ' 
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BW  SJKS  OF  INTEGRATED  SLOPE  N  =  !M  K  =  11547.  MX  -  19959.9  RANG E  s  S72S.8 
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I  6.67  13.33  29.99  “'IecONDS3'  ” 

Event  34  -  East  of  Severnoya  Zemlya  82. ON  114. 2E  10/21/84 
Origin  Time:  18:57:  55.6  Depth:  14km  +  4  M,  : 4 . 3 
Residual  Error:  -0.3  “ 


vw’.hw: 
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88  6.67  13.33  20.00  ^'SECOND?' 0'” 

Event  36  Komandorsky  Islands  Region  55. 7N  165. OE  10/23/84 
Origin  Time:  08  :04:46.1  Depth:  20  km  +3  M,  : 4 . 8 
Residual  Error:  0.3 
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v  wraw-V-Yi 


RANC  sue  OF  INTEGRATED  SLOPE  N  =  t«  K  =  12195.  M  =  IKM  «*C£  =  55*5. a 
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li 
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START  TINE  I  1984  297  17.52:44.83  S429717S3.413  LOCAL  SOLED  8Ti  1.8 


II  , 


13.33  28.ee 


26.67  33.33 
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48.ee  46.67  53.33  68.86 
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THE  CURIOUS  CASE  OF  THE  MISSING  EXPLOSION 


Eugene  Herrin 
Geophysical  Laboratory 
Southern  Methodist  University 

"Is  there  any  point  to  which  you  would  wish  to  draw  my 
attention?" 

"To  the  curious  incident  of  the  dog  in  the  night-time." 

"The  dog  did  nothing  in  the  night-time." 

“That  was  the  curious  incident,"  remarked  Sherlock  Holmes. 

(“Silver  Blaze",  A.  Conan  Doyle) 

This  narrative  could  properly  be  classified  as  historical 
fiction.  Most  of  what  is  reported  here  actually  happened,  but 
some  of  the  events  occurred  only  in  the  imagination  of  the 
writer.  We  begin  by  considering  the  capabilities  of  a  regional 
network  designed  to  monitor  an  area  of  thick  salt  deposits  in 
the  western  portion  of  the  Permian  Basin  of  Texas  and  New  Mexico. 
The  stations  in  the  network  are  at  Lajitas,  Texas;  Hobart, 
Oklahoma;  and  Winnemucca,  Nevada  (see  Figure  3).  Noise  levels  at 
the  three  stations  are  based  upon  actual  noise  observed  at  these 
sites  and  at  similar  sites.  The  minimum  background  noise  at 

Lajitas  is  the  lowest  ever  observed  in  the  frequency  band  of  5  to 

40  Hz.  This  minimum  noise  level  is  shown  in  Figure  1.  The 

minimum  levels  reported  for  NORSAR  (NORESS  site)  are  somewhat 

higher,  by  about  10  dB  or  more.  Measurements  at  Hobart  in  the 
frequency  band  1  to  4  Hz  show  noise  levels  similar  to  NORSAR. 
The  background  level  at  Winnemucca  based  on  early  measurements  is 
between  the  Lajitas  and  the  Hobart  values.  The  Lajitas  and 
Winnemucca  stations  have  state-of-the-art,  three-component,  short 
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period  instruments  in  50  ft.  boreholes.  At  Hobart,  a  15  element, 
short  period  array  of  surface  instruments  is  located  on  basement 
rock.  This  array  is  identical  to  the  NORESS  array  without  the 
outer  ring. 

Background  levels  shown  in  Figure  1  are  the  minimum  noise 
levels  under  ideal  conditions  which  exist  only  a  small  fraction 
of  the  time.  The  major  source  of  background  noise  is  the  effect 
of  wind  at  the  site  as  is  shown  in  Figure  2.  NORSAR  has  noise 
levels  similar  to  those  at  Lajitas  under  high-wind  conditions. 
We  expect  to  observe  the  same  effect  at  the  Hobart  array. 

Figure  3  shows  the  location  of  a  number  of  events  as  well  as 
the  locations  of  the  stations.  GNOME  was  a  3  kt  explosion  tamped 
in  salt  which  was  actually  observed  at  Lajitas,  Hobart  and 
Winnemucca.  SLEUTH  is  a  planned  3  kt  decoupled  shot  in  the 
Sal  ado  formation  in  the  same  general  area  as  the  GNOME  event. 
Figure  4  shows  the  strati graphi c  units  in  the  area.  The  Salado 
salt  provides  the  depth  and  thickness  needed  to  decouple  a  3  kt 
nuclear  explosion. 

The  southeastern  corner  of  New  Mexico  is  an  area  dotted  with 
potash  mines  as  shown  in  Figure  5.  Mining  potassium  bearing 
minerals  from  the  salt  and  anhydrite  units  is  accomplished  using 
the  room-and-pi 1 1 ar  method  with  the  separation  of  the  ore  being 
done  on  the  surface  near  the  working  shafts.  The  area  is  almost 
a  wasteland  covered  with  mounds  of  discarded  evaporites  and 
dessication  ponds.  Once  the  mining  has  proceeded  as  far  away 
from  the  working  shaft  as  is  practical,  the  pillars  are 
systematically  removed  allowing  the  mine  to  subside.  This 


collapse?  leads  to  obvious  surface?  effects  over  the  area  which  has 
been  mined.  After  this  procedure  is  completed,  a  new  shaft  is 
dug  and  the  operation  is  repeated.  Travellers  crossing  thi-s 
region  must  beware  of  dangerous  scarps  which  develop  in  the 
highways  above  collapsing  mines. 

In  the  midst  of  this  region  under  one  of  the  mined  areas  a 
cavity  with  a  40  meter  radius  has  been  constructed  at  a  depth  of 
1000  meters  in  the  Salado  salt.  The  salt  removed  in  this  process 
represents  only  a  small  addition  to  the  wastes  already  present  on 
the  surface.  A  3  kt  nuclear  device  is  placed  in  the  cavity  ready 
for  the  decoupled  test  which  has  been  code-named  SLEUTH. 

The  F'n  signal  levels  from  a  3  kt  shot  fully  tamped  in  salt 
(GNOME)  observed  at  the  three  stations  in  the  monitoring  network 
are  given  in  Table  1,  along  with  the  distances  to  the  stations. 
The  very  low  signal  level  at  Winnemucca  resulted  from  the  high 
attenuation  of  F'n  across  the  Basin  and  Range  province. 
Propagation  to  Lajitas  and  Hobart,  however,  is  as  expected  in  the 
mid-continent.  We  assume  a  decoupling  ratio  of  100  at  1  to  2  Hz 
for  SLEUTH  compared  to  GNOME,  thus  we  can  accurately  predict  the 
F'n  amplitude  levels  (1-2  Hz)  at  the  three  stations  for  the 
decoupled  shot. 

On  Tuesday,  25  June  19B5,  preparations  were  being  made  for  a 
4  kt  HE  shot  (MINOR  SCALE)  at  the  White  Sands  test  site  (see 
Figure  3).  A  cold  front  was  crossing  Colorado  and  Utah  at  that 
time  as  shown  in  Figure  6.  The  pattern  of  fronts  moving  west  to 
east  shown  in  this  map  is  typical  of  the  weather  pattern  in  this 
region  in  the  spring  and  early  summer  and  again  in  the  fall.  The 
frontal  movements  and  wind  patterns  are  highly  predictable.  On 


Wednesday , 


26  June  (Julian  day  177)  ,  the  cold  -front  had  passed 
into  New  Mexico,  but  had  not  yet  begun  to  a-f-fect  the  wind, 
patterns  at  Lajitas.  That  afternoon  an  earthquake  of  about 
magnitude  3  occurred  west  of  Amarillo,  Texas,  and  was  recorded  on 
the  high-frequency  (sample  rate  250  per  sec)  system  operating  at 
the  Lajitas  station.  The  distance  to  Lajitas  from  the  epicenter 
is  670  km.  Figure  8  shows  the  signal  and  the  si gnal -to-noi se 
spectrum  for  this  event.  There  is  good  si gnal -to-noi se-rati o  to 
frequencies  greater  than  15  He.  This  event  was  located  within 
the  16-element  seismic  network  operated  by  Stone  and  Webster 
Engineering  as  part  of  the  nuclear  waste  disposal  survey  in  the 
Texas  Panhandle;  therefore,  we  were  able  to  compute  an  accurate 
epicenter  far  the  event  using  the  network  records.  Digital  data 
were  available  for  a  station  48  km  from  the  epicenter  so  that  a 
good  displacement  spectrum  could  be  computed.  This  spectrum 
showed  a  clear  corner  at  6  He ,  a  constant  level  at  lower 
frequencies  and  a  roll-off  above  the  corner  frequency  of  60 
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dEt/decade  (f  ).  Using  this  spectrum  and  the  digital  record  at 
Lajitas,  we  were  able  to  produce  a  good  estimate  of  the  apparent 
Q  for  F'n  along  this  path.  The  value  of  Q  was  246  which  is 
consistent  with  Q(Fn)  reported  from  northeast  of  Moscow  along  a 
line  from  the  Volga  River  to  Vorkuta  (Yegorkin  and  Run, 
I z vest l ya ,  1978,  Vol .  14,  No.  4,  262-269). 

On  Thursday,  27  June  (Day  178)  1985,  the  cold  front  had 
m aved  through  Oklahoma  and  much  of  Texas  (Figure  9).  Winds  at 
Hobart  and  Lajitas  were  20  to  30  mph  from  the  north.  That 
morning  MINOR  SCALE  was  fired  at  White  Sands,  and  thirty  seconds 


later  the  decoupled  nuclear  shot,  SLEUTH,  was  -fired.  Figure  10 
shows  the  signal  from  MINOR  SCALE.  The  amplitude  of  F'n  from  this 
event  was  about  4  times  as  large  as  from  the  Amarillo  earthquake 
the  previous  day,  as  would  be  expected  based  on  the  yield  of 
MINOR  SCALE.  The  wind  at  Lajitas  was  20-30  mph  at  the  time  of 
the  event  (Figure  10),  whereas  conditions  had  been  nearly  calm 
during  the  recording  of  the  Amarillo  earthquake  on  the  previous 
day.  Even  though  the  vertical  instrument  was  located  at  a  depth 
of  100  meters,  wind-noise  was  a  major  problem  on  day  178.  This 
effect  is  clearly  seen  in  Figure  11,  where  the  si gnal -to-noi se 
spectrum  falls  to  zero  dB  at  6  Hz  for  MINOR  SCALE.  This  result 
can  be  contrasted  with  the  effective  bandwidth  of  16  Hz  seen  for 
the  more  distant,  significantly  smaller  event  recorded  on  the 
previous  day  when  the  wind  was  nearly  calm. 

SLEUTH  was  not  detected  at  Lajitas.  Pn  was  below  the  noise 
level  and  the  Lg  wave  train  was  swamped  by  the  Lg  signal  from 
MINOR  SCALE.  Figure  12  shows  the  predicted  displacement  spectra 
for  GNOME  and  SLEUTH  based  on  the  corner  frequencies  for  3  kt 
tamped  and  decoupled  events  given  by  Archambeau  and  the  Q-value 
obtained  from  the  Amarillo  earthquake.  The  displacement  spectra 
of  the  background  noise  on  days  177  and  178  are  also  shown.  From 
this  figure  we  see  that  if  SLEUTH  had  been  fired  on  day  177,  Pn 
would  have  been  detected  at  Lajitas,  in  agreement  with 
Archambaau's  predictions.  By  picking  the  right  time  to  fire 
SLEUTH,  based  on  the  weather  patterns  and  the  known  time  of  MINOR 
SCALE,  the  most  sensitive  high-frequency  station  in  the  network 
was  made  incapable  of  detecting  the  event.  At  Hobart,  Oklahoma, 
the  wind  was  15  to  20  knots  (around  25  mph)  from  the  north 
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(Figure  9).  The  surface  array  there  could  be  expected  to  produce 
si gnal -to-noi se  improvement  over  a  single  surface  instrument  of 
about  3.5  at  1  to  2  Hz  and  about  2  at  10  Hz.  The  background 
noise;  however,  could  be  expected  to  be  higher  by  at  least  a 
factor  of  5  because  of  the  high  winds.  Thus  the  Hobart  array 
would  fail  to  detect  Pn  from  SLEUTH.  Again,  Lg  would  be  lost  in 
the  coda  of  MINOR  SCALE,  and  could  not  be  pulled  out  by  array 
processing  because  of  the  similarity  of  azimuths  for  the  two 
events  relative  to  the  array.  Because  of  poor  propagation  across 
the  Basin  and  Range  Province,  F'n  from  SLEUTH  could  not  be 
detected  at  Winnemucca  no  matter  what  the  noise  conditions  were 
at  that  station.  Thus  we  see  that  an  excellent  reginal  networl 
designed  to  monitor  an  area  with  salt  deposits  failed  to  detect  a 
3  kt  decoupled  shot. 


"The  dag  did  nothing  in  the  night-time." 


The  failure  of  the  network  occurred  because  the  evader  could 
pick  the  most  advantageous  time  to  fire  the  clandestine  test. 
SLEUTH  is  not  one-of-a-kind.  Predictable  weather  conditions 
similar  to  those  on  27  June  1985,  occur  in  West  Texas  and 
Oklahoma  several  times  each  year.  Shots  of  1/2  kt  of  HE  or 
larger  are  not  uncommon  at  White  Sands  Test  Site.  Cavity  mining 
could  go  on  year  after  year  in  southeastern  New  Mexico  completely 


5ked  by  the  potash  mining  in  the  area. 


By  waiting  for  the 
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same  sequence  o-f  events  that  occurred  on  27  June,  several 
decoupled  nuclear  tests  a  year  could  be  carried  out  with 
virtually  no  chance  o-f  detection. 

Two  -factors  are  required  for  a  successful  clandestine  test 
of  the  kind  described  in  this  paper.  The  first  is  a  degradation 
of  detection  capability  at  the  critical  stations  because  of  hiqh 
wind.  The  second  is  a  legitimate  HE  test  under  the  control  of 
the  evader  which  can  be  used  to  mask  the  decoupled  shot.  The 
second  factor  must  be  regulated  by  legal  means.  Control  of  the 
first  factor  depends  upon  our  ability  to  protect  instruments  from 
the  effect  of  wind-induced  seismic  noise.  Until  these  problems 
are  salved,  a  clever  evader  can  pick  the  right  time  to  fire  a 
decoupled  shot  with  little  risk  of  detection. 
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Figure  1.  Minimum  Seismic  Noise  Levels  at  Lajitas  and  Norsar 
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Figure  3.  Map  showing  locations  of  seismic  stations  and  events. 
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Figure  4.  Stratigraphic  section  of  GNOME/SLEUTH  site. 


Figure  5.  Map  showing  potash  mining  activity 
in  the  GNOME/SLEUTH  vicinity.  Dots 
indicate  location  of  potash  mines. 
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Figure  8.  Amarillo  earthquake,  High-frequency  recording 
(250  samples  per  second)  from  330  ft.  Z 
component  at  Lajitas  and  signal-to-noisc  spectrum 
for  this  record. 


Figure  9.  Weather  map  for  Thursday,  27  June  1985 
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Figure  10.  Minor  Scale:  Record  from  330  ft.  borehole 
instrument  and  plot  of  wind  speed  during 
P-wave  arrival  time. 
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